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The Critical and Two-Phase Flow of Steam’ 


The results of a digital computer and analytic study of the critical flow of a compressi- 
ble fluid are presented in this paper. The expanding flow of a fluid in a single-phase 
region as well as the expansion of a fluid to a two-phase region is considered and de- 
scribed by analytic expressions relating choking velocity, critical pressure ratio, and 
flow per unit area characteristics. A comparison ts made of the analytic results which 
assume a constant value of the isentropic expansion exponent, with the digital computer 


results using the actual properties of steam. All analyses assume the fluid to be in 
thermodynamic equilibrium. 

A skeleton Mollier diagram is presented for steam showing the exponent in the wet and 
superheated regions. The choking velocity is presented in plot form as a function of the 
inlet conditions as well as state point conditions; critical pressure ratio 1s presented as a 
function of inlet conditions. The critical flow per unit area ts presented in the form of a 
factor K plotted versus inlet conditions; this factor K when multiplied by inlet 


Introduction 


be critical flow of compressible fluids has been con- 
sidered and discussed in the past by many authorities [1, 2, 3]?; 
the topic is of intense interest to those concerned with fluid 
handling apparatus that may operate at or near the critical flow 
condition. The additionally imposed problem, of the fluid ex- 
panding to a region where two phases of the fluid exist under 
equilibrium conditions, must also be grappled with. For ex- 
ample, this latter aspect is one that is of some importance to the 
designers of turbo-machinery and water cooled nuclear reactors. 
The following presentation will be concerned with the one-dimen- 
sional, isentropic flow of a compressible fluid at the point of 
minimum area in a converging or a converging-diverging passage. 


Discussion and Computer Study 


A flow per unit area maximization study was made on an IBM 
704 digital computer in which a general problem was set up to 
accept the arbitrary expansion line inlet conditions of pressure 
and enthalpy. An exit pressure was then estimated at which 
flow per unit area would maximize following an isentropic ex- 
pansion? (i.e., reach choking conditions). An iterative scheme 


1 Awarded Melville Prize Medal for Original Work, awarded for the 
best original paper on a mechanical engineering subject by a member 
of ASME presented during 1959. 

2? Numbers in brackets designate References at end of paper. 

3 Flow per unit area G equals the product of p times V, where p 
is the fluid specific weight at the expansion line end point and V is the 
corresponding velocity, 

V = [2Jg(hr — Netatic) ]'/2 

Contributed by the Power Division and presented at the Annual 
Meeting, Atlantic City, N. J., November 29—-December 4, 1959, of 
THe AMERICAN Soctety OF MECHANICAL ENGINEERS. Manuscript 
received at ASME Headquarters, August 26, 1959. Paper No. 59— 
A-223. 


pressure produces the desired value of critical flow. 


is then used to obtain the final exhaust pressure; velocity, exhaust 
specific volume, exhaust pressure, and flow per unit area values as 
well as inlet conditions are printed out. This choking velocity in 
the superheated steam region corresponds to the velocity of 
propagation of a small disturbance, commonly defined acoustic 
velocity. As the wet region is a region of co-existence of steam 
and water droplets which are assumed to be homogeneously sus- 
pended, the definition of the velocity of small disturbance is 
nebulous although the computed velocity at maximum G may be 
still considered that of a theoretical choking value; i.e., if one 
considered the wet region to be in equilibrium and its pressure- 
density characteristics those of a single phase substance, the 
computed velocity at maximum G will be that of small disturbance 
propagation with respect to those equilibrium relationships be- 
tween pressure and density. The computer study was made us- 
ing the properties of steam as presented in the ASME Trans. 
[4]. 

An analytic study was then undertaken to substantiate the 
computer results; these expressions are presented in the following 
section. 

The analytic study produced two analyses: The first con- 
siders the isentropic expansion exponent to be a constant, the 
second attempts to consider the variation in the exponent. As 
the computer results inherently include the variation in the 
exponent (being close approximations to the true steam proper- 
ties) they provide the base to which all numerical comparisons 
are made. 

For convenience in the analytic study a complete and con- 
sistent knowledge of the isentropic expansion exponent was neces- 
sary. A skeleton Mollier diagram has been included as Fig. 1 
showing the exponent in the superheated and wet regions. As 
can be seen, the exponent’s characteristics in the two domains have 
completely different trends. Values for the superheated steam 
region were obtained from Fig. 8 in the steam tables [5] and wet 
region values from Fig. 4 of a paper by W. H. Wiebe [6]. 


Nomenclature 
c = velocity of propagation of a small p = specific weight, #/ft* D = superheated steam region (dry) 
disturbance, ft/sec A = area, ft? S = saturated vapor line conditions 

g = acceleration of gravity, ft/sec? G = weight flow per unit area, #/ft?-sec T = total conditions 
h = enthalpy, Btu/# V = velocity, ft/sec W = wet steam region 
p = pressure, #/ft? W = weight flow, #/sec ( )s = along an isentrope 
s = entropy, Btu/#-deg R Subscripts 
t = temperature, deg F 2 = exhaust conditions in two phase Special Terms 

= s i 3 y2 
hr = total enthalpy, hr = hetatie + 5— 
Y = isentropic expansion exponent 0 = expansion line inlet conditions 2/9 
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ENTROPY-BTU/LB-°F 
Fig. 1 Mollier diagram showing the isentropic expansion exponent of 
p 


steam: pv’ = constant; y = — {| 
Pp Op 


Pk ts of choking velocities are presented in Figs. 2 through 5. 
Pigs. 2 and 3 present choking velocity as a function of inlet pres- 
sure and enthalpy, the former figure being concerned with the 
wet acu transition regions and the latter figure, the superheated 
region. Figs. 4 and 5 present choking velocity as a function of 
state point pressure and enthalpy, the former figure presenting 
the we* region and the latter figure, the superheated region. The 
critie | pressure ratio Pexhaust/Pinlet for both the superheated 
and vet regions are presented in Fig. 6. The numerically higher 
ratios tucrespond to the wet region. The critical flow factor K is 
ineluced in Fig. 7 as a function of inlet conditions. 


Analytic Study 


The square of the velocity of propagation of a small disturbance 
in a compressible fluid is expressed as 


This -relocity corresponds to the velocity attained by a gas at 
choking conditions (maximization of its flow per unit area charac- 
teristi¢s, a continuous function) in a converging passage following 
an isentropic process (adiabatic and reversible), as well as the 
maximum velocity at the end of a long constant area pipe, in 
which the flow process is adiabatic and irreversible (Fanno-type 


- flow). 


Flo v per unit area is defined as 


G=-— = pV (2) 


then 
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INLET ENTHALPY - BTU/LB 
Fig. 2 Choking velocity in the wet and transition regi as af 
of inlet conditions 


and at its maximum, dG/dp = 0, therefore 


d 1V 
(3) 


dp V dp 


Considering the process to be isentropic (as mentioned above), 
equation (3) may be rewritten as 


__ p fav 


As the fluid velocity at maximized flow per unit area is expressed 
by equation (1), equation (4) is rewritten as 


2 d 2 
Op/, 9g \op/, 
Equation (5) is the general equation used for determining density 
change with respect to pressure at constant entropy, and at chok- 
ing conditions. _ 

After determining appropriate expressions for (0V/dp), and 
p, the resulting expression for (0p/dp), is then introduced into 
equation (1) yielding the desired velocity. This method is fol- 
lowed for the two considerations of single and double phase ex- 
pansion, with constant and variable isentropic expansion ex- 
ponents. 
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Fig.3 Acoustic velocity in the superheated steam region as a function of 
inlet conditions 
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STATE POINT ENTHALPY-BTU/LB. 
Fig. 4 Choking velocity in the wet region as a function of state point STATE POINT ENTHALPY-BTU/LB. 
conditions Fig. 5 Acoustic velocity in the superheated region as a function of state 


point conditions 


CRITICAL PRESSURE RATIO FOLLOWING AN ISENTROPIC EXPANSION 
OBTAINED FROM COMPUTER STUDY 


Single Phase Region 
L END POINT | Constant Isentropic Expansion Exponent. The velocity expression 
IN SUPER HEATED REGION ——-—~ determined in Appendix A is the well-known expression 

Variable isentropic Expansion Exponent. The velocity is expressed 
(see Appendix B) as a function of inlet conditions as well as 
ry TS (d7/dp),, where is the exponent of eitber the wet or dry region. 
Bey 1 1 The complete analytic expression is determined by the substitu- 
2 1300 tion of the following expression for (O0p/Op), in equation (1). 
1100 +—+—+ 4 1 
— EXPANSION LINE END PONT (=) (=) 
| 1000 REGION vw \y +1 — 1) op/, 
54 55 57 58 59 60 61 

(26) 
CRITICAL PRESSURE RATIO- £-EXHAUST op/, 


P INLET 
Fig. 6 Critical pressure ratio as a function of inlet conditions 


Transactions of the ASME 


148 / aprit 1961 


|. 
: 
beg 
ais, 


ow 


TLE 
TT 


| 


THEORETICAL CRITICAL FLOW OF STEAM 
PER SQUARE INCH OF NOZZLE AREA 
VS. INLET PRESSURE AND ENTHALPY 


FOR 100 % FLOW COEFFICIENT 
BASED ON KEENAN & KEYES 1936 STEAM TABLES 
FLOW = ABSOLUTE INLET PRESSURE x FACTOR “K"(LB/HR-IN®) 


| 


+4 


INLET ENTHALPY, BTU/LB. 
Fig.7 The critical flow factor K as a function of inlet conditions 


Comparison By equation (14); constant y c = 1817.1 ft/sec 
Superheated Inlet and Exhaust By equations (26) and (1); variable y c = 1819.6 ft/sec 
po = 50 psia By IBM 704 (flow maximization) c = 1820.0 ft/sec 
ho = 1300 Btu/# 
vo 11.7166 ft?/# Wet Inlet and Exhaust 
Yo 1.301 = 10 psia 
27) 1100 Btu/# 
= —0.2 ft? 
(2 36.73 ft?/# 


By equation (14); constant y 1751.8 ft/sec 1.132 
By equations (26) and (1); variable y 1754.2 ft/sec —r 
By IBM 704 (flow maximization) 1757.1 ft/sec sie aa li 


Superheated Inlet and Exhaust By equation (14); constant y 1344.2 ft/see 
Po 10 psia By equations (26) and (1); variable y 1342.7 ft/sec 
t = 600 deg F By IBM 704 (flow maximization) 1342.3 ft/sec 
ho 1335.1 Btu/ * 
63.03 ft*/# 
ms = 1.3009 Double Phase Region 


(2%) —1.6 X 10-5 ft2/# Constant Isentropic Expansion Exponent. The velocity expression 
op/, derived in Appendix A is as follows: 


2gporo YD Yo(Y" — ‘a Yo(¥w — 
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t. The velocity is deter- 


Variable isentropic Exp i Exp 
mined as in the previous section on single phase with variable 
exponent. Two exponents are encountered, yp and yy, the 


superheated and wet region exponents, respectively. 


(°°) + 2CD (2) 


op 2(A + B) Uo op 


OY w 210) (22) 
2CE {| 2DE {| — — D? — 


oYw \? 
+m 


where the terms A, B, C, D, and E are defined in Appendix B. 


Comparison 
Superhected Inlet and Wet Exhaust 
po = 10psia Computed by 
ho = 1,160 Btu/# IBM 704 
= 40.586 ft'/# = 5.8203 psia 
So = 1.81258 Btu/#-deg R = 63.0036 ft?/# 
p, = 7.3614 psia 
Yo = 1.3195 
Yw = 1.1328 
(Oyp/dp), = —0.555 X 10> ft?/¢ 
(Oyw/dp), = 1.045 x ft?/4 
By equation (22); constant y’s P2 = 5.8045 psia 
By equation (24); constant y’s c = 1387.0 ft/sec 
By equations (27) and (1); variable 7's c = 1384.0 ft/sec 
By IBM 704 (flow maximization) c = 1384.96 ft/sec 


Transition Region 


From the results of the IBM 704 flow maximization study it 
was seen that severe changes in choking velocity were encountered 
as inlet enthalpy entered a narrow range of values in the super- 
heated steam region. The range of these severe velocity changes 
is called the transition region, where closely spaced points were 
computed for precise definition of this region. The span of this 
region was found to lie between 60 and 70 degrees superheat as 
plotted on a Mollier diagram (covering a pressure range from one 
inch Hg absolute to 2,000 psia). In analyzing the transition re- 
gion data the end point in every case was for all practical pur- 
poses on the saturated vapor line. 

Constant isentropic Expansion Exponent. The choking velocity will 
then be that generated by the isentropic drop from the inlet point 
in the transition region to the saturated vapor line (see Diagram 


3), i.e., 


2apes ( 
1 Po 


Variable isentropic Expansion Exponent. As the wet region is not 
actually entered by this process, the expression for (0p/dp), for 
the double phase expansion with variable isentropic expansion 
exponent may be modified to appear as follows, 


op /, 


2B v6 
The choking 


(6) 


(7) 


where B, C, and D have been previously defined. 
velocity is found by equation (1). 
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Comparison 
po = 10 psia Computed by IBM 704 
ho = 1175 Btu/# Pp2 = 5.6637 psia 
Vo = 42.548 ft?/# = 65.4208 ft3/# 
8) = 1.8340 Btu/#-deg R 
Pp, = 5.662 psia 
Yp = 1.3185 
= —0.555 X 10> ft*/# 
op /, 
By equation (6); constant y c = 1446.7 ft/sec 
By equations (7) and (1); variable y c = 1446.74 ft/sec 
By IBM 704 (flow maximization) c = 1446.97 ft/sec 


Results and Conclusions 


As has been explained, the computer study maximized flow per 
unit area, following an isentropic expansion from some arbitrary 
inlet state point. Fig. 2 shows a portion of these results, choking 
velocity versus inlet enthalpy with inlet pressure parameters. 
The sudden changes occurring in the previously defined transition 
region are quite apparent, appearing as nearly vertical lines on 
the chosen grid. Figs. 4 and 5 show no such transition region, as 
these plots present state point characteristics [see equation (1)] 
and are not dependent on any process relationships. 

The plot of critical pressure ratio versus inlet conditions in- 
dicates two regions. That region of lower numerical magnitude 
contains the expansion line end point in the superheated steam 
region, while the other contains its end point in the wet region. 
Connecting the various pressure parameters would be a line (not 
shown) for the transition region; the means of obtaining the locus 
of this line is explained in the transition region section. 

The analytic expressions presented closely check the values of 
velocity obtained from the digital computer study. Although 
the variable exponent expressions are somewhat more accurate 
than the fixed exponent expressions, results within 0.5 per cent 
may be obtained using the latter equations. 
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APPENDIX A 


Constant isentropic Expansion Exponent 


Single Phase Region, The consideration of an isentropic expan- 
sion in a homogeneous domain with no property discontinuities 
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constitutes a single phase expansion. The energy equation (as- 
suming the process adiabatic, workless, at constant elevation, 
and the initial velocity to be zero, i.e., dhy = 0) may be written 


as 


See Diagram 1. 


Diagram 1 Single phase expansion diagram 


The incorporation of the relationship pu? = constant (where 


7 is the isentropic expansion exponent) in the inviscid motion 
_ (dp Vav 
equation {| — + —— = 
p 


0) yields an expression for the square 
g 
of the velocity 


V2 = 29 ( 


upon differentiation, dV /dp is obtained. 
1 
g vo (” 
dV L2 Po Pp 


dp 
( 
Po 


The well-known critical pressure ratio (i.e., that ratio at which G 
is a Maximum) is 


which when substituted into equation (9) yields 


1 


dV gw (yt1\]4 f(y +1 
(11) 
dp 2 Po | 2 
Also, as pv? is constant along an isentrope, p at the critical pres- 


sure may be expressed as 


where p = 1/v. Upon the substitution of equations (11) and 
(12) into equation (5), (0p/Op), is determined to be 


(2) (2 + ') 
Op], Y 


The choking velocity is then determined by equation (1) to be 


(12) 


(13) 
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c= (- Povo (14) 
which is the commonly encountered equation as a function of inlet 
conditions. The same results may be obtained by the substitu- 
tion of equation (10) into equation (8). 

Double Phase Region. The consideration of an isentropic ex- 
pansion with property discontinuities (i.e., from one homogeneous 
domain to a second homogeneous domain, where separation is by 
a common boundary line) constitutes a two phase expansion. 
The energy equation may be written as 


ho =h,+ 


See Diagram 2. 


ve 


Diagram 2 Double phase expansion diagram 


we 


Diagram 3 Transition region expansion diagram 


The square of the velocity V, may be expressed similarly 
to equation (8) where we now use p, and yp (defining p, to 
be the intersection of the isentrope and the saturated vapor line, 
and ¥Y p to be the isentropic expansion exponent in the dry region) 


Also, h, — may be expressed using p,, and (where 
is the isentropic expansion exponent in the wet region) as 


Yw Ps P2 
— hy {— (2 (17 
) J ) 


Since 
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The density p: is expressed as 
1 1 


Po Ps 


Upon differentiation of equation (18) the following expression is 
obtained 


At this point, some expression for the critical pressure ratio is 
needed [say, something similar to equation (10)]._ We proceed to 
write the expression for G as follows . 


Ps 


By taking the derivative of equation (21) with respect to 
(p2/p,) and setting this equal to zero, an expression for the de- 
sired critical pressure ratio is obtained [(po/p,) is known and as- 
sumed to be a constant]. 

Then 


= 0 


and the critical pressure ratio is 


Yor! Yw-1 
2 YD 27 
( + | _ tre (22) 
+t Ps Yw(¥w + 1)(¥o — 1) 


Upon the substitution of equation (22) into equation (20) and 
then substituting this result and equation (19) into equation (5), 
the final expression for (0p/dp), is as follows 


ra) 1 

2povo E _ — Yo(Yw uf 
(yw + 1) : (yp — 1) 


P2 


Po 


The choking velocity is then determined by equation (1) to be 


(yw + 1) 


Po (Yp 1) (Yp 1) 
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(21) 
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The same result for c may be obtained by the direct substitution 
of equation (22) into equation (18). 


APPENDIX B 
Variable Isentropic Expansion Exponent 


As part of the analytic study, a somewhat irregular investiga- 
tion was undertaken yielding results that are of surprisingly 
greater accuracy than those obtained from the constant exponent 
analysis. The procedure followed was based on the assumption 
that the expressions for velocity (equations (8) and (18)) were 
such that all of the quantities expressed by symbols (i.e., p, p2, Y:, 
Yp, and Yw) are indeed variables, i.e., all are subject to differentia- 
tion. 

Single Phase Region. Following the same general procedure as 
that in the body of the report, the derivative of the velocity ex- 
pression, equation (8), is taken with respect to pressure, and 
with the use of equation (10), we obtain 


(y- + 1) 2 op 
(25) 


By substituting equations (25) and (12) into equation (5), 
(0p/dp), is determined to be 


1) 


2 


2 


ty 
oy ) 2 
(26) 
(3), 
The choking velocity is again determined by equation (1). In 
the comparison, values of (Ov/Odp), are obtained from Fig. 1. 
Double Phase Region. Again the derivative of the velocity ex- 
pression (this time equation (18)) is taken with respect to pres- 


sure, (0p/dp), is determined from equation (5) using equations 
(19) and (22) to be as follows: 


2 2 


\ po Ps 


2(A + B) 


ovo) 


where 


- [1-( 
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The choking velocity is again determined by equation (1). 
In the comparisons, values of (Oyp/Op), and (OYw/Op), are ob- 
tained from Fig. 1. 


APPENDIX C 


It became necessary in the study to obtain values of choking 
velocity at state points near the saturation line and at various 
other locations for check purposes and for the precise definition 
of the trend of curves presented in Figs. 2 through 5. The 
derivative (0p/dp), was computed by the following equations 
obtained from NACA TN2214 [7]. Consider Diagram 4 where 
the values of p at equally spaced values of p have been com- 
puted along an isentrope. 


oP 


Diagram 4 Numerical differentiation diagram 


The slopes are then for the corresponding points: 


Oth point 
(2) = [—2.0833p0 + 4.09: — 3.0p2 + 1.339; — 0.25p,] 
op /, Ap 
2nd point (mid-point) 
op 


= = [0.0833 — 0.669: + 0.669; — 0.083sp,] 


4th point 
[0.25p0 — 1.33p1 + 3.0p2 — 4.0p3 + 2.0833p,] 
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Error terms corresponding to the foregoing numerical ap- 
proximations may be found in the cited report. 


DISCUSSION 
Gordon J. Van Wylen* 


Two questions can be raised regarding this interesting and 
stimulating paper. All of the analytic study is based on the as- 
sumption that the steam, both in the superheat and two-phase 
regions, is an ideal gas with constant specific heat. This is a 
reasonable assumption in the superheat region at low pressures, 
but the validity of this assumption can be questioned both at high 
pressures in the superheat regions and in the two-phase region. 
Perhaps it is significant that the comparisons between the IBM 
704 program and the analytical results which are cited are all at 
low pressures. It is noted that the agreement between the IBM 
704 program and the analytic study is less at 50 psia than at 10 
psia. How do these compare at pressures above 1000 psia? 

The second question involves the basic validity of taking an 
expression which was derived assuming constant specific heat 
(and therefore constant y) and differentiating it to find the effect 
of variable y. This raises the question of the general applicability 


4 Professor and Chairman, Department of Mechanical Engineering, 
University of Michigan, Ann Arbor, Mich. Mem. ASME. 
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of this approach and whether or not the rather good agreement 
between the IBM 704 program and the analytical studies is some- 
what coincidental. 


Author’s Closure 


I would like to express my appreciation to Professor Van 
Wylen for his interest in this paper and for his several appropriate 
questions. With regard to the fluids’ adherence to the perfect 
gas relationships during the processes under consideration, I 
have found excellent agreement both in the superheated and wet 
steam regions. At the higher pressure levels in the superheated 
region, a comparison of several points at pressures up to 1000 psia 
has shown close agreement with the several plotted relationships. 
To quote one point, consider it to be at low entropy with an 
inlet pressure of 1000 psia and at a temperature of 740 F, equation 
(14) indicates that c = 1826.38 ft/sec, the computer results 
yield 1828.80 ft/sec. 

The general concern over the differentiation of an expression 
with respect to a term considered constant is certainly justified; 
acknowledgment of this fact may be found in the Introduction 
and Appendix B. However, the point remains as to whether 
this procedure may yield correct trends when examined for 
changes in this pseudo-constant. The affirmative view is sup- 
ported by the results of the paper’s comparison points. 
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GEORGE W. FORMAN’ 
NEIL W. KELLY’ 


Introduction 


te MECHANICAL DRAFT COOLING TOWER has become 
an important equipment component in chemical plants, petroleum 
refineries, electric power plants, and large air-conditioning in- 
stallations. Lack of natural water in adequate quantities for 
cooling purposes often dictates the use of evaporative water cool- 
ing in order that available water may be recirculated within the 
plant. Numerous publications have dealt with the thermal per- 
formance of cooling towers, but little has been presented concern- 
ing the performance of the large axial-flow fans which have come 
into common cooling tower use. Misconceptions exist concerning 
the effects of differences in fan diameter, number of blades, blade 
shape, fan rpm, the design of air discharge cylinders, stacks, and 
so forth. This paper is presented to clarify some of these miscon- 
ceptions and to show a method of predicting the performance of 
cooling tower fans which requires a minimum of dimensional and 
descriptive data. 


Performance Requirements 


The performance requirement of a given cooling tower fan ap- 
plication depends upon the thermal duty of the cooling tower 
served. The volume of air required may range from a few 
thousand to nearly a million cubic feet per minute, and the static 
pressure against which the fan must operate will ordinarily be 
less than */y-in. water gage. The fan operates in a horizontal 
plane at the top of the cooling tower in a warm, moisture-satu- 
rated air stream. The cooling tower designer’s problem is to select 
from a fan series the fan diameter, number of blades, and rpm 
which require the minimum horsepower consistent with available 
space, available reduction gearing, and first cost of the over-all 
cooling tower installation. 


1 Assistant Professor of Mechanical Engineering, University of 
Kansas, Lawrence, Kansas. Mem. ASME. 

2 Manager of Manufacturing, J. F. Pritchard and Company of 
California, Kansas City, Mo. 
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Nore: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those of 
the Society. Manuscript received at ASME Headquarters, March 16, 
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| Cooling Tower Fan Performance 


Fundamental Considerations 


Theoretical analysis of the performance of a fan does not yield 
precise information, and tests are necessary in order to obtain 
exact performance data. However, considerable insight may be 
gained through a simple and straight forward analysis by con- 
sidering the fan blade as an airfoil. The forces acting on the 
blades (and hence the forces exerted by the blades on the air 
stream) may be estimated by considering an element of the blade 
as shown in Fig. 1. 

An element of blade of length dr at radius r will operate 
essentially as does a section of an airplane wing in flight. 

Referring to Fig. 2, the forces acting on a wing may be calcu- 
lated from 


L = p/2 X VA (1) 


p/2 X (2) 
where 
L = lift force, acting perpendicular to the direction of 
relative wind, pounds 
D = drag force, acting parallel to the direction of 
relative wind, pounds 


p = mass density of the air, pounds X sec?/ft* 


Nomenclature 


lift force on an airfoil acting per- Cp 
pendicular to direction of rela- 
tive wind, lb 

drag force on an airfoil acting 
parallel to direction of relative 
wind, lb 

= mass density of air, lb sec?/ft* 

axial air velocity, ft/sec 

velocity of relative wind, ft/sec 

airfoil lift coefficient, dimension- 
less 


less 
= area, sq ft 


radius, ft 
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airfoil drag coefficient, dimension- 


airfoil angle of attack, deg 


angle between direction of relative 
wind and plane of fan rotation, 
deg 
angle whose tangent is D/L 
= fan efficiency 


resultant force on airfoil, lb = total pressure rise, lb/sq ft 
= component of force in the plane of 
fan rotation, lb = 
component of force perpendicular 
to plane of fan rotation, lb 
fan speed, revolutions/sec 


= velocity pressure, lb/sq ft 

fan solidity, the decimal fraction 
of fan circumference at radius r, 
occupied by blades 

total blade width at radius r, ft 
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A = wing plan area, ft? 
C,andCp = coefficients deranding upon the angle of attack a 
and the airfoil shape 
a = angle of approach of wind relative to the airfoil 
V = velocity of relative wind, ft/sec 


The efficiency of the blade element at radius r will be deter- 
mined in the following analysis. Referring again to Fig. 1, the 
lift and drag forces for the blade element may be calculated, and 
the resultant air force dR on the element found as the vector sum 
of the lift and drag forces. By resolving the resultant force into 
components parallel to the axis of rotation of the fan and in the 
plane of the fan, the useful output or working force dP and the 
torque component dQ may be determined. From these, the ef- 
ficiency may be found. From Fig. 1 


tan @ = V/2mrN, (3) 
where N = speed, revolutions per second. 
Also, from Fig. 2 


p /2 
pV?C,A/2 


tane = D/L = = Cp/C, (4) 


The useful power output is the product of the force component 
parallel to the shaft axis and the velocity of the sir 
Output = VdP (5) 


The power input is the product of the torque component dQ 
and the velocity at which it moves 


Input = 24rNdQ 
The efficiency of the blade element will be 
_ Output dP 


Efficiency = = Input ~ = dQ tang (7) 
From Fig. 1, 
dP = dR cos (¢ + €) (8) 
dQ = dR sin (6 + €) (9) 
Thus by substitution 


In this expression for efficiency, the angle ¢ depends on the air- 
foil lift and drag: If the airfoil is ‘‘good” ¢€ will be small. The 
angle ¢ depends on the ratio of the air velocity through the fan to 
the tangential velocity of the blade element. Since both € and @ 
usually vary from hub to tip of the fan blade, an integration is 
necessary to obtain the over-all efficiency of the fan accurately. 
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Fig. 3 Total fan efficiency versus air velocity to tip speed ratio 


Fortunately, the efficiency of the fan as determined by an in- 
tegration from blade root to tip is approximated by the efficiency 
of the blade element three quarters of the distance from the center 
of rotation to the blade tip, i.e., at 75 per cent radius. Fig. 3 
shows a plot of fan total efficiency calculated from Equation (10) 
and relates efficiency to the ratio of air velocity through the fan to 
fan tip speed, which is proportional to tan ¢. The family of 
curves represents airfoils of different lift/drag ratio. 

It is apparent from Fig. 3 that total efficiencies in the order of 
ninety per cent are possible, provided the airfoil lift/drag ratio is 
about 20:1 and the helix angle @ (as indicated by the ratio air 
velocity through the fan divided by fan tip speed) is within cer- 
tain limits. Total efficiency, however, is not a criterion of the 
relative power required by different fans in handling the same air 
flow against the same static head. A more meaningful indication 
is given by “static efficiency,”’ defined by 


Nt = static efficiency = Not 
static pressure 
static pressure + velocity pressure 


(11) 


Note that the static efficiency of the fan is always less than 
its total efficiency, but that the static efficiency approaches the 
total efficiency as the discharge velocity approaches zero. Thus 
the fan is charged with the discarded kinetic energy, and the fan 
with the highest static efficiency will require the minimum input 
horsepower to deliver a given air volume against a given static 
pressure. 

Fig. 4 shows the relationship between efficiency and air velocity 
through the fan for various static pressures. This plot is based on 
a constant fan tip speed of 10,000 feet per minute and an airfoil 
lift/drag of 15. Current practice is herein exemplified: Com- 
mon tip speed values range from 8000 to 13,000 ft per minute, 
and an L/D of 15 represents a “good” airfoil application. 

For fans operating at tip speeds other than 10,000 fpm, a cor- 
rection factor has been calculated (Fig. 5) to be applied to the 
efficiency indicated by Fig. 4. Similarly, a correction for an air- 
foil lift/drag ratio other than 15:1 must be applied. The airfoil 
correction factor is shown in Fig. 6. The corrected efficiency 
may be found from ’ 


Corrected efficiency = Efficiency, Fig. 4 X F, X Fy (12) 


A brief study of the information thus far presented indicates 
that fan static efficiencies of about 50 to 60 per cent are possible 
for fans operating against static heads of approximately 0.5 in. of 
water. Higher efficiencies would require lower than common tip 
speeds, or airfoils that are considerably better than the ones in use 
today. 
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Fig.4 Total and static efficiency at constant tip speed of 10,000 fpm and 
an airfoil lift/drag ratio of 15 
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Fig. 5 Correction factor for Fig. 4 for tip speed other than 10,000 fpm 


Figs. 4, 5, and 6 can be used for predicting the horsepower re- 
quired for any cooling tower fan application if the operating L/D 
of the fan is known or can be determined. 


Number of Blades and Blade Width 


The lift/drag (L/D) ratio at which the airfoil operates depends 
on the physical shape of the airfoil and the angle of attack. 
Many systematic wind-tunnel studies of airfoil contours have 
been conducted, and the findings are readily available in reports 
of The National Advisory Committee for Aeronautics. Fig. 7 
shows a plot of the lift coefficient C, and the lift/drag ratio as 
functions of the angle of attack for a Clark Y airfoil. This par- 
ticular airfoil has found wide use in axial flow fans and aircraft 
propellers and can be taken as typical of a good airfoil. 

By taking the characteristics of the Clark Y airfoil as typical 
of common cooling tower fan designs, the value of L/D can be de- 
termined for given operating conditions for a given fan “solidity.” 
The procedure is as follows: 

Referring to Fig. 1, the force dP is equal to the total air load on 
the element of fan disk area at radius r, i.e., is equal to the total 
pressure differential across the fan disk times the area of the ele- 
ment dr. In equation form 


dP = p,2mr dr (13) 


where p, is the total pressure rise, pounds/foot.? Also, 
dP = dR cos (o + €) 
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Fig. 6 Correction factor for Fig. 4 for airfoil lift/drag ratio other than 15 
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Fig.7 Airfoil characteristics for Clark Y airfoil 


Combining Equations (13) and (14) 
dR = p,2mrdr/cos (@ + €) (15) 
But 
dR = dL/cos € = C,WdrpV?/2 cos € (16) 


where W = the total blade width in feet at radius 7, i.e., the prod- 
uct of the number of blades and the width of each blade. 
Setting the two expressions for dR equal to each other 


C,WdrpV?2/2 cos € = 2mrdrp,/cos (¢ + €) (17) 


p.2/pV? = C,W cos + €)/2mr cos € (18) 


Let W/2rr = S = solidity, i.e., the decimal fraction of the 
fan circumference at radius r occupied by blades, then 


p,2/pV? = CS cos (@ + €)/cos € 
V = V/sind 


(19) 


From Fig. 1, (20) 


Substituting 


p.2/pV? = C,S cos (@ + €)/cos € sin? d (21) 


since pV2/2 = p,, we can write 
p./pS = Cz cos + €)/cos € sin? d 


In Equation (22), C, and € are associated with the performance 
of the airfoil. This equation may be used to determine the L/D 
at which a given fan must operate to produce a particular volume 
of flow against a known static head. (Performance at the */, radius 
is taken as the “average” performance of the entire blade.) Fig. 8 
shows a plot of Equation (22) based on the Clark Y airfoil. Fig. 9 


(22) 
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. 8 Chart for determining airfoil lift/drag ratio required for given 
requirement 


presents a graph from which the value of p,/p, can be easily read 
for known values of air velocity and static pressure. It is ap- 
parent from the last derived equation that a required fan solid- 
ity may be provided by either a few wide biades or several nar- 
row blades. The method of fabrication of the blades controls the 
design choice of number of blades: The economics of cast blades 
calls for a relatively large number of narrow blades; with fabri- 
cated construction fewer and wider blades are used. r 

By using Fig. 8 in conjunction with Figs. 4, 5, and 6, the horse- 
power required for any fan can be predicted if the following in- 
formation is known: 


1 Air volume, cfm 

Static pressure 

Fan diameter 

Fan rpm 

Number of blades 

6 Blade width at */, radius 


orm to 


To illustrate the method of performance calculation consider 
the following example: 


A six-bladed, 19-ft diameter fan with a 63-in. diameter hub is 
to handle 372,000 cfm against a static pressure of 0.400 in. water 
(actual values), and the air has a density of 0.0705 lb/ft?. The 
blades are 18*/;,in. wide at the */, radius, and the fan is to operate 
at 175 rpm. What horsepower is required? (Note: This set of 
numbers has been chosen for illustration, since under these con- 
ditions a field test gave a motor output of 50.8 BHP.) 

Solution: Fan tip speed = 7 X 19 X 175 = 10,450 ft/min 
Air velocity through fan = 372,000 XK 4/7 X 

{192 —'(63/12)?] = 1420 ft/min 
Air velocity/tip speed = 1420/10,450 = 0.136 
The static pressure that would exist in ‘standard air’’ would be 
Standard static pressure = 0.400 X 0.075/0.0705 = 0.425 in. 
From Fig. 4, the uncorrected static efficiency = 0.56. 
From Fig. 5, the tip speed correction factor = F, = 0.995 

To find the airfoil correction factor F, from Fig. 6, we must 
first find the L/D at which the airfoil will operate. Assume the 
fan to have a Clark Y airfoil. (The actual characteristics will not 
be greatly different from the Clark Y for a fan of reasonable de- 
sign.) 

6 X 18.19/12 


= ss ————— 658 
rX19X*/ 


From Fig. 9, for V = 1420 and p, = 0.425, p,/p, = 4.40 
Thus, p,/p,S = 4.40/0.203 = 21.7 
From Fig. 8, L/D = 19 
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Fig. 9 Chart for determining total pressure/velocity pressure ratio 


From Fig. 6, F, = 1.05 

The static efficiency = 0.56 0.995 X 1.05 = 0.585 

BHP = air horsepower/static efficiency 

cfm X static pressure, lb/ft? 
33,000 

372,000 X 0.400 X 62.5/12 
33,000 

Horsepower required at the fan shaft = 23.5/0.585 = 40.2 

If the drive mechanism has a 90 per cent over-all efficiency, the 
motor horsepower required = 40.2/0.90 = 44.7, which is within 
12 per cent of the 50.8 experimental value. 

In this example, the agreement of the calculated value of the 
power required and the experimental value is not as precise as de- 
sired, but about as good as normally obtained. Considering the 
difficulties inherent in field measurements, it is believed that the 
method of calculation here presented is of considerable practical 
value. 


Air horsepower (static) = 


= 23.5 


Air horsepower = 


Fan Performance by Tests 


No entirely satisfactory instrument has been developed for use 
in measuring velocities in the air stream of cooling tower fans. 
Conventional vane anemometers have been used, but the ac- 
curacy of calibration as used in the warm, humid air stream that 
contains water droplets has not been established. The pitot 
tube is perhaps the most reliable instrument, but accurate meas- 
urements here would require the use of a precision manometer 
which is not well-suited to field test conditions. The air flow 
on the discharge side (which is the only conveniently accessible 
location) is not strictly parallel to the fan axis, and correction 
must be made for the rotational velocity of the discharge air 
stream. In general, a volume measurement in full scale testing is 
subject to an error in the order of plus or minus ten per cent. 
The power input to the fan is impossible to determine, and is 
generally taken as the power input to the motor driver less as- 
sumed values for motor and transmission efficiencies. Tests 
under field conditions are often influenced by prevailing winds, 
and a complete range of static pressures cannot be conveniently 
established. The static pressure is determined by measurements 
taken in the plenum chamber below the fan, and as such does not 
include the friction loss due to flow entering and leaving the fan 
cylinder. (The method of calculation presented in this paper is 
based on a static pressure difference across the fan disk, and as 
such must include all static pressure losses throughout the entire 
system, including the fan cylinder. This difference in static pres- 
sures probably accounts for part of the discrepancy between 
calculated and measured horsepower in the example.) 
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Fan manufacturers’ performance data are usually obtained 
either through accurately conducted tests of scale models, by 
theoretical calculation, or through the use of limited full-scale 
field testing. Well-made models under laboratory conditions pro- 
vide accurate data which are finding increased use among fan pro- 
ducers, although the basis of published ratings is often vague. 


v/s 


FAN BHP 


% 100 200 300 400 500 


CFM STANDARD AIR, THOUSANDS 
Fig. 10 Performance characteristics of 19-ft diameter 6-bladed fan with 


63-in. diameter hub based upon 3-ft diameter scale model tested in 
standard flared wood cooling tower fan cylinder 


Test data (or other basis of rating) may be presented in many 
ways, one of which is shown in Fig. 10. Note that this plot is for 
a specific fan as tested in a fan throat of a particular configuration, 
and for a tip speed of 10,000 feet per minute. The static pressure 
drop due to air flow through the fan cylinder is here charged 
against the fan, and as such the fan and its cylinder are rated as a 
unit. Fig. 10 was obtained from a 36-in. diameter model of the 
19-ft fan discussed in the example and indicates that 58 horse- 
power is required to handle 372,000 cfm at 0.425 in. water stand- 
ard static pressure. Thus the power required at an actual air 
density of 0.0705 Ib/cu ft, 


HP = 58 X 0.0705/0.075 = 54.5 


and if the power transmission efficiency is 90 per cent, the motor 
horsepower required is 54.5/0.90 = 60.5. Agreement with the 
full scale field test value of 50.8 is not entirely satisfactory. As in 
this example, the efficiency of models is usually somewhat lower 
than full scale fans, due to the difficulty of maintaining tip clear- 
ances that are accurately to scale, and the fact that the model 
operates at a lower Reynolds number. 


Fan Cylinders and Stacks 

As previously mentioned, the theoretical method of per- 
formance estimation set forth in this paper is based on the static 
pressure difference across the fan disk. In using the method of 
calculation, it is necessary to estimate the static pressure loss in 
the fan cylinder in order to obtain best accuracy. The fan cylin- 
der static pressure loss may be taken as some percentage of the 
velocity pressure at the fan disk. The factor may be as low as 2 
per cent, if the cylinder is of well designed proportions with a 
flared inlet and may be as high as 200 per cent for a sharp edge 
orifice mounting. The usual commercial design of cooling tower 
fan cylinder will have a static pressure loss in the order of one 
quarter of a velocity head at the cylinder. 

Ordinarily, extended stacks as are sometimes used to assist in 
discharging the humid air at a distance of 20 or 30 feet above the 
fan will have a negligible effect on the fan performance. If the 
stack is flared outward some kinetic energy recovery will result, 
which may be of practical significance under conditions of very 
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low static pressure (say in the order of !/, in. of water or less). 
The method of fan performance estimation presented here may be 
used with reasonable accuracy, and the static pressure reduced by 
the amount of velocity head recovery. It is doubtful that more 
than about one half of the kinetic head can be recovered eco- 


nomically. 


Summary 


A simple method of predicting the performance of cooling tower 
fans has been presented and compared by example with results of 
both model and full scale tests. The method may be used to ad- 
vantage by purchasers of cooling tower fans to determine the 
validity of performance data presented by prospective suppliers 
of cooling tower fans. Although the calculated horsepower may 
be in error by as much as 20 per cent, the predicted air flow will 
be within approximately six per cent (the cube root of 1.2), which 
is within field test accuracy. Manufacturers of cooling towers 
and cooling tower fans have, over the years, indicated accuracies 
of fan performance and horsepower prediction that are incon- 
sistent within the state of knowledge of the art. 


DISCUSSION 
J. D. Holmberg’ 


This paper is an excellent presentation of a method of calculat- 
ing the maximum capability of a cooling tower fan. It clearly 
shows the development of the simple blade element theory which 
has generally been used in cooling tower fan design. This is the 
first time the theory has been presented in this form, which 
greatly simplifies calculations through judicious choice and use of 
curves. The authors are to be commended for the effectiveness 
of their approach. 

It should be kept in mind, however, that the calculated 
performance represents a performance “‘ceiling.’’ The actual per- 
formance is always something less, partly because of losses which 
oceur in the application of the fan and partially because the 
theory is based on several assumptions which are not precisely 
true. These basic assumptions are discussed, not as a criticism 
of the authors but to emphasize the fact that the calculations do 
establish a maximum performance which the fan can, at best, only 
approach. 

This theory assumes that the discharge air flow is truly axial; 
actually it always has a rotational component. Airfoil drag co- 
efficients are higher and maximum lift coefficients are lower at the 
lower Reynolds numbers corresponding to actual fan operation 
than at the Reynolds numbers for which the NACA has published 
airfoil characteristics. The method of calculation also assumes 
that the performance of the fan as a whole is represented by that 
of the blade elements at 75 per cent radius. This is most nearly 
true at design blade pitch and at design static pressure, and it 
becomes less accurate as an approximation as either of these is 
changed. 

Calculations from the theory are based upon operation of the 
fan under ideal conditions. In the actual application, losses will 
occur because of back flow at the hub and blade tips, nonuni- 
form velocity distribution of air approaching the fan, and/or 
angularity of inlet air flow. These losses are influenced by the 
design and arrangement of other cooling tower components, 
particularly the fan cylinder. In applying a cooling tower fan, 
the ultimate goal is to obtain fan performance very nearly equal 
to that calculated by the method of the paper. How closely 
that goal is approached can be determined only by test. 

The rating of a cooling tower is no more accurate than the 
rating of the fan. The fan designer, guided by limitations learned 
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from experience, can successfully use the theory in the fan design. 
One fan design can be compared with other fan designs by thor- 
ough tests of an accurate scale model. Performance ratings, 
however, are acceptable only after having been confirmed through 
field tests of the full-scale fan as applied on the cooling tower. 
The paper points out the difficulty of conducting accurate field 
tests. It must be agreed that it is a difficult, painstaking task. 
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It is only after numerous tests and gradual refinement of testing 
methods that we have reached a point of being able to obtain 
accurate field test results. Volume measurements consistently 
agree with the fan rating within three per cent. Differences 
greater than that would call for a thorough investigation as to the 
cause. It has been through investigations of past differences 
that correct procedures for applying, rating, and testing cooling 
tower fans have been developed. 
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GEZA VERMES 


Engineer, Steam Turbine Department, 
Allis-Chalmers Manufacturing Company, 
Milwaukee, Wis. Mem. ASME 


A Fluid Mechanics Approach to the 
Labyrinth Seal Leakage Problem 


The paper describes investigations of labyrinth seals carried out recently; derives new 
theoretical and semitheoretical formulas for computation of the leakage which agree 
within 5 per cent with the tests for three different types of seals; off-design performance 
of the seals is treated theoretically and experimentally. 


I, was more than 50 years ago that Prof. H. M. 
Martin first presented a formula in the British periodical Engi- 
neering [1]! for the computation of the rate of flow through a 
labyrinth packing, but the problem still keeps haunting mechani- 
cal engineers. An impatient discusser [2] referred to the problem 
some time ago as one among those which fascinate the mathe- 
matically minded—a polite way of calling it a useless play 
of the mind. Nevertheless, research work continued—a book- 
let has even been written on the subject and, surprisingly 
enough, published in war-torn Germany in 1943 [3]. The interest 
has not subsided even up to the present time [4]. 

Two reasons might be cited to explain this maintained activity. 
The rapidly growing size of turbomachinery makes the investiga- 
tion of relatively small losses worth while. In addition, the re- 
cent possibility of leakage of contaminating nature requires a 
more accurate calculation procedure. These factors constitute 
the first reason for the work done in this field. 

Another impetus for activity might be derived from the very 
nature of the gland process. It can be seen from Fig. 1 that the 
processes in an ideal labyrinth are of a rather orderly nature; 
acceleration of the fluid in Zone 1 with subsequent dissipation 


1 Numbers in brackets designate References at end of paper. 
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of the kinetic energy generated in Zone 2 (the very low velocity 
in Zone 3 does not contribute to the energy processes). These 
two processes are repeated throughout the labyrinth, and they 
were computable in their simplest form in Professor Martin’s 
time. As flow processes became more amenable to calculation, the 
possibility of understanding more complicated gland processes 
became apparent. This growing possibility of understanding 
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Fig. 1 Ideal and straight seal 
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= virtual number of sealing 


points 


= Reynolds number 


distance of sealing points, in. 

pressure, psia 

pressure drop across a regu- 
lar sealing point, psf 

total pressure drop across an 
elbow, psf 

gas constant, 53.3 ft/deg R 

width of exhausting jet, in. 

velocity in a constriction, fps 


= weight flow, lb/sec 


axial distance, in. 

radial distance, in. 

virtual seal factor 

tip to clearance ratio of seal- 
ing point, see Fig. 1 
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gland factor based on virtual 
number of seals 
= specific weight, pef 
see Fig. 1 
see equation (21) 
pressure ratio 
density, lb sec?/ft* 
constant 
= elbow loss factor 


® 


fr QD Aa or 


Subscripts 
= center line 
geometrical 
long 
atX =S 
short 
in the X-direction, see Fig. 1 
in the Y-direction, see Fig. 1 


Q 
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gland processes might be the second reason for the sustained 
activity concerning the subject. 

Most of the results to be presented in this paper are based on 
some findings in fluid mechanics which permit the application of 
Martin’s classical approach to more complicated seal configura- 
tions. Throughout the following discussions gland velocities will 
be assumed subsonic unless expressly stated otherwise. The flow 
from one constriction to another will be considered incompressi- 


bie. 


Martin’s Formula 


Martin’s well-known formula gives the mass flow for an ideal 


gland as follows 
1 ( Pn y 
Po Po 


W = 5.68 KA, 


(See Nomenclature for notations.) A gland is considered ‘‘ideal”’ 
when the reconversion of kinetic energy into heat is virtually 
complete; i.e., velocities in Zone 3, Fig. 1, are negligible. Martin 
assumes in his derivation that the process as a whole is isothermal 
(To is approximately equal to T,,) and velocities are everywhere 
below sonic velocity. The formula was rechecked without as- 
suming an isothermal over-all process (as Martin did), and the 
reheat factor too has been taken into account. The rather in- 
volved formula obtained has been compared with that of Martin 
for low-pressure air with pressure ratios ranging from critical to 
0.95 and for various numbers of sealing points from 3 to 34. It 
was found that by taking 5.76 instead of 5.68 in equation (1), 
Martin’s formula gave results within 3 per cent of those obtained 
with the more rigorous derivation. The dependence of the 
formula on k = C,/C, is such that with decreasing k the agree- 
ment between the more exact formula and that of Martin gets 
better. Writing 


= f(N, 


Ot 0.2 0.4 06 081.0 


é 


Fig. 2 Clearance factor of annular orifices 


2.0 40 6.08010 
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and taking 5.76 as the constant, the formula 


W = 5.76 KA, 8 (1a) 


should give sufficiently accurate results when K is known. As 
it is pointed out in Martin’s original work, K is the clearance fac- 
tor of an annular orifice. The results of an excellent paper [5] on 
annular orifices are redrawn in Fig.2.. As A = f(L/C, Nre) and as 
Nre is a function of the unknown mass flow, the necessary first 
approximation can be made with K = 0.67. Equation (2) is 
shown in Figs. 7 and 8 of reference [6]. 

As is well known the difficulties with Martin’s formula start 
when (a) Zone 3 cannot develop because the pitch is too short, or 
(b) inflow and outflow in some (or all) of the constrictions is of a 
more complicated nature than in Fig. 1, or (c) both difficulties 
prevail. In other words, Martin’s formula is satisfactory as far as 
thermodynamics are concerned; the problems to be solved are of 
a fluid mechanics nature. The possibility (a) is the case of a 
straight seal, Fig. 1, when the pitch is relatively short, (b) can 
be found, e.g., in a “stepped seal,’ Fig. 3, and (c) is the case of a 
“combination” seal, Fig. 4. 


Straight Seal 


It has been shown [4] that in case of a relatively short-pitched 
straight seal, Martin’s formula is still applicable provided the 
“residual energy’’ is known. The expression “residual’’ refers to 


D= 10" 
| 
Fey 
L=0.015" 
psf 


Fig. 3 Stepped seal 


Fig. 4 Combination seal 
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the kinetic energy of the jet emanating from a constriction that is 
not converted into heat when acceleration into the next constric- 
tion begins. 

Based on [4] equation (1) can be written 


A, Po 
W = 516K (3) 


where @ is the relative amount of kinetic energy at the beginning 
of Zone | (i.e., at the end of Zone 2, Zone 3 nonexisting). An ade- 
quate appraisal of a@ can be made with the aid of the boundary- 
layer theory (see Appendix 1); it was found that 


(4) 


Stepped Seal 

The seal shown in Fig. 3 has two types of constrictions; i.e.5 
the one denoted by A is similar to those of the straight seal, 
but the sealing point denoted by B is different. Whereas the A 
constrictions should have a constriction factor K that can be ob- 
tained from Fig. 2, the B constrictions will have a different 
clearance factor AK’. Denoting the pressure drop through the 
regular A type constriction with Ap, and through the B type 
sealing point with Ap, A’ can be defined as 


dp _(AK\_ 

Ap,  \4,K’) ~ \K’ 

because the pressure drop is proportional to the square of the 

velocity and velocity is inversely proportional to the flow area. 
Writing 

(6) 


where 


8.0 


5.0 60 7.0 


Fig. 5 Loss factors for stepped seals 
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either K 
K, = f(n) = (7) 
or K’ 


the gland can be considered as having constrictions of variable 
area. This problem was treated in reference [7]; the resulting 
equation (in terms of this paper) is 


1/3 


Po 


W = const 
PoAo 


n=0 
Assuming AjAy ~ (A,K)?, (Ao/Ay) © 1 and const = 5.76 as in 
(1a) equation (8) can be written 


(9) 


W = 5.76 KA, 


(10) 


B’ = f(N', 


(11) 


N=n 
N’ = f (K/K,,)?dn 
n=0 


Comparing equations (1a) and (9) it seems appropriate to call NV’ 
“virtual number of seals.’”?” When N’ is known the stepped seal 
yields to the same equation as the ideal seal. 

Equation (11) can be written approximately as 


N+4 


w~ 


2 

as half of the sealing points are regular and the other half have 
the same (K/K’)? factor, the points being identical among them- 
selves. Evaluating equation (12) in Appendix 2 the type B seal- 
ing points will be considered as a sequence of two elbows with loss 
factors £ and (C/t)? both of which are functions of x/C, Fig. 5. 
Using Fig. 5 to obtain values of & and (C/t)? equation (12) can 
be computed by writing 


C\2 
(13) 


Combination Seal 

The concept “virtual number of seals” introduced with the 
stepped seal can be used for solving the problem of the combina- 
tion seal, Fig. 4. It can be seen from the sketch that the longer 
seals can be treated as the type B sealing points of the previous 
section, the upstream shorter points as “‘ideal’’ seals, and the 
downstream shorter ones as sealing points of a straight seal. 

Equation (9) can be used again for obtaining the mass flow, but 
equations (12) and (13) have to be rewritten. It can be seen 
from Fig. 4 that the combination seal can be considered as con- 
sisting of N* groups of seals each composed of three different 
sealing points. (K/K’)? for the long points can be taken from 
Fig. 5; it is 1.0 for the upstream short points; and—based on 
equation (3)—it is 
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a 
= 
where 
= 
and 
= 
x 
C 
xX 
| 50 \ Als (14) 
40 = fo (4) (14a) 
20 30 


K 
(=) = K —\ =l-a (15) 
(1 — a)” 


for the downstream shorter seals. Therefore 
N’ = +&+ +1+(1- a) | = N*y” (16) 


and 
y"=y'+1+(1-a) (17) 


Using equations (15), (16), and (17) it can be written 
8" = f(N", (18) 


and thus equation (9) can be used for this type seal too. It should 
be noted that with the type of seal in Fig. 4 


N* = = (19) 


Off-Design Performance of Stepped Seals 


It is a question of practical importance in turbomachinery work 
that leakage through labyrinth seals can be predicted when the 
relative position of shaft and housing changes; e.g., different 
parts of the machine undergo a different thermal expansion. 

In case of the stepped seal (Fig. 3), no basic difference will 
occur in computing N’ as y’ consists of two members anyway, 
equation (13), and—according to Appendix 2—they are inde- 
pendent of each other. For a nonsymmetric position of the sealing 
point, (1 + &) and (C/t)? have to be taken from Fig. 5 at different 
values: X, for the upstream coefficient (1 + &) and X, for the 
downstream coefficient (C/t)?. Equations (9) and (10) through 
(14a) can be used as before. 

Off-design performance of the combination seal can be computed 
too. The long sealing point contributes the member y’ in equation 
(17); it should be computed with X, and X, as in case of the 


‘stepped seal. The short sealing points may change their 


geometry. In Fig. 6 (relative downstream shift of seal) C for 
the downstream short sealing point changes to C’ therefore 


y” = y' + (C/C’)%1 — a’) +1 (17a) 
where 
9 
a’ = = (4a) 
+ 7.23 


In case of an upstream shift of the seal (Fig. 7), the upstream 
short sealing point changes its geometry; therefore instead of 
1.0, (C/C’)? has to be written as its contribution in equation (17) 


+ (<) - (176) 
and a’ is taken from equation (4a) and Fig. 6. Thus obtaining 
y”, 8” can be computed with equations (16) and (18). 


Sonic Velocity in the Gland 


It is shown in [2] that the pressure distribution in the gland 
(p,/Pe-1) is a monotonously decreasing function of n. There- 
fore when the over-all pressure-ratic 7 is lowered, the last (Nth) 
constriction only will develop sonic velocity; and equation (1a) 
still will be applicable to the N — 1 subsonic sealing points 
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W =5.76K (1b) 


It can be shown that the relative magnitude of (py-i/po)? being 
small no substantial error will be caused by taking (py/po) in- 
stead of (py-1/po). With this assumption 


(Beritical)y (20) 


A dashed line in Fig. 6 of [8] indicates the region where equation 
(20) should be applied. The approximate error is, e.g., 3 per 
cent for 10 seals. 


Comparison With Tests 

The previously derived formulas were checked with air in static 
tests. An outline of the test fixture can be seen in Fig. 8; the 
fixture was preceded by a standard measuring duct connected to 
a supply of compressed air. The maximum initial pressure, po, 
used was 70 psia. Air flow was measured with standard ASME 
orifices; the apparatus discharged into open air. 
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Fig. 7 Upstream off-design position of combination seal 
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Fig. 8 Test fixture 
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The three types of glands that were tested and the significant 
dimensions of the seals are shown in Figs. 1, 3, and 4. The tests 
were conducted with radial clearances ranging from 0.004 to 
0.060 in. and in case of the stepped seal up to 0.100 in. Theory 
and measurements are compared in 


Fig. 9 for the straight seal. 
Fig. 10 for the stepped seal. 
Fig. 11 for the combination seal. 


It can be seen that, the smallest clearances excepted (C < 
0.015 in.), agreement between theory and measurement is within 
5 per cent; in case of the combination seal, it is 4 per cent for all 
the clearances. It should be mentioned here that the experi- 
ments with the small radial clearances for the straight and 
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stepped seals were conducted at the beginning of the program 
when the procedure for establishing the exact radial clearance was 
not quite adequate. Therefore, the accuracy of the theory might 
be considered as 5 per cent. 

An important assumption in the derivation of Martin’s formula 
is that the number of sealing points be very large. It is, there- 
fore, of interest to investigate whether the correlation with tests 
prevails in case of a labyrinth with only a few sealing points. A 
straight seal (C = 0.040 in.) with generally the same dimensions 
as in Fig. 1 but with every other sealing point removed has been 
tested, thus giving information on the performance of a labyrinth 
with five sealing points. Results, as shown in Fig. 12, maintain 
the same correlation as those with the larger number of seals. 

Some of the off-design performance measurements for stepped 
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Fig. 12 Correlation for straight seal, C = 0.040 in., 5 sealing points 
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seals are shown in Fig. 13 (C = 0.004 in.), Fig. 14 (C = 0.021 in.), 
and Fig. 15 (C = 0.100 in.); for combination seals in Figs. 16 (C 
= 0.005 in.), 17 (C = 0.020 in.), and 18 (C = 0.050 in.). These 
correlations give less favorable results than those with the seals 
- in symmetrical position, the general trend being the same; large 
clearances give results closer to the calculated values. The maxi- 
mum deviation is 15 per cent; in most cases it is 7-8 per cent. 


Conclusions 

1 Relatively simple considerations of a fluid-mechanics nature 
allow the extension of the validity of Martin’s formula for ideal 
seals to seal configurations in industrial use and give results within 
5 per cent of those obtained from measurements using air. 

2 The formulas describe the off-design performance of the 
seals mostly within 8 per cent but at least within 15 per cent of 
the test results. 
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Fig. 13 Stepped seal off-design performance, C = 0.004 in. 
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3 The same type of approach might be useful for seals of not 
too different configurations from those tested. 
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APPENDIX 1 
Residual Energy in a Flat Jet 


In order to compute the residual-energy factor @ in a seal 
chamber, the following assumptions will be made: 


(a) D (Fig. 1) is very large in comparison with C and E. 

(b) The amount of energy lost through friction along the quasi- 
flat smooth surface d-d is negligible in comparison to the mixing 
losses with the stagnant fluid in the chamber. 

(c) The point Xo where @ has to be computed will be the “be- 
ginning” of the inflow process into the next constriction. Zone 3 
will be considered very short and the flow pattern of the expand- 
ing jet supposed to prevail up to Xo. It will be assumed that the 
laws of the potential inflow process start to exert a remote effect 
on the flow at 2 aC where the kinetic energy is 2 per cent of the 
velocity head in the constriction in case of inflow from an in- 
finite container with fluid at rest. This condition does not prevail 
in the actual case and serves only to establish the limit of the 
potential flow in Zone 1. Concerning the choice of the value of 
2 per cent it has to be considered that 1 per cent change in 
measured values was the limit of the instrumentation used and 3 
per cent was the desired accuracy of the method. 

Formulas were derived on basis of [9] for the calculation of a 
and it was found to be 1.26; therefore 2 aC = 2.52 C. 

(d) Based on assumptions (a), (0), (c), the flow in a gland- 
chamber should be represented by one half of a flat symmetrical 
jet originating from an infinitesimally narrow slot (in Fig. 1, the 
slot is at X = 0). The flow of a jet of this type has been described 
by various researchers; the following considerations are based on 
{10 and 11]. 


As kinetic energy is proportional to the square of the velocity, 
the amount of residual kinetic energy at Xo can be defined if the 
average of velocity squares in the X-direction, ux?, is known there. 
(It can be shown that uy? is very small in comparison with wat) 
Defining a new parameter 


€= x’ (21) 


where o = constant: 


ux = (ucr)s () (1 — tanh? e) 


X (22) 


The foregoing equation is based on equation 23.45, p. 499, in 
reference [10] and (uc1)s is the velocity at the center line of the 
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| N | 
0.7 
ae 


two-dimensional jet at some characteristic axial location X = S. 
The constriction was chosen as the characteristic location. a@ can 
be written as 


(23) 
us? 
where 
Y=+o 
uy*dm 
um? (24) 
f dm 
Y=0 


by definition, dm being the mass flow passing through a cross 
section of unit depth and dY width in the elementary two-dimen- 
sional jet of p density: 


dm = puy 1dY (25) 


Combining equations (22), (24), and (25) gives 


Y=+o 
(1 — tanh? €)*dY 
= (uct*)s (26) 
(1 — tanh? e)dY 


from equation (21) (taking X = Xo) 
dY = de 
o 


and from formulas 679.90 and 679.40 in reference [12] 


2 
(28) 
Xo 15 
As equation (28) refers to the square of the maximum velocity, 
(uei?)s, whereas equation (23) to the average of the squares of 
the axial velocity, a relationship has to be established between 
(uct)s? and us?. According to reference [10] loc. cit., the 
kinematic momentum per unit length is constant in the jet and it 
can be written 


M= uy*dY (29) 
for the half jet. 
As 
2 S 
M = — (uct)s* — (30) 
3 


(loc. cit.), which in turn has to equal 


M = (31) 
therefore 
(wate = (32) 
3:48 
From (23), (28), and (32) 
a=0.80¢0 (33) 


In determining the ratio C/X» the angle 6 is essential, Fig. 1. It 
can be seen from equation (22) for uy and later from the formula 
for uy that the flow in a two-dimensional jet is similar along ¢ = 
constant lines. These lines originate from the point X = Y = 0, 
Fig. 1. It can be assumed that uy = 0 along e*. 

As 
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(reference [10] loc. cit.) it can be written 


1 
2(coth «* — tanh 


(35) 


e* 


A graphical solution of the foregoing equation gave e* = 1.083. 
o is constant deriving from properties of turbulent flow; its 
value is o = 10.64 [11].? It can be seen from Fig. 1 that 


=e* 


C 
= 


Therefore from equation (21) 


and 


= = 9.75 (38) 


X> can be computed from statement (c) and equation (36) accord- 
ing to Fig. 1 


X=S+P-—2aC-L 


(39) 


X S P-L 
C —2a (40) 


Taking equations (38) and (40) and substituting in (33) gives 


(41) 
at + 9.75 —2a 


Substituting 2.52 for 2 a from statement (c), leads to equation (4). 


APPENDIX 2 
Loss Coefficients for Stepped Seal 


(K/K’')? in equation (12) was defined in equation (5) as the 
ratio between the pressure drop of a “regular” (A-type) and a B- 
type constriction, Fig. 5. Consider the B-constriction as a se- 
quence of two miter elbows with infinite aspect ratio and assume 
that velocities in the constriction itself will not differ materially 
from the velocities in the regular constriction. Denoting the 
total-pressure drop in this first elbow with Ap’ it can be written 


Ap u? 
(42) 
é 29 

where u? ~ us? in Appendix 1 and ¢ is the elbow loss factor. 

No total-pressure losses will be taken into account in the second 
elbow; those losses are incorporated in the dissipation of the 
kinetic energy that takes place after the flow leaves the constric- 
tion. 

The static pressure will decrease in the second elbow in the 
amount of (C/t)?(u?/2 g), Fig. 5. Therefore the over-all pressure 
drop across the B-type constriction can be written as 


2 Reference [10] changes the co-ordinate system in course of the 
derivation whereas it was not changed here; hence ¢ = 7.67 in [10]; 
{11] shows the original test data. ; 
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(44) 


Ap 
(45) 

Equations (45) and (5) enable the evaluation of equations (12) 
and (13) if & and (C/t)? are known. 

The square of the acceleration coefficient C/t has been com- 
puted from Fig. 10 in [13]. Even though the curve was obtained 
on a potential-flow basis and this particular case is not verified 
with experiments in [13], nevertheless, other tested and rather 
similar cases reported in the same reference show good agreement 
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between theory and experiment. The curve was used, therefore, 
as a basis for the calculation of (C/t)?, Fig. 5. 

It can be seen from Fig. 5 that C = 0.100 in. is close to the case 
of a miter elbow with equal outlet and inlet areas (X/C = 1.0) 
and infinite aspect ratio. The loss coefficient £ for such an elbow 
should be around 2.0 according to Fig. 14 of [14]. Recomputed 
from tests, £ came to 2.5. Reference [15] indicates that & for the 
discussed case might be higher than 2.0. Therefore & has been 
taken as 2.5 for X/C = 1.0. 

On the other hand, when inlet velocity head compared with 
the exit head can be neglected, total pressure losses based on 
outlet velocity head are zero. In case of X/C = 6.0 the inlet 
velocity head would be around 3 per cent of the outlet velocity 
head. Therefore when X/C increases from 1.0 to 6.0, & drops 
from 2.5 to 0. As the course of this function should indicate an 
asymptotic approach, a hyperbola like curve was drawn in agree- 
ment with what was mentioned above and the curve (1 + &) 
shown in Fig. 5 constructed accordingly. 
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Ap ur C\? u? 
29 (48) 
As 
Ap, 
29 
so 
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Introduction 


= us advanced power reactor concepts envision a 
“mobile-fuel’’ consisting of a molten alloy, slurry, or solution con- 
taining the fissionable material. In fast reactor designs such as 
LAMPRE-II the mobile fuel remains within the core whence heat 
is removed by a pumped coolant. The fuel side of the core may 
be a single, connected vessel with the coolant in separate tubes 
passing through it; or the fuel may be contained in separate 
sealed ligaments with the coolant on the “shell side’’ [Fig. l(a 
and })]. A variant is a ligament design with the ligaments con- 
nected at one or both ends into headers [Fig. I(c)]. Economic 
considerations for a fast power reactor limit the inactive fissiona- 
ble material inventory and require the maximum utilization of the 
active inventory. Hence very small cores and high heat fluxes 
result, and the use of a circulating fuel system is prohibitive. 

Thermal power reactor concepts such as the LMFR also in- 
volve a mobile fuel. In this case, however, the core is sufficiently 
large and the critical loading sufficiently low that a circulating 

Contributed by the Nuclear Engineering Division of THe AMERICAN 
Socrety or MecHANICAL ENGINEERS and presented at the Nuclear 
Engineering and Science Conference, New York, N. Y., April 4-7, 
1960. Manuscript received at ASME Headquarters, January 19, 
1959. Paper No. 58. 
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system with external heat exchanger seems most desirable [Fig 
1(d)]. 

In all of these systems natural convection influences to a sig- 
nificant extent the velocity and temperature profiles and the heat 
flux distributions. The externally cooled, LMFR-type design 
exhibits a combined foreced-natural convection within the core 
where the heat is generated and within the fuel system in general. 
This has been treated in Reference [1].1. The internally cooled 
arrangement, typical of the fast reactor designs, where the liquid 
fuel is within a closed vessel and circulation is motivated solely 
by the internal heat generation due to fission is the subject of this 
paper. [Figs. l(a, b, or c)]. 

There are two types of liquid-metal natural convection of in- 
terest for internally cooled cores. These are: 


1 the flow and heat transfer within a single, vertical, small- 
bore ligament [Fig. I(a)], in which heat is generated by fission 
within the fuel and removed to the coolant through the ligament 
wall; 

2 the over-all flow and heat transfer within a larger vessel 
consisting of vertical, parallel, single-ligament passages which are 
either connected through headers at the ends or continuously 
connected [Figs. 1(b and c)]. 


'‘ Numbers in brackets designate References at end of paper. 


Nomenclature 


a = radius of test section aga AT max ¢. to wall) , l tial between centerline at any 
c, = specific heat = VK el given axial position and wall at 
: . Subseri dlies 
g = acceleration of gravity ' bottom Su script &o appli 
h fil Sy Sear f Ra, = Rayleigh number based on length to differential between center- 
= film coefficient for heat transfe ; : 
. : and maximum temperature line at top and wall at bottom. 
k = thermal conductivity differential; Subscript R applies to differen- 
! = length of test section tial from wall to fluid at any 
Nu, = Nusselt’s number based on ra- a agl*( Tmax €. to walt) radius r and at any axial posi- 
dius; ha/k VK tion z. 
Pr = Prandtl number; v/k i velocity in axial direction 
: : ; R,r = dimensional and nondimensiona u = nondimensional velocity in axial 
q, = nondimensional volumetric heat co-ordinates in radial direction direction = a?U/kl 
g T = temperature 
source ———— temperatun v = volume 
¢ = nondimensional temperature dif- = dimensional and nondimensional 


Q, = volumetric heat source—energy 
per unit volume 


ferential = ————-; without 


aga‘A7T co-ordinates in axial direction 


a = coefficient of volumetric expan- 


Qs: 2 
Ra, 


wall heat flux—heat flow/unit 
surface area 

wall heat flux at mid-point of tube 

Rayleigh number based on radius 
and maximum temperature dif- 
ferential; 
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subscript nondimensional tem- 
perature differential, wall to 
fluid centerline, at any given 
axial position. Subscript o ap- 
plies to differential between wall 
and centerline at top of tube. 
Subscript € applies to differen- 


sion 
nondimensional core thickness. 
1 — nondimensional 
boundary-layer thickness. 
thermal diffusivity = k/pc, 
kinematic viscosity 
density 
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The flow pattern pertaining to the second configuration is a 
superposition of the pattern found for a closed-loop with hot and 
cold leg [Fig. 1(d)] and that found for a single ligament (item 1 
above). The closed loop can be simply approximated if knowledge 
of the effective friction factor and density change with tempera- 
ture is available [2]. 

This type of behavior would be encountered if the rate of heat 
generation did not properly match the coolant flow at all core 
radii. The heat transfer and fluid flow behavior for the case of the 
single, sealed ligament is the primary subject of this paper, in 
which it is attempted to provide sufficient information for engi- 
neering evaluations. [The case is of application in itself and is 
necessary to the analysis of item 2.] 

Typically, the proposed passage diameters are very small (on 
the order of 50 mils). Although the volumetric heat source is 
very large, it is shown that the temperature differences between 
passage centerline and wall are approximately those predicted on 
the basis of pure conduction; i.e., considering the fuel to be a solid 
rod. Nevertheless, knowledge of the actual temperatures and 
velocities is important since, as shown, the following factors exist: 

1 There is substantial motion of the fuel (velocities 100 feet/ 
hour) for sealed ligament and even greater motion for parallel 
ligaments. Hence mass transport of the container material 
between hot and cold regions is possible. 

2 The axial distribution of wall heat flux is substantially dif- 
ferent from that of the volumetric heat source. Since wall heat 
flux is the limiting factor in many designs, realistic knowledge of 
its magnitude is mandatory. 


COOLANT 


LIQUID FUEL 


REACTOR CORE —~ 
PRESSURE VESSEL 


Fig. 1(b) Reactor with fuel 
contained in separate closed 
tubes 


Fig. 1(a) Reactor with coolant 
flow passing through tubes in 
the core 
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A 
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Fig. 1(c) Reactor with fuel contained in separate 
tubes connected at top and bottom by headers 
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3 Because of (1) disposition of the fission gases is affected. 
The fluid velocities are sufficient to overcome the effects of 
static diffusivity, to cause a possible hold-up of gas in the fluid 
since the downward velocities are similar in magnitude to the 
rising velocity of small bubbles, and perhaps to affect the possible. 
bubble growth on the ligament walls. 

4 For somewhat larger passage diameters, there is a signifi- 
cant effect upon the temperature distributions. The parame- 
ter delineating the effectiveness of the natural convection in- 
volves the sixth power of the diameter [3]. Thus small increases 
in diameter rapidly become of importance. 


Limiting Cases 

No analytical or experimental data are presently available for 
the case of interest; i.e., a liquid metal contained within a vertical, 
sealed ligament of arbitrary length to diameter ratio within which 
heat is generated (steady-state being maintained by the removal 
of this heat through the ligament walls to a coolant). However, 
certain limiting cases have been explored both experimentally and 
analytically and serve to provide a guide for the understanding of 
the more general case. 

A Boundary-Layer Soluti Aq Fluids. Analytical predic- 
tions and approximate experimental verification of the tempera- 
ture and velocity profiles for an identical configuration, but for 
aqueous fluids (Prandtl number about one),? have been given in 
previous papers by one of the present authors [3, 4]. The analy- 
sis was based upon the integrated relations for conservation of 
mass, momentum, and energy for a disk of differential thickness 
between planes normal to the axis. Temperature and velocity 
profiles (Fig. 2) were assumed in such a way that boundary-layer 
thickness is a function of axial position. The boundary layer as 
here used includes some upward velocity components. The 
general effect of low Prandtl number, typical of liquid metals, is 
to reduce the temperature differentials and velocities. 

There is a rising core of fluid along the centerline which has 
approximately constant temperature and velocity at a given axial 
position. Along the ligament walls there is a descending boundary 
layer wherein the temperature and velocity both vary rapidly 
as a function of radius. The temperature differential between 
wall and centerline is a maximum at the top and falls to ap- 
proximately zero at the bottom. The boundary-layer thickness 
generally grows as the distance from the top increases. It reaches 
a maximum near the bottom before becoming, theoretically, zero 
at the bottom. Asa result of these temperature and boundary- 
layer thickness variations, the wall heat flux is a maximum at the 
top even though the rate of heat generation for the vessel is uni- 


2 It was shown in Reference [5] that the assumptions are not justi- 
fied for Prandtl number below about 0.5. 
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Fig. 1(d) Thermal reactor with external heat ex- 
changer 
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form and the wall temperature constant. Typical wall heat flux 
distributions [4] are shown in Fig. 3. 

It was found that the over-all temperature differences could best 
be presented in a plot of nondimensional volumetric heat source 
versus nondimensional temperature differential. A typical plot 
[4] is shown in Fig. 4. Wall temperature distribution (assumed 
linear in the axial direction and constant around the circumference 
at each axial position) and axial heat source distribution are the 
parameters. The results are theoretically applicable to any 
length to diameter ratio. 

The curves terminate, depending on wall temperature distribu- 
tion, at a value of g, between 10? and 10%. This termination 
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Fig. 2 Test section nomenclature schematic 
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Fig. 3. Normalized wall conduction versus nondimensional axial posi- 
tion, constant wall temperature, uniform heat source distribution 
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Fig. 4 Nondimensional heat source versus over-all temperature dif- 
ferential, experimental and calculated data 


signifies the lower limit of applicability of the boundary-layer 
solution used [4]. At this g, value the boundary layer, at some 
point along the axis, has grown to the extent of filling the entire 
vessel. The assumed profiles then no longer apply. Also, the 
boundary-layer thickness can no longer be a function of axial 
position. 

If a method of analysis completely similar to the original pro- 
cedure [3] is followed, except that the boundary-layer thickness is 
fixed equal to the vessel radius, a difference equation relating the 
nondimensional temperature differentials, heat source, and axial 
position can be derived. The steps beyond those of the analysis 
of Reference [3] are in the Appendix (Long-Tube Analysis). The 
entire analysis is too lengthy to repeat here. 

For those cases where the nondimensional heat source is suf- 
ficiently low so that the boundary layer fills the entire vessel at 
some point, it is necessary to use a combined solution to obtain 
results for the entire vessel. The original solution [3] can be used 
for that portion of the vessel near the top where the boundary 
layer is growing. The “long tube’’ solution is then used for the 
remainder of the vessel. Such a combination assures that con- 
servation of mass, momentum, and energy are observed for any 
plane normal to the axis through the vessel at any axial position. 
To obtain consistency at the point of joining of the solutions it is 
necessary that the original boundary-layer solution give a dimen- 
sionless temperature differential of g,/4 at the axial position 
where the boundary-layer thickness becomes equal to the vessel 
radius. The implementation of the boundary-layer solution [3] 
involves a numerical procedure which has been programmed for 
a digital computer. This analysis is also limited, as was the 
original [3], to Prandtl number near unity. It may be noted that 
the situation is somewhat similar to that of pipe flow where the 
boundary layer grows in the entry section (analogous to the liga- 
ment top) until it occupies the entire pipe. Henceforth it is 
“fully developed’’ and no longer grows with axial distance. 

A typical plot of nondimensional temperature difference be- 
tween centerline and wall as a function of axial position is shown 
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Fig. 5 Temperature versus axial position, constant wall temperature 


in Fig. 5 for two heat source strengths at constant wall tempera- 
ture. If a variable wall temperature is assumed, the wall to 
centerline temperature differential in the center section will be 
unchanged and the top and bottom will be changed only slightly 
from the constant wall temperature case, provided axial heat 
conduction remains negligible. It is noted that the dimension- 
less temperature differential is zero at the bottom and becomes 
asymptotic to q,/4, approaching this value very rapidly for small 
g». The q,/4 value of the nondimensional temperature differen- 
tial corresponds to the case of static conduction. 

B Infinite Length/Diameter Case—Aq or Liquid Metal. At the 
opposite extreme from the boundary-layer case [3, 4] is the case of 
infinite length to diameter ratio. An examination of the parame- 
ter grouping which forms q, discloses that this case corresponds to 
zero nondimensional heat source. This case has been discussed 
in detail by Murgatroyd [6]. A short examination of the physical 
situation will disclose the significant features. 

If the flow is laminar and the tube infinitely long, heat transfer 
normal to the axis is only by conduction. There are no convective 
effects either from turbulent mixing or radial transport. This 
case corresponds to ‘rod flow’’ (if axial temperature gradients 
are relatively small) for which the Nusselt number based on 
diameter and mixed mean temperature is 8 [7]. It is shown in 
the Appendix that the Nusselt number corresponds to a non- 
dimensional temperature differential between wall and center- 
line equal to q,/4 at a given axial position. The fully developed 
solution was asymptotic to this value for low q,. Hence we have 
an independent verification. 

The assumption of laminar flow is well justified for the fast 
reactor applications. The transition [4] from a generally laminar 
condition to a generally turbulent one appears to occur at a 
Rayleigh number based on a vessel radius of about 4 X 10’, cor- 
responding to q, of about 108 for the tests conducted. The 
corresponding Rayleigh number and q, for a typical fast reactor 
application are about 50 and 10, respectively. 

C Anticipated Behavior of Liquid Metals as Compared With Aqueous 
Fluids. No analysis or experiment applicable to the boundary- 
layer regime has been made for fluids with low Prandtl number 
in the applicable geometry. However, the infinite length analysis 
applies regardless of Prandtl number. It seems certain that the 
general type of flow behavior observed and predicted for aqueous 
fluids with high q, would also be observed with liquid metals. 
However, low Prandtl number should result in the temperature 
difference for a given heat source being reduced and the tem- 
perature gradients extending further into the fluid. As a result 
of the small temperature differences, small coefficient of volu- 
metric expansion, and high conductivity, the Rayleigh number 
would be much lower. Also, the forces motivating natural con- 
vection velocities would be reduced. Nevertheless, the asymptote 
reached as 9g, approaches zero is the same. On the basis of these 
qualitative considerations, it seems likely that the general type of 
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behavior would be the same, but that the fully developed flow 
condition would be reached at a somewhat higher q, since the 
temperature gradients would extend further into the fluid. 


Application of Results to Small-Bore Ligament With 
Liquid Metals 

Utilizing the limiting cases discussed in the previous section, it 
is possible to delineate the significant results to be expected from 
natural convection in a small-bore ligament filled with liquid 
metal fuel in a typical fast power reactor design. 

A Temperature Differentials—Centerline to Wall. Centerline to 
wall temperature differentials for aqueous fluids within the range 
of the boundary-layer solution (g, above 10%) have been given pre- 
viously [3, 4, 8, 9]. These are not directly applicable to liquid 
metals but the type of behavior for high qg, is believed typical. 
The present liquid metal fast reactor concepts, however, exhibit 
4, values within the fully developed range. Fig. 5 shows the non- 
dimensional temperature differential as a function of axial position 
and q, for constant wall temperature. The significant numerical 
results are listed in Table 1. It is noted that the dimensionless 
temperature differential is equal to approximately g,/4 (pure 
conduction value) for most of the vessel length, but shows an in- 
crease near the top and a decrease at the bottom. The increase 
is a function of g, and reaches a maximum of about 1.4 X q,/4 at 
9g, = 10%. This increase is very local (occurring within 6 per cent 
of the top) and may not be of real significance since a small axial 
heat flow in the tube wall could relieve any resultant over-heating. 
As far as the central portion of the tube is concerned, the results 
apply directly to any axial wall temperature distribution provided 
the axial gradient is small compared with the radial gradient. The 
results at the tube ends apply strictly only to constant wall tem- 
perature. However, they are typical of any distribution since the 
substitution of different temperature configurations produce dif- 
ferences only of degree. Also, the end results apply only to aque- 
ous fluids, while the central portion is of general applicability. 
However, fluids of low Prandtl number should exhibit similar be- 
havior with end effects somewhat reduced in magnitude and axial 
extent. 

B Wall Heat Flux. The axial wall heat flux distribution for 
uniform heat source, constant wall temperature, and aqueous 
fluids for the fully developed condition is shown in Fig. 6 and 
Table 1. Fig. 3 showed similar results for the boundary-layer 
solution at high q,. In both cases there is a decrease of wall heat 
flux at the bottom, a generally constant region in the central 
portion, and a sharp rise toward the top. These effects become 
very much more localized as q, is reduced into the fully developed 
region. The central portion of the curves are applicable for any 
fluid and wall temperature distribution, again with the qualifica- 
tion that the axial gradients are relatively small. The end por- 
tions apply strictly to aqueous fluids and constant wall tempera- 
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Table 1 Nondimensional temperature differential and normalized wall ture. However, a local rise of somewhat similar magnitude is to 
conduction for fully developed flow be expected with liquid metals and also with various wall tem- 

F Normalized wall* Boundary-layer perature distributions. Again, this increase of heat flux may not 

é conduction thickness/radius yp physically significant because of its very local character. Local 
qe = 10% y over heating so encountered may well be relieved by axial heat 
flow. The problem is of potential importance because the ther- 
mal stresses due to the wall heat flux are the limiting factor in 
many designs. 

C Velocity. Predicted and observed velocities for aqueous 
fluids in the boundarylayer regime have been given [3, 4, 8, 9]. 
Reference [6] shows analytical predictions for laminar and turbu- 
lent flow in vessels of infinite length. The maximum velocity to 
be expected in the laminar fully developed regime as a function of 
axial position and q, is shown in Fig. 7. As mentioned, the cases 
of interest seem limited to laminar flow. Except for the extreme 
ends, the values apply to any fluid and approximately to most 


0.996 

0.984 

0.920 

0.826 

0.456 
0 


* Temperature normalized to q,/4 and wall conduction normalized 


. to mid-point of the tube. 
** Fully developed flow from this point on. 


Normalized wall* Boundary-layer 
conduction thickness/radius 
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0.308 
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0.889 
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0.969 
0.975 
0.981 
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1.000 


* Temperature normalized to q,/4 and wall conduction normalized 0 = 
to mid-point of the tube. 
** Fully developed flow from this point on. 
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Fig. 6 Normalized wall conduction for constant wall temperature and 
uniform heat generation—fully developed flow 
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Fig. 7 Nondimensional boundary layer and core velocity versus axial position, 
constant wall temperature 
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wall temperature distributions. The end values should be correct 
in the direction of shift and its order of magnitude. Under the 
assumed velocity profiles, the maximum velocity in the fully de- 
veloped case is the ascending velocity along the centerline. How- 
ever, this may be merely a result of the assumptions. In the high 
q, boundary-layer cases the descending velocity adjacent to the 
wall was the maximum. The order of magnitude of velocity is 
believed to be essentially correct. This value is of prime im- 
portance from the viewpoints of erosion and mass transport. 


Conclusions 


No directly applicable theoretical or experimental results are 
available for natural convection in vertical liquid-metal filled 
ligaments encountered in mobile-fuel fast power reactor con- 
cepts. However, it is possible to delineate the general nature of 
behavior expected with respect to temperature, velocity, and 
wall heat flux profiles from limiting analyses and available ex- 
periments. These data are illustrated in the paper. It is shown 
that natural convection may be of importance from the view- 
points of perturbation of wall heat flux*distribution (a limiting 
design condition), and motivation of velocities important as a 
possible mechanism for mass transport, and hence is also an in- 
fluence on fission gas disposition. It does not appear that natural 
convection will substantially affect temperature differentials in 
the cases of interest. These can be adequately estimated as- 
suming pure conduction. 
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APPENDIX 
Long Tube (Fully Developed Flow) Analysis 


In order that the solutions will be consistent at the transition 
points, the same arbitrary velocity and temperature profiles are 
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assumed as in Reference [8]. However, it is assumed that the 
boundary-layer thickness is no longer a function of z but is con- 
stant and equal to the tube radius. 

From Reference [8], for constant wall temperature, the as- 
sumed temperature and velocity profiles are: 


for l1>r>8B 


tg =t for B>r>0 (1) 


+ Or — | for 1S 6 


(2) 
u=-y for B>r>0 
but in this case 8 = 0, therefore 
tr = r?) (3) 
—y[l — + — (4) 


= 


Substituting equation (4) into the continuity relation (Ref- 
erence [8]) gives 


1 1 Ps 
ruar = 0 = rar rdr 
rar +8 ff rart (5) 


from which 6 = —5. Equation (4) can then be written as 
u = —y(5r? — Gr? + 1) (6) 


0, = —y¥(15r? — 12r) = 0 


= 
dr 

Now, substituting the velocity and temperature profiles, Equa- 
tions (3) and (4) into the momentum equation (Reference [3] ), 


ou 
= ( 
+ ( ) (8) 


and the velocity, Equation (6) becomes 


(10) 


= 


t 
“42 (5r3 — 6r? + 1) 


Repeating the same procedure with the integrated energy 
equation (Reference [8] ) 


of ot 
- (=) 
or 0 Or 2 


1 
rutdr = 
or Jo 
or rel 


(11) 
we have 


(t?) 


12 
336 Or 


(13) 
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| 
and at 
1 t 
2 
we get 
t 
(9) 
- 
= 


and Equation (11) reduces to 


(14) 
336 Or 2 
— 
axial convection wall heat 
conduction source 


for constant wall temperature. 

From Equation (14) it can be seen that ¢ approaches q,/4 as a 
limit as the tube length becomes infinite (0/drz +0). The axial 
convection term is a function of the axial temperature and the 
wall conduction term is a function of the wall to centerline dif- 
ferential temperature. For variable wall temperature, it can be 
shown by a similar analysis that 


1 
336 ar & +e 2 ° 


(15) 


where 


(16) 


*@=t+t, 


Therefore 


(¢ + — (t + 
Ar 


+ 6723 + 1689, = 0 (17) 


and 


i = ty + ty 
ty 


Solving by quadratic formula: 
= —(tww-1 + 168Ar) + + 168Az)? 
+ (ty + toy)? — 336Azty + 1689,Azr — toy, 


(18) 


where Az and ¢ are positive numbers 


Natural Convection Laminar Flow in Infinite Length Tube 
The heat flow in the axial direction must be negligible com- 
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pared with flow in the radial direction for the case of pure con- 
duction, and the heat lost through the wall in a given axial section 
must equal the heat generated in that section 
Then 
Q,r? 


——— + C, Inr + C, (Reference [2], page 662) (19) 


T = 
4k 


= (R- 1) (20) 
4k 

For such a parabolic distribution in a cylirider it is well known 

that the mean temperature differential between fluid and wall is 


Twa: 
1/, the maximum differential, i.e., McAdams, 


Reference [7], page 233, shows that for an infinitely long tube 
with laminar rod-like flow and uniform heat flux the Nusselt 
number based on the diameter and the mean temperature 
approaches 8 away from the entrance section. 


ha 


mean k 


(21) 


hD 
NuD. = 7% = Sor Nua.r 


Thus based on the radius and the temperature difference from 
wall to centerline 
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Ia the present analysis using the nondimensional heat source 
term 
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as may be shown by simple algebraic substitution for g, and ¢ in 
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Selection of the Hydraulic Turbines 
for the Mammoth Pool Project 


Mammoth Pool Project has a normal static head of 1100 ft which is within the operating 
range of either the impulse or reaction type hydraulic turbine. An economic analysis 
is made for using two 88,000 hp turbines, either of the multinozzle impulse or of the reac- 
tion type. Power generation and motoring operations of each type are evaluated with the 
final evaluation justifying the selection of two reaction turbines for this powerhouse. 


Introduction 


Ta Southern California Edison Company owns and 
operates 25 hydroelectric plants which have an effective operating 
capacity under optimum conditions of approximately 637,000 kw. 
This is about 18 per cent of the total effective operating capacity 
of the combined hydroelectric and steam electric plants owned 
and operated by the Southern California Edison Company. 
Seven of these hydroelectric plants having an effective capacity 
of 540,000 kw comprise the Big Creek Hydroelectric Develop- 
ment located in Central California, 248 miles north of Los 
Angeles. 

The static heads on these hydroelectric plants vary from a 
maximum of 2418 ft to a minimum of 230 ft. The highest head 
plants have horizontal shaft impulse turbines of the single-jet, 
double-overhung runner type. The highest head reaction type 
hydraulic turbines have a 827-ft static head. 

The Mammoth Pool Project is the newest addition to the Big 
Creek Hydroelectric Development being scheduled to be com- 
pleted early in 1960. This project is in the west portion of the 
Big Creek Development, being located on the upper San Joaquin 
River in Madera and Fresno Counties, approximately 80 miles 
northeast of the City of Fresno, Calif. The powerhouse is on the 
San Joaquin River, at the upstream end of the existing afterbay 
of the Edison Company’s Big Creek Powerhouse No. 8. Water is 
supplied to the powerhouse from a 123,000-af reservoir through an 
8-mile long tunnel and a 2072-ft long penstock. The electric 
energy generated will be integrated with that produced by the 
existing Big Creek Hydroelectric Plants. 

The hydraulic design criteria for the Mammoth Pool turbines 
listed in Table 1 show that the maximum normal static head on 
the turbines is about 1100 ft and the normal effective operating 
head, based on one-third draw-down of the reservoir, is ap- 
proximately 950 ft. From the head standpoint, the turbines 


Table 1 Hydraulic design criteria 


Water surface elevation 
Mammoth Pool Reservoir - ft 3330 
Static head - ft 
Operating head - ft 94 


Discharge for 2 turbines,cfe 1980 


Expected maximum plant 


output - hp 204,000 1B4000 146,000 


Contributed by the Prime Movers Subcommittee of the Hydraulic 
Division and presented at the Gas Turbine Power & Hydraulic 
Conference, Houston, Texas, March 6-9, 1960, of THe AMERICAN 
Society or MECHANICAL ENGINEERS. 

Note: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those of 
the Society. Manuscript received at ASME Headquarters, December 
3, 1959. Paper No. 60—Hyd-10. 
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of this project could be either of the impulse or reaction type. 
This paper will not include the selection of turbine size, speed, 
nor other details such as the turbine setting, but the economic and 
operating characteristics of each of the two types will be re- 
viewed and an evaluated annual cost of each will be estimated for 
use in the Mammoth Pool Project. 

The following is a brief description of the two types of turbines 
that were considered: 


Impulse Turbine Arrangement 

The impulse turbine studied is of the vertical shaft, 6-nozzle 
type operating at 200 rpm and developing 88,000 hp at 940-ft 
head. Fig. 1 shows the powerhouse which would be required for 
the impulse turbines. The turbine spiral case is set with the 
horizontal center line of the case and the jets at elevation 2238 
which is 8 ft above the normal tailwater elevation of 2230. The 
runner is an alloy steel casting and of the integral type. 
The runner has an approximate pitch diameter of 130 in. and 
weighs about 30 tons. It is believed that a runner of this size 
made up of bolted buckets, if not entirely impracticable, would 
not give as satisfactory service as the integrally cast type. Each 
nozzle is equipped with a cast steel jet deflector capable of de- 
flecting the jet in less than two seconds. As it may be necessary 
to operate the turbines during flood conditions when the tail- 
water is 13 ft higher than normal, that is, at elevation 22438, 
which is 5 ft above the center line of the runner, the turbines 
are equipped with a pneumatic tailwater suppression system. 


Reaction Turbine Arrangement 

The reaction turbine studied is of the vertical shaft type operat- 
ing at 360 rpm and developing 88,000 hp at 950-ft head. Fig. 2 
shows the powerhouse which would be required for the reaction 
turbines. The turbine spiral case is set with the horizontal center 
line of the case at elevation 2225 which is 5 ft below the normal 
tailwater elevation of 2230. The pressure relief valve is of the 
cylindrical sliding sleeve, fixed cone type. This valve is at the 
same elevation as the spiral case and capable of discharging 100 
per cent of the turbine flow. For reasons which will be con- 
sidered later, the turbine is equipped with a pneumatic suppres- 
sion system for depressing the water level in the draft tube regard- 
less of the elevation of the outside tailwater. 


Operating Conditions 
To satisfy the requirements of this hydroelectric plant, the 
hydraulic turbine-generator units must be capable of operating 
efficiently as either a base load, a peaking, or a regulating plant. 
Three methods of operating the turbine-generator units for 
generating electrical energy are considered, namely: 


1 “Regulated stream flow” 
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Fig. 1 Mammoth Pool Powerhouse—plan and section required for impulse turbines 
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Fig. 2 Mammoth Pool Powerhouse—plan and section required for reaction turbines 


2 Peaking, using two units 
‘ 2000 - 
3 Peaking, using either one or two units. a | | | | 


| 
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In addition, two conditions of “motoring” the turbine-generator ! 


units are considered, namely: REGULATED MONTHLY 


MEAN DISCHARGE 


Fig. 3 shows the flow-duration curve for the release of stored 
water for the Mammoth Pool Project based on a 30-year water 
study and an assumed optimum storage regulation. The criteria 
followed for “regulated stream flow’? was based on using the 
water quantities shown in Fig. 3 and operating the plant cor- 
respondingly at varying output conditions. For total water 
flow rates of about 1000 cfs and higher, both turbine-generator 
units were considered to be in operation. For plant flow rates 
less than approximately 1000 cfs, only one turbine-generator unit 
was considered to be operating. 

Operating the plant for peaking operation using two turbine- wma T 
generator units considered the same flow duration curve given by 
Fig. 3 with respect to the total quantity of water released by the 
plant. The operating time of the units, however, was shorter; 
the minimum total flow rate for the two units being 1800 cfs. The ; 
operating time chosen for operating the turbine-generator units, } 
was an amount that gave the same total quantity of water re- | 
leased as represented by Fig. 3. { | | | | | 

Operating the plant for peaking operation using one or two o5 : 
turbine-generator units was similar to the previous case except 
that, when the block of water released by the regulated stream PER CENT OF TIME EQUALED OR EXCEEDED 
flow method equaled that required by one turbine-generator unit, Fig 3 Mammoth Pool Project—mean monthly flow-duration curve 
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then this peaking method considered only one unit in operation. 
The minimum water flow rate for one turbine was taken as 900 
cfs. 

Fig. 4 shows the turbine performance characteristics for each 
type turbine for two head conditions. As noted in Fig. 4, each 
turbine has the same maximum efficiency of 90 per cent. How- 
ever, the point of maximum efficiency is different for each type. 
For the normal design conditions it is apparent that, for turbine 
operation of about 74,000 hp or more, the reaction type tur- 
bine will develop more electrical energy than the impulse turbine 
for the same quantity of water. Similarly, for turbine outputs 
below 74,000 hp, the impulse turbine will be capable of generating 
more electrical energy than the reaction turbine for the same quan- 
tity of water. The multinozzle impulse turbine is capable of main- 
taining high efficiencies at low powe* output by decreasing the 
number of jets in use. This is shown in Fig. 4 where the change 
from six jets to four jet operation is made at about 50,000 hp 
when operating at 940-ft head. 
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Fig. 4 Turbine perf impulse and reaction turbine-efficiency 
versus power output 


It is readily seen that the total electric energy generated from 
the same rate of water flow by either the impulse or reaction tur- 
bine is going to be dependent on the method of turbine operation. 
That is, when considering the turbine performance curves of Fig. 
4 for heads of 940 to 950 ft at power outputs above 74,000 hp the 
reaction turbine has a greater output. Similarly, at power out- 
puts below 74,000 hp the impulse turbine has a greater output. 

Detailed computations were made for each of the afore-men- 
tioned three methods of plant operation to determine the amount 
of electrical energy that could be generated when the generators 
were driven by either the impulse or reaction turbines. The re- 
sults of these computations are given in Table 2. In comparison 
with the reaction turbines, the impulse turbines are capable of 
generating the greatest annual amount of electrical energy totaling 
575,000,000 kw-hr, when the plant is operated as a regulated 
stream flow plant. Conversely, either of the two methods of 
operating the plant on similar peaking conditions indicate that 
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the reaction turbine will generate more electrical energy than the 
impulse turbine. Operation of the plant for peaking and using 
either one or two turbines, indicates that 573,000,000 kw-hr will 
be generated by the reaction turbines. This is the best method 
of operation for the reaction type turbines. 

In summary, from the standpoint of capability, Table 2 shows 
that, if the plant is operated on the regulated stream flow basis, 
the impulse type turbines will produce 9 million kw-hr more than 
the reaction turbines. However, if the plant is. operated as a 
peaking plant, the reaction turbines will produce from 3 to 5 
million kw-hr more than the impulse turbines. 


Motoring Operation 


Current Company operating policy requires that some of the 
hydraulie turbine-generator units in the Big Creek area at times 
“motor”’ either as spinning reserve or as synchronous condensers 
dependent on system requirements. The method of operating the 
turbines for each motor condition is similar although there are 
slight variations. When a turbine-generator unit is motoring as 
a spinning reserve, operational requirements demand that the 
unit be capable of generating electrical energy within a time 
interval not exceeding several seconds. If, however, the turbine- 
generator unit is motoring as a synchronous condenser, the only 
requirement is that the unit add reactive kva to the electrical 
transmission system, practical instantaneous load pick-up not 
being required. 

Past experience indicates that the units can be expected to 
motor about 25 per cent of the time. The motoring time will be 
about equally divided between the two motoring conditions. 
However, on regulated stream flow operation, the generating time 
of the units is of such a duration that the remaining time available 
for motoring will be only 19 per cent. 

Impulse Turbine. The impulse type turbine has several ad- 
vantages over the reaction type turbine for “motoring”? opera- 
tions at this particular hydroelectric plant. The runner of the 
impulse turbine will always be revolving in air while the generator 
is motoring and the windage losses of the runner will be low. It 
is normal to have full penstock pressure at the impulse nozzles, 
and, with the needles in the closed position, the water leakage is 
negligible. In addition, an impulse turbine does not cal! for any 
water comparable to that required by seal rings of a reaction tur- 
bine. With full penstock pressure against the closed needles the 
impulse turbine unit is capable of generating electrical energy in 
a matter of seconds while being motored either as a spinning re- 
serve or as a synchronous condenser. 

Reaction Turbine. Reaction type turbines for motoring opera- 
tions are successfully used in existing plants of the Big Creek 
Hydroelectric Development. However, it is not as advantageous 
to use the reaction turbine for these operating conditions as its 
losses are greater than those of the impulse turbine. 

The ability of a reaction-turbine-generator unit, while motoring 
as a spinning reserve, to produce full output in a minimum of 
time necessitates full penstock pressure against the turbine 
wicket gates while the gates are in the closed position. A. definite 
water loss exists through the wicket gates as it is not practical 
to have the wicket gates drop-tight under the high operating 
head of the Mammoth Pool Powerhouse, the water loss varying 
as the square root of the head maintained in the scroll case. 
To minimize such losses through the wicket gates, the wicket 
gates are cast of 13 per cent chrome steel and machined to close 


Table 2 Annual electrical energy generated in million kw-ht 
Impulse Turbines Reaction Turbines 


"Regulated stream flow" 575 S66 


Method of Operating Rydroelectric Plant 


Peaking - always using two units 559 


Peaking - using either one or two unite 568 573 
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clearances. It is possible to operate the turbine with the turbine- 
shut-off valve (at the entrance to the spiral case) in a closed 
position. However, operating the turbine in this manner would 
not meet the system requirements as the delay in opening the 
turbine shut-off valve and the control problems would be un- 
satisfactory for this plant. 

Another source of water loss on the reaction type turbine is 
due to the upper and lower seal rings. The only way the water 
loss in the seal rings can be maintained at a low value is by keep- 
ing the seal ring clearances at a practical minimum, in this case 
in the order of 0.015 to 0.020 in. 

Motoring the reaction turbine for synchronous condenser opera- 
tion is identical in all respects to the foregoing spinning reserve 
method of operation, except water pressure is maintained in the 
turbine spiral case at a lower head. In this method of operation 
it is possible to close the turbine shut-off valve. Pumped water 
from the draft tube maintains the spiral case pressure at 5 to 10 
psig. With such low pressures the water losses through the 
closed wicket gates will be materially less than the corresponding 
losses when the turbine is in the spinning reserve operation. 


Motoring Losses 


As the Southern California Edison Company has had no pre- 
vious operating experience on either an impulse or reaction tur- 
bine of the size being considered for the Mammoth Pool Power- 
house, it was necessary to estimate and extrapolate existing field 
data to determine turbine-generator losses for the motoring con- 
ditions corresponding to both spinning reserve and synchronous 
condenser operation. 

Impulse Turbine. ‘The measured power losses for motoring a 
40,000 kw single-jet, double-overhung runner, hydraulic turbine- 
generator unit at a Big Creek Hydroelectric Plant having a static 
head of 2418 ft amount to 963 kw. The motoring losses for the 
Mammoth Pool impulse turbine were, therefore, assumed to be 
about the same, or 970 kw. Other studies indicate that the 
evaluated incremental energy cost at the Mammoth Pool Plant 
for comparative purposes is about 4 mills per kw-hr (by incre- 
mental energy cost in this case is meant the equivalent steam re- 
placement based on the cost of fuel only). The cost of electric 
energy for each hour one impulse turbine is motored equals, there- 
fore, 970 K 0.004 = $3.88. If the two turbine-generator units 
are motored ten per cent of the time, the annual cost of electric 
energy for motoring is valued at 3.88 X 2 K 876 = $6797.76, 
or, say, $6800. There are no additional energy requirements or 
water losses to be considered for the impulse turbine evaluation. 
There may be periods where the tailwater will be high enough 
to drown out the turbine runner unless a pneumatic tailwater sup- 
pression system be used but the electrical energy used by such a 
system is negligible and, therefore, is not given further considera- 
tion in this comparison. 

Reoction Turbine. In the case of reaction turbine-generator unit 
the losses when motoring the unit are higher than the comparable 
losses when motoring an impulse turbine-generator unit. The 
reaction turbine losses increase as the distance from the bottom of 
the turbine runner to the water surface in the draft tube de- 
creases and become quite large when the runner is in contact with 
the water. Tests on dther reaction turbines in the Big Creek 
project indicate that the water level in the draft tube should be 
several feet below the bottom of the turbine runner in order to 
minimize the losses. To determine these losses tests were per- 
formed on a reaction turbine-generator unit having an effective 
capacity of 42,000 kw at a speed of 257 rpm and a static head on 
the turbine of 416 ft. Lowering the water in the draft tube of 
this turbine from an elevation of 1.2 ft below the bottom of the 
runner to 5.2 ft decreased the power consumed by the turbine- 
generator unit from about 3000 kw to 1500 kw. Venting the draft 
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tube to the atmosphere at the lower water level decreased the 
motoring losses only by another 60 kw. When the water in the 
draft tube was closer than 1.2 ft from the runner bottom the 
motor losses increased proportionately. 

Motoring losses on two other reaction type turbine-generator 
units at Big Creek have been measured at 1790 kw and 1738 kw. 
Each of these units is rated at 25,000 kw and a static head of 827 
ft. 

Based on the foregoing three units, the motoring losses for the 
reaction turbine-generator units considered for the Mammoth 
Pool Project are taken as 2000 kw. On the basis of an evaluated 
incremental energy cost of 4 mills per kw-hr the cost of motoring 
the turbine-generator unit per hour of motoring equals then 2000 
x 0.004 = $8.00. 1° the two turbine-generator units are motored 
ten per cent of the time, the annual cost of electric energy for 
motoring is valued at 8.00 X 2 X 876 = $14,016, or, say, $14,000. 

As the turbine runner is below the normal tailwater elevation, 
it becomes necessary to depress the water in the draft tube to an 
elevation of at least five feet below the runner or the motoring 
losses of the turbine-generator unit will substantially exceed the 
estimated value of 2000 kw. The 20-hp air compressor and air 
storage tank system used to depress the water in the draft tube 
for motoring operation has an annual evaluated power consump- 
tion of less than $100. This small cost is negligible and not con- 
sidered further in this evaluation. 

Instead of using a pneumatic method of tailwater suppression 
for motoring operations, it would be possible to lower the draft 
tube gates and maintain the draft tube in an unwatered condition 
by use of the sump pumps. This type of operation would not 
comply with the rapid load pickup requirement of spinning re- 
serve operation. However, this method could be followed for the 
synchronous condenser operation except that the required con- 
trols make the method somewhat complicated and are therefore 
not considered. 

The value of the water lost when motoring the reaction turbine- 
generator unit is an appreciable amount. As previously men- 
tioned, it is necessary to have water lubrication for the seal rings 
and, in addition, there is a loss of water through the turbine 
wicket gates. Based on operating experience it is estimated that 
these water losses may total as much as 28 cfs for each reaction 
turbine with full penstock pressure against the closed turbine 
wicket gates. Of this total less than one cfs is required by the 
seal rings while the major part, or 27 cfs, is being lost through 
the closed gates. It is believed that the 28 cfs water loss for a 
reaction turbine is high as the writer has information from other 
sources indicating that these losses are in the order of one-quarter 
of this amount. 

On the basis of average head and turbine efficiency conditions, 
the evaluated annual energy corresponding to the 28 cfs water 
loss amounts to approximately $710 per turbine for motoring one 
per cent of the time, or $14,200 for both turbines motoring ten 
per cent of the time. 

With the turbine shut-off valve closed and the spiral case pres- 
sure maintained below 10 psig by means of pumped water, it is 
estimated that the total water loss will be about 5 cfs. This water 
is evaluated to cost (5/28) * 710 = $126.78, or about $130 per 
turbine for motoring (as synchronous condenser) one per cent of 
the time, or $2340 for both turbines motoring nine per cent of the 
time. The annual energy used by the 15-hp pump which main- 
tains water pressure in the spiral case during synchronous con- 
denser operation is evaluated at less than $50 per turbine. 
This small amount is neglected in the final evaluation of the two 
turbines. 


Maintenance Costs 


An evaluation of the two types of turbines would not be com- 
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Fig. 5 Mammoth Pool Powerhouse—plan and section as built showing reaction turbines 


plete without consideration being given to the maintenance 
costs. Based on judgment and experience on operation of similar 
hydraulic turbines the following maintenance cost values are 
assigned to these turbines. 

Impulse Turbine. Assuming that the vertical shaft, multinozzle 
impulse turbine would not differ greatly in maintenance costs 
from the single jet, double-over-hung-runner, impulse turbine, 
and giving due regard to the increased number of jets, the annual 
maintenance costs per turbine are estimated to be $1950. This 
amount is considered necessary to maintain the turbine in an ef- 
ficient operating condition and includes runner repairs but no al- 
lowance for runner replacement. 

The impulse turbine runner will be an integral casting of carbon 
steel. The Company service experience with integral cast wheels 
is limited when compared with its experience with bolted bucket 
type runners. However, it is estimated that the life of the inte- 
gral runner will be about 15 years. The approximate cost of one 
runner is $94,000, corresponding to an annual charge of $6300 
per turbine. The total annual maintenance cost for both impulse 
turbines is then $16,500. 

Reaction Turbine. In comparison with the impulse turbine it is 
more difficult and costlier to maintain the original efficiency of the 
reaction turbine. The maintenance costs of the reaction turbine 
increase with the higher heads and the necessity of having a 
minimum water loss through the closed turbine wicket gates. It 
is estimated that the annual maintenance cost per reaction tur- 
bine will be $3600. This amount includes turbine runner repair 
but no allowance for runner replacement. The annual turbine 
overhauls include careful inspection of the turbine runner and 
corrective measures are taken, such as welding, making slight 
modifications to runner contours, and overlaying cavitation areas 
with stainless steel. Experience indicates that when this type 
of attention is given to a turbine runner its life may exceed 35 
years. Turbine runners on existing reaction turbines in the Big 
Creek Development have never been replaced except as dictated 
by changing the operating frequency of the turbine-generator 
units from 50 to 60-cycle operation. Consequently, the annual 
replacement cost of the reaction turbine runner is not given 
further consideration in this analysis. 


Powerhouse Structural Costs 


The final item to be evaluated is the cost of the powerhouse 
structure and major equipment. Fig. 1 shows the general over-all 
size of the powerhouse if impulse turbines are used. Fig. 2 shows 
similar features for the powerhouse if reaction turbines are used. 
Comparing Figs. 1 and 2, it is evident that the impulse power- 
house is the larger of the two. The increased size of the impulse 
powerhouse is due to the larger spiral case, 43 ft diameter, as 
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compared to the reaction turbine spiral case diameter of 25 ft. 
In addition, adequate space or clearance must be provided for re- 
moving the needle assembly from the impulse turbine nozzles. 
However, the reaction turbine due to its lower setting requires a 
deeper excavation. 

Table 3 shows actual turbine and generator costs for each type 
of powerhouse. The estimated powerhouse structural costs and 


* the auxiliary equipment costs are also shown. For the items in 


Table 3 there is a net saving of $595,555 in initial construction 
cost if the reaction turbines are used. Table 3 also shows the 
annual fixed charges for these items. For this evaluation the an- 
nual charges are considered to be 11.5 per cent of the fixed 
costs. This percentage represents the charges for interest, de- 
preciation, taxes, and insurance for this project. 

Fig. 5 shows the powerhouse as built and includes still addi- 
tional savings not considered in this analysis which makes the 
reaction turbine installation even more attractive from a cost 
standpoint. 


Annual Cost Summary 


Table 4 summarizes the annual cost for regulated stream flow 
operation of the Mammoth Pool Project with the impulse and 


Table 3 Comparative turbine, generator, and powerhouse costs 
$2,178, 700 


Items Reaction Turbines 


Turbines and spare parts $1,670,810 
2, 311,000 
1,005, 500 

$5,495,200 


2,337,160 


891,675 
$4,899,645 


Generators 


Powerhouse and auxiliary equipment 
Total Fixed Cost 


Saving ip Construction Cost favoring 


Reaction Turbine Installation $ 595,555 
Annual Fixed Charges (0.115 x Total Fixed Cost) $ 563,450 $ 631,09 


Anmial Saving in Construction Costs favoring 


Reection Turbine Installation $ 68,500 


Table 4 Summary of annual costs for regulated stream flow plan 
operation 


Items Reaction Turbines 


$563,450 


Annual Fixed Charges $631,590 
Credit for Higher 
Output at & mille/kw-hr 


Motoring of time, Comei.ting of 
Spinning Reserve - 10% of time 
Energy Lovses 
Water Losses 
Synchronous Condenser - v6 of time 
Energy Losses 
Water Loeses 


nance 
Totel Evaluated Annual Cost $625, 27 
Saving in Total Evaluated Annual 


Stream-Flow Plant Operation Cust ‘favoring 
Reaction Turbine In.tallation 
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Table 5 Summary of annual costs for peaking plant operation using 
turbine-generator units 


Items Reaction Turbines se Turbines 
Annual Fixed Charges $553,450 $631,990 
Credit for Increased Energy Qutput 
et & aills/kw-br -12,000 
Motoring - 25% of time, consisting of 
Spinning Reserve - 13% of time 
Energy Losses 18,200 8,840 
Water Losses 18,500 ° 
Synchronous Condenser - 12% of i. 
Energy Losses 16,800 8,160 
Water Losses 3,120 ° 
Maintenance 7,200 16, 500 
Total Bvaluated Annual Cost $615,270 $665,450 


Saving in Total Bvaluated Annual 
Peaking Plant Operation Cost, 
Pavoring the Reaction Turbine Installation 
$ 50,180 


reaction turbines. From Table 2 the energy generated by the 
impulse turbine-generators is 9 million kw-hr more than that 
generated by the reaction turbine-generators, therefore, Table 4 
shows the value of this energy differential, or $36,000, as a credit 
to the impulse turbines. The other values in Table 4 have been 
discussed in preceding paragraphs. The total evaluated annual 
cost of the reaction turbine-generator arrangement has the lower 
cost by an amount of $11,580. 

Table 5 summarizes the annual cost for the peaking method of 
operation, based on the use of two turbine-generator units. As 
the energy generated by the reaction turbine-generator units is 
higher by 3 million kw-hr, the reaction turbines are credited with 
$12,000, which is the evaluated cost of this electrical energy. 
The remaining cost items have been covered previously. Taking 
all the cost items into account, the total evaluated annual cost 
of the reaction turbine-generator installation is less by an amount 
of $50,180. The other peaking condition considered will also 
show that the reaction turbine annual cost advantage is even 
greater than the peaking operation considered by Table 5. 


Conclusions 

This analysis of the two types of turbines for the Mammoth 
Pool Project shows definitely that the annual cost of ownership 
and operation of the reaction turbines is less than the comparable 
cost of the impulse turbines. Another definite advantage in us- 
ing the reaction turbines is their lower initial project cost. It 
is very probable on the other hand that, under different operating 
conditions requiring energy generation at low power outputs, 
multijet impulse turbines could have a definite annual cost ad- 
vantage. 

On the basis of lower evaluated over-all annual costs and lower 
first costs the reaction type hydraulic turbines were selected for 
the Mammoth Pool Project. 


DISCUSSION 
G. D. Johnson! 


Since the writer is employed by the company that manu- 
factured the two Francis-type turbines installed in the Mammoth 
Pool powerhouge, it is obvious that we agree with the author’s 
selection in this ease. In our proposal letter to the author’s com- 
pany, we oversimplified the problem as follows: “.... We 
believe that Francis-type turbines are more desirable for this 

' Chief Hydraulic Engineer, Hydraulic Division, Allis-Chalmers 
Manufacturing Company, York, Pa. Mem. ASME. 
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installation. They are more economical due to lower costs for 
turbines and generators and savings in powerhouse space re- 
quirements. Turbine efficiencies are higher on a net head basis, 
while on a gross head basis there is the further advantage of 
utilizing the entire head differential to the tailwater surface, 
whereas Impulse units would sacrifice 8 feet of head when the 
tailwater is at its normal elevation. During times of high tail- 
water, normal operation of the Francis turbines is not compli- 
cated by the necessity for depressing the water in the discharge 
pit.”” Although the author’s analysis would differ from ours in 
a number of details, his conclusions are the same. 

The static head of approximately 1100 feet would be suitable 
for either Francis-type or multijet vertical-shaft impulse tur- 
bines, depending upon local water and operating conditions. 
A 6-jet impulse turbine was selected about ten years ago? for 
the Bridge River Plant of the British Columbia Electrie Com- 
pany Limited with a static head of approximately 1200 feet, pri- 
marily because of the silt content of the water which has a par- 
ticularly destructive effect on Francis turbine runner seals. It 
was felt that maintenance of the impulse runners, nozzle seat 
rings, and needle tips would be much more convenient and less 
expensive. Some twenty years ago,’ a Francis turbine was 
selected for the Nantahala Plant of the Nantahala Power and 
Light Company in North Carolina with a static head of ap- 
proximately 1000 feet. Maintenance has been reasonable as 
stainless steel is quite durable in clean water, even at the high 
velocities due to the relatively high head. 

At the time when the order for these turbines was awarded, ours 
was a relatively small company specializing in hydraulic tur- 
bines, pumps, valves and related equipment and, because of our 
convictions as quoted above, we had not built any 4 or 6-jet 
vertical impulse units. For the past year, we have been operat- 
ing as the Hydraulic Division of one of our former large com- 
petitors in this field. Now, our combined backlog of experience 
includes a number of successful installations of all types of im- 
pulse turbines, as well as Francis, for the “controversial” head 
range above and below 1000 feet. As in the past, we will con- 
tinue to recommend the type of equipment that appears to be 
most suitable and economical for each particular installation. 

The writer wishes to commend the author for his presentation 
of this very interesting paper. We need more good papers from 
engineers who purchase and operate hydraulic machines to sup- 
plement those written by engineers who design and manufacture 
the equipment. This résumé of the various economic considera- 
tions and operating conditions that must be studied by a power 
company in order to decide upon the type of equipment for a 
particular project. should assist those responsible for future in- 
stallations (1.e., both users and manufacturers) to specify and/or 
to offer the most suitable machines for each installation. 


Author’s Closure 

After the preparation of the original analysis it has been pos- 
sible to check the assumption that the motoring losses of the 
turbine-generator units would be 2000 kw. Tests indicate that 
the electrical motoring losses for one unit on spinning reserve 
operation are 2300 kw and about 865 kw when the unit is operated 
as a synchronous condenser. 


2W. F. Boyle and 1. M. White, ‘62,000-Hp Vertical Six-Nozzle 
Impulse Turbines for the Bridge River Hydrodevelopment,’’ Trans. 
ASME, vol. 73, 1951, pp. 289-296. 

3 J. P. Growdon, R. V. Terry, and H. H. Gnuse, Jr., ‘‘Nantahala 
Turbine,’’ Trans. ASME, vol. 68, 1946, pp. 687-700. 


Transactions of the ASME 


; 
A 
ct 
= 
“y 


On the basis of actual test data there are additional annual 
savings applicable to Tables 4 and 5. The energy losses of the 
reaction turbines for spinning reserve operation for 10 per cent 
of the time become 2300 0.004 K 2 X 0.1 & 8760 = $16,118, 
say, $16,200. For operating the units as synchronous condensers 
for 9 per cent of the time the energy losses become 865 0.004 
xX 2 X 0.09 X 8760 = $5456, say, $5500. The annual energy 
losses for motoring the units 19 per cent of the time amount to 
$16,200 + $5500 = $21,700, which is $4900 less than the origi- 
nal value shown in Table 4. On this basis the saving in the 
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total evaluated annual stream-flow plant operation cost, favoring 
the reaction turbine installation, becomes $16,400, instead of the 
originally given value of $11,500. 

In a similar manner new values for Table 5 energy losses are 
computed as $28,400 for operating the reaction turbines 13 per 
cent of the time as spinning reserve and 12 per cent of the time 
as synchronous condensers. The saving in total evaluated annual 
peaking plant operation cost, favoring the reaction turbine in- 
stallation, becomes $56,780, instead of the originally given value 
of $50,180. 
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Reduction of Turbine Runner 
Vibration and Noise 


A résumé of hydraulic turbine vibration experience in the U. S. Army Engineer 
Division, Southwestern, Corps of Engineers, hydroelectric stations and methods used 
in reducing the vibration and noise to a more satisfactory level. These hydroelectric 
stations are scattered over a rather large area extending from Southern Missourt to 
Central Texas. All of the hydroelectric installations are pari of multiple pur pose projects 
primarily designed for flood control. This provides an operating head variation of 
nearly 50 per cent of the turbine design head at some projects. There are 18 units 
presently in operation, all of the Francis type, with design heads varying from 59 to 190 


F. TAYLOR 


Chief, Power Section, 
Construction-Operations 

Division, Southwestern Division, 
Corps of Engineers, Dallas, Texas. 


ft. 


aa or objectionable vibration and noise have 
been experienced in 12 of the 18 hydraulic turbines installed in 4 
of the 7 hydroelectric stations in the U. S. Army Engineer 
Division, Southwestern, Corps of Engineers. As a rule, if a 
turbine does not have an objectionable vibration and noise level 
to the “feel’’ and “ears,’”’ it is considered satisfactory, and no 
further attention is paid to these characteristics. If the vibration 
and noise level are objectionable, then some effort is usually made 
to improve the operating conditions. The lack of a well-defined 
method of measuring the turbine vibration and noise level for use 
in field acceptance tests adds to the complication. Contract 
specifications for the purchase of hydraulic turbines are, in gen- 
eral, very detailed in their requirements with the exception of the 
vibration and noise level of the turbine. 

In the spring of 1944 two hydroelectric stations, Denison Dam 
on the Red River near Denison, Texas, and Norfork Dam on the 
North Fork of the White River near Norfork, Ark., were com- 
pleted with initial installation of one unit each. Due to low 
water conditions at the time, operations were on a strict emer- 
gency status until the following spring. During this period the 
Norfork turbine appeared to have a smooth flow of power and 
to be practically free of any objectionable vibration or noise. 
The Denison turbine appeared to have a smooth flow of power, 
but was very rough and noisy from about 50 to 75 per cent gate 
opening. Since both of the units were operating well below the 
bottom of their respective power pools, the turbine operating 
characteristics were not considered too seriously at the time. 
With the increase in the lake water elevations in March, 1945, 
and subsequent normal operation of the two units, it was ap- 
parent that the Norfork turbine had very desirable vibration and 
noise characteristics, while the Denison unit had very unde- 
sirable vibration and noise characteristics as well as a 1.5 to 2.5 
MW load swing at 50 to 60 per cent gate opening. 

Four additional hydroelectric stations were completed in the 
Southwestern Division in 1952-1953, and in three stations the 
turbines had undesirable vibration and noise levels. 

The following is a station-by-station résumé of operating ex- 
periences with turbine vibration and noise, work required to de- 
termine the source of the trouble, and corrective measures taken 


Contributed by the Prime Movers Subcommittee of the Hydrau- 
lie Division and presented at the Gas Turbine Power and Hy- 
draulic Conference, Houston, Texas, March 6-9, 1960, of THe Amert- 
caN Society oF MECHANICAL ENGINEERS. 

Note: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those of 
the Society. Manuscript received at ASME Headquarters, Novem- 


ber 23, 1959. Paper No. 60—Hyd-. 
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in an effort to reduce the vibration and noise. In order to mini- 
mize the unit outage time, practically all corrective work on the 
turbines was accomplished during periods of major maintenance 
from 1945-1957. 


Denison Dam 

2-35 MW units; hp, 56,000; design head, 102.5 ft; 
operating head, 85-135 ft. 

The Denison turbine vibration and noise were extremely 
severe. There was a heavy oscillating rumble at all heads of 90 
to 135 ft with wicket gate openings of 55 to 75 per cent. The 
vibration would actually aerate the lubricating oil in the turbine 
bearing oil reservoir to the point that the low oil alarm switch 
float would sink. (The float was set to be stable under normal 
operation with counterweights but was not excessively buoyant. ) 
The oscillating rumble would go in and out of phase with the 
generator’s 60-cycle hum, as when synchronizing generators. 
Vibration tests made with an RCA Sound Analyzer and crystal 
pickup indicated a rather high amplitude vibration of the turbine 
bearing and head cover near 120 cycles and at several other fre- 
quencies that were not harmonics of 120 cycles. 

In the early stages of the turbine vibration work at Denison 
Dam a great deal of effort and time was expended in making 
vibration studies to determine the frequencies and amplitude 
of vibrations on various parts of the head cover, turbine bearing, 
and other parts outside the turbine. Due to the complexity of 
frequencies and low amplitude of vibration, little helpful informa- 
tion was gained. In general, when vibration and noise were re- 
duced to a level that was not objectionable to the feel and 
ears and harmonic vibrations in cover plates, powerhouse 
windows, hand rail, ete., were sharply reduced, a definite drop in 
amplitude of some frequencies would be indicated while other 
frequencies would show little or no change in amplitude. 

As noted, the Denison turbine vibration was extremely severe, 
and it may be of interest to know the technique used in determin- 
ing the source of vibration and noise to be the turbine buckets. 
This was done by wedging oak struts or blocks between the 
buckets at about midpoint between the runner hub and the 
shroud ring near the discharge edge of the buckets. The oak 
struts would stay in place long enough to make one or two test 
runs from 0 to 100 per cent gate. The struts stopped the oscillat- 
ing roar completely and changed the vibration and noise at the 
draft tube doors from a running-gravel-type sound to the swish- 
ing sound of running water. This did not determine the cause of 
the vibration and noise, but indicated that bucket vibration or 
bucket trailing edge flutter was the principal cause of the unde- 
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sirable vibration and noise. This also made it possible to devise 
methods of reducing the vibration and noise to a satisfactory 
level. 

Seven methods have been used in efforts to reduce vibration 
and noise in the Denison Dam turbines to date. 


Unit No. 1 


1 25 to 50 cfm compressed air injected into the scroll case at 
the low pressure Winter Kennedy piezometer tap. Results were 
very effective at all heads and gate openings. 25 to 50 cfm com- 
pressed air injected at the top net head piezometer tap on the 
turbine side was found to be very effective around normal heads 
and all gate openings. Air pressure used was just high enough to 
force the air into the scroll case. Abnormally high heads require 
the higher volume of air to obtain the same results. 

2 Sharpening the entrance edge of the runner buckets was 
very effective around normal heads at some wicket gate openings. 
The sharp entrance edge was accomplished by welding !/2-in. 
key stock, on edge, to the entrance edge of the buckets, filling 
voids with weld metal, and grinding smooth (Fig. 1). This caused 
some cavitation on the back of the buckets at the entrance edges, 
requiring repair with 18-8 stainless steel inlays (Figs. 3 and 4). 

3 3-in. extensions were added to the trailing edge of the 
buckets from hub to shroud ring, using a shaped cross section 
from the bucket edge thickness to approximately 1/, in. at the 
discharge edge. This was very effective at some gate openings 
when operating around design head, but appeared to amplify the 
vibration and noise at the best gate openings of 70 to 80 per cent. 

4 Following removal of the 3-in. extensions from the trailing 
edge of the buckets, the bucket edges were ground flat or at 90 
deg to the bucket surfaces at the discharge edge (Fig. 2b). No 
improvement was noted at any head or gate position. 

5 Struts were installed between the buckets using 18-8 stain- 
less steel with a hydrofoil-shaped cross section (Fig. 5). This was 
very effective at all heads experienced and at all wicket gate 


Removed fo reduce 
cow lation 


key stock 


Fillet weld 


Pressure side 


Fig. 1 Turbine bucket modification; reshaping leading edge of buckets 


side 


a. Origina/ SAape 
b. 90° angle 
c. 45° angle 
d. Half huller-nose 
Fig. 2 Turbine bucket modification; reshaping trailing edge of buckets 
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Fig. 4 Denison turbine runner bucke? entrance edge showing cavitation 
damage 


openings; however, the struts caused cavitation ‘‘downstream”’ 
from the struts on the back of the buckets (Fig. 6). Following 
approximately 1 year’s service the struts began to fail, breaking 
off near the bucket surfaces. Replacement of one to three struts 
was made annually for 3 or 4 years; then lost struts were left 
out as long as adjacent struts were in place. This provided 
bracing between groups of three buckets. 


Unit No. 2 


6 Arunner bucket bracing ring was installed using a 17/,-in. 
mild steel bar of round cross section welded to the bucket edges 
with extra heavy welds. This was very effective at all heads ex- 
perienced and at all wicket gate openings. Following approxi- 
mately 1 year’s service the ring began to fail by cracking apart, 1 
to 3 in. from the bucket welds. One to three failures were found 
at each annual inspection until the mild steel ring was replaced 
with an 18-8 stainless steel ring. The first annual inspection re- 
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4 a 
Fig. 3 Denison turbine runner bucket entrance edge sharpened to knife 
: edge 


heads the bracing struts, in addition to the modified bucket 
discharge edges, provided better operation. 


Whitney Dam 
2-15 MW units; hp, 20,700; design head, 91.5 ft; rpm, 128.6; operat- 
ing head, 88 to 149 ft. 

The Whitney Dam turbines had a rather heavy vibration and 
noise level during the initial starting operations at below-normal 
operating heads. The condition became worse with increasing 
head, and bucket vibration was suspected. Tests were made 
with wooden blocks driven between the buckets about midway 
between the runner hub and the shroud ring at the discharge edge 
of the buckets. The results of the tests indicated that the vibra- 
tion and noise could be reduced to a satisfactory level by bracing 
the buckets. As in the results of the turbine No. 1 tests at 
Denison Dam, the vibration and noise at the draft tube doors 
were changed from the flowing-gravel-type sound to the swishing 
sound of flowing water. 

25 to 50 cfm compressed air injected into the scroll case through 
the low-pressure Winter Kennedy tap reduced the vibration and 
noise to a satisfactory level. 

A wrought-iron ring was installed on the turbine No. 1 runner 
similar to the one on Denison turbine No. 2; and, in addition, 


Fig. 6 Denison turbine runner bucket showing cavitation damage 
repairs below struts 


Fig. 7 Denison turbine runner bucket bracing ring 


vealed five places where this ring had parted. The 18-8 stainless 
steel ring was replaced with a wrought iron ring (Fig. 7). This 
ring gave good service for approximately 2 years; then it began 
to crack apart near the bucket welds, and rather severe cavitation 
occurred on the downstream side of the ring. Sections were 
finally removed from the ring, leaving three buckets tied together 
in each group (Fig. 8). This apparently provided ample bracing 
because operation of the turbine continued to be satisfactory; 
however, the remaining sections of ring showed additional cracks 
at the next annual inspection and were removed. 

7 Following removal of the wrought iron bracing ring, the 
trailing edges of the turbine buckets were reshaped to a half- 
bullet-nosed cross section, with a 45-deg angle as shown in Fig. 
2(d), extending from the fillet at the hub to the fillet at the shroud 
ring. This bas given excellent results at all gate openings except 
at abnormally high heads, when it is necessary to inject com- 
pressed air into the scroll case through the Winter Kennedy low 
pressure piezometer tap in order to maintain a satisfactory vibra- 
tion and noise level. In 1957 the unit No. 1 turbine bucket dis- 
charge edges were reshaped with the 45 deg angle as shown in 
Fig. 2(c). There was no apparent difference in the operation of 
the two turbines around normal heads. At excessively high Fig. 8 Denison turbine runner bucket bracing ring secti 
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t Fig. 5 Denison turbine runner bucket bracing struts 
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Denison Dam - 1 September 1959 


Unit Number 
Gross Hours Operation 
Gross Hours Generation 
Hours Generation Since Turbine Bucket Edge Mod. 
MWH Gross Generation 
MWH Generation Since Turbine Bucket Edge Mod. 
Unit Capacity Factor Since Turbine Bucket 
Edge Mod. 


Whitney Dem - 1 September 


No. 1 
22,210 
13,001 
13,001 

153,807 

153,507 


Unit Number 
Gross Hours Operation 
Gross Hours Generation 
Hours Generation Since Turbine Bucket Edge Mod. 
MWH Gross Generation 
MWH Generation Since Turbine Bucket Edge Mod. 
Unit Capacity Factor Since Turbine Bucket 
Edge Mod. 1% 


Fort Gibson Dam - 1 September 1959 


No. & 
31,716 
17,741 


8,587 
181,785 


No. 3 
32,196 
18,041 


No. 1 No. 2 
25,829 21,655 
14,196 13,515 


5,694 5,914 8, 366 
140,924 140,700 180,412 


61,435 64,278 89,564 
9% 8% 


Unit Number 

Gross Hours Operation 

Gross Hours Generation 

Hours Generation Since 
Turbine Bucket Edge Mod. 

MWH Gross Generation 

MWH Generation Since 
Turbine Bucket Edge Mod. 

Unit Load Factor Since 
Turbine Bucket Edge Mod. 


92,683 


Bull Shoals Dem - 1 September 1959 


No. 3 
42,710 
2, 604 


9,398 
801, 312 


328,156 


No. 1 No. 2 
52,447 50,684 
30,162 30,434 


8,224 
9TT, 701 


505,583 


Unit Number 

Gross Hours Operation 

Gross Hours Generation 

Hours Generation Since 
Turbine Bucket Edge Mod. 

MWH Gross Generation 

MWH Generation Since 
Turbine Bucket Edge Mod. 

Init Capacity Factor Since 
Turbine Bucket Edge Mod. 75% 654 874 


Fig. 9 Unit operation data 


1,659 
960, 056 


49, 57% 


the discharge edges of the buckets were modified with the half- 
bullet-nosed shape as shown in Fig. 2(d). Turbine No. 2 bucket 
discharge edges were given the same treatment with no 
bucket bracing added. Operation of the turbines for approxi- 
mately 3 years was completely satisfactory. On the annual 
inspection following the third year of service, the bracing ring 
on turbine No. 1 had cracked apart in one place and was removed 
from the runner without apparent change in the vibration and 
noise level of the turbine. 

During the year following removal of the bracing ring, ab- 
normally high headwater elevations were experienced for the first 
time. The turbine vibration and noise returned to an objectiona- 
ble level and it was necessary to inject compressed air into the 
scroll case to maintain a satisfactory operating condition. 


Fort Gibson Dam 

4-11.25 MW units; hp, 16,000; design head, 59 ft; rpm, 100; operat- 
ing head, 56 to 88 ft. 

The Fort Gibson Dam turbines were considered very rough 
and noisy at 50 to 65 per cent gate opening in addition to having 
an uneven power flow or load swing. At gate openings of 70 to 85 
per cent the units caused an excessive amount of vibration and 
noise in the powerhouse structure. With all four units operating, 
harmonic vibrations made the windows between the generator 
room and the control room vibrate enough to see double images 
in the reflections. It was an irritating condition, causing a great 
deal of comment from the operators. A conference room adjacent 
to the control room was too noisy to use with any degree of 
success. 

In the initial operating tests, the turbine vacuum valves were 
blocked open at all gate positions. This reduced the load swing 
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to a fairly reasonable amount and reduced the powerhouse 
vibration and noise level to some extent. At the time the noise 
was not considered to be caused by bucket vibration. Due to 
low inflow conditions in the area and the pressure of putting two 
new stations into operation, little effort was made toward correct- 
ing the trouble for a period of nearly 2 years. 

Following the successful reduction of vibration and noise in 
the Whitney turbines it was decided that the vibration and noise 
in the Fort Gibson turbines and powerhouse structure might 
possibly be caused by turbine bucket vibration. Bracing the 
buckets of one of the turbines with wooden blocks indicated 
the buckets to be the source of the trouble. 

The bucket discharge edges on all four turbines were reshaped 
with a 45-deg angle as shown in Fig. 2(c). The results were com- 
pletely beyond all expectations. The vibration and noise in the 
wheel pits and at the draft tube doors were practically eliminated. 
The powerhouse structure’s vibration and noise were reduced to 
a normal level. The turbine vacuum valves were returned to 
service, and at present a !/;-in. opening, above 50 per cent gate, 
cushions the load swing to a reasonable amount. 


4-40 MW units; hp, 54,000; design head, 190 ft; rpm, 128.6; operat- 
ing head, 144 to 246 ft. 

The Bull Shoals turbines are all identical; however, unit No. 3 
had a rather severe vibration and noise level at 40 to 65 per cent 
gate opening. The noise was a rather deep, steady roar that ap- 
peared to be in synchronism with the 60-cycle hum of the gen- 
erator. The unit was operated out of this critical range most of 
the time. The vibration and noise at the draft tube door were of 
the flowing-gravel type. Above 65 per cent gate opening, the 
vibration and noise level were not considered to be excessive 
except at the draft tube door. Annual inspections had shown 
cavitation at the discharge edges of the runner buckets on all 
four turbines. The cavitated areas were not consistently in any 
one location on the bucket edges, nor on every bucket. 

In 1956 the No. 2 runner bucket discharge edges were shaped to 
the half-bullet-nosed cross section as shown in Fig. 2(d). The 
results were gratifying in reducing the noise and vibration at the 
draft tube doors. Improvement at other locations was not too 
impressive; however, the results were considered well worth the 
effort. The remaining turbines were given similar treatment 
during their scheduled maintenance outages. The vibration and 
roaring noise in turbine No. 3 were reduced to a satisfactory level. 
Excessively high heads have not been experienced since complet- 
ing the bucket modification on these turbines. 


Norfork Dam 


2-35 MW units; hp, 42,000; design head, 160.0 ft; rpm, 128.6; op- 
erating head, 140 to 208 ft. 

The Norfork turbines have operated over the complete range of 
operating head with a very satisfactory level of vibration and 
noise. Cavitation damage has occurred on some areas of most 
of the bucket discharge edges. This has been severe enough to 
require some repair work with 18-8 stainless steel inlays. Al- 
though these turbines are relatively free of objectionable vibra- 
tion and noise, it is proposed to reshape the bucket edges to the 
45 deg angle as shown in Fig. 2(c) in order to reduce the cavita- 
tion action on the discharge edges of the buckets. 


Comments 

The action of the compressed air injection in reducing the 
vibration and noise caused by bucket vibration is not understood. 
As the air is injected above the wicket gates, it is assumed that 
the air bubbles flow through the turbine in the upper part of the 
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No. 1 No. 2 
103,702 39,093 
68,404 32,219 
5,298 8,415 
1,945,139 847,420 
141,059 227, 256 
No. 2 
18,745 
11,725 
11,725 
136, 986 
136,986 
18% 
a = Bull Shoals Dam 
90% 
No. & 
38, 44S 
19,407 
2,901 
664,419 
85,796 
Tht 


water passage. Considering the air bubbles to form a cushioning 
effect in the upper half of the runner, it could be possible that the 
source of the trouble is near the runner hub. With this in mind, 
reshaping of the upper half of the bucket discharge edges might 
serve to reduce the vibration and noise to a satisfactory level. 

Cavitation damage to the bucket discharge edges has been re- 
duced to a very minor amount on all runners with the reshaped 
edges. During the 1959 annual unit inspection, it was noted that 
a hard black scale was forming on most of the bucket edges. In 
Fig. 9, the Unit Operation Data show the unit service time and 
MWH generated since reshaping the bucket edges compared to 
the total unit time and MWH generated. The capacity and load 
factors shown give a general idea of unit load conditions since re- 
shaping the bucket edges. 

The failure of the bracing ring installed on the Whitney Dam 
turbine runner indicates that the bucket vibration has not been 
completely eliminated by reshaping the bucket discharge edges. 

The return of the turbine vibration and noise at abnormally 
high heads above design head may indicate that a different con- 
figuration of the bucket discharge edges would be more effective 
for all operating conditions. 
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DISCUSSION 
|. Swiecicki* 


The author comments that the action of compressed air injec- 
tion in reducing the vibration and noise caused by bucket vibra- 
tion is not understood. An explanation can be found in the paper 
“Acoustic Method of Investigation of Cavitation in Hydraulic 
Turbines,” by S. B. Stopsky, Gidroturbostroenie, Moscow, 1957. 

8. B. Stopsky describes a simple method which allows deter- 
mining incipient cavitation, its growth and sigma break equally 
wel! on a model as on a prototype. It consists of recording the 
amplitudes of pressure fluctuations in a narrow ultrasonic band of 
frequency. Describing his method Stopsky points out the great 
effectiveness of gaseous bubbles as pressure wave absorbers in 
liquids. 

When the bucket discharge edges shed von Karman vortex 
streets, the pressure waves travel from these edges in all direc- 
tions and are reflected from the water passage walls, coming back 
and superimposing on each other in a most irregular manner. 

The water itself is a very bad damping medium; to absorb al- 
most all of the input of the oscillatory energy coming from the 
impulses generated at the bucket discharge edges, the runner, 
casing, and the adjacent structures must develop a relatively high 
amplitude of vibration. 

When gaseous bubbles are present anywhere in this zone of 
intensive pressure wave travel, they effectively absorb the 
energy of the waves. According to Stopsky, this can be ex- 
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plained first of all by the flow of heat from the bubbles to the 
surrounding water as well as by the dispersion and transfer of the 
sound wave energy to the bubbles. When a part of the energy 
input is absorbed by the bubbles, only the remaining part must 
be dissipated by vibration of the turbine and the adjacent strue- 
tures. Therefore, the amplitude of vibration decreases and there 
is less noise. 

According to this analysis, air admission and cavitation, both 
supplying bubbles, should be effective in reducing the vibration of 
the unit and the noise level. Furthermore, it can be theorized 
that it is not necessary to admit air exactly at the places where 
the trouble arises because vibration occurs not only due to the 
original waves created at the discharge edges of the buckets, but 
also in connection with the waves reflected from the steel walls. 
The wall reflecting a wave absorbs only a small part of the wave 
energy, the remaining energy travels back with the reflected 
wave which must bounce back and forth several more times 
before being completely dissipated. 

In such a case it is immaterial where the vibration damper, con- 
sisting of a certain number of bubbles, is introduced as long as 
it is located in the space where the intensity of the wave travel is 
high. 

If a given volume of air is going to be admitted, the best effect 
may be expected from its admission somewhere upstream of the 
runner. In this case it remains in the zone of intense pressure 
pulsation longer before it is washed away through the draft tube. 

According to Stopsky, the size of the air bubbles undoubtedly 
is also of importance. Unfortunately, there is no known tech- 
nique for controlling the size of bubbles. The natural frequency 
of the bubble oscillation depends on its diameter, and the bubble 
absorbs the maximum amount of pressure wave energy when its 
natural frequency coincides with the frequency of the pressure 
waves. 

It appears that the several phenomena described by the author 
confirm the conclusions from the theory described in the forego- 
ing, or can be explained by it. For instance, the sharpening of 
the entrance edges of the Denison runner buckets undoubtedly 
resulted in some local cavitation at certain heads and gate open- 
ings. Therefore, the vapor bubbles were present and, as the 
author observes, this sharpening of the entrance edge was very 
effective in reducing the noise and the turbine case vibration. 
The writer knows also that when testing models in modern high 
head laboratories the turbine often produces some high pitched 
noise at a high cavitation coefficient, sigma. At lower sigmas, but 
well above the sigma break, this noise diappears. This seems to 
indicate that (in this case) the incipient cavitation, occurring 
well above the sigma break, is sufficient to prevent audible noise. 


Author’s Closure 


The author wishes to thank Mr. Swiecicki for his discussion of 
the effect of compressed air injection in reducing the vibration 
and noise level in turbine runners. The possible phenomena 
caused by sharpening the entrance edges of the buckets and its 
relation to compressed air injection into the water passage in 
reducing bucket vibration are very interesting. With this in 
mind the sharpening of the entrance edges of the buckets does 
not appear as desirable as originally considered. 

The reshaping of the bucket trailing edges to the half-bullet- 
nose or 45-deg angle, as described in this paper, is not a cure-all 
for turbine vibration; however, experience has indicated that 
the modification can be useful to the operating engineer in 
economically reducing objectionable vibration and noise caused 
by bucket vibration to a more satisfactory operating level. 
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Experiences With Hydraulic 
Prime Mover Controls 


This paper is devoted to the prime mover control aspects of some problems encountered 
by the U. S. Bureau of Reclamation with hydrogenerating units. The problems include 
governor instability, remote control of governor dashpot compensation, and control 
associated with overspeed following load rejections. 


First symptoms were often com- 


plicated by contributions from voltage regulators, relays, oil pumps, or tie-line control 


apparatus. 


Instability 


0... of the seemingly largest problems paradoxically 
developed with one of the smallest units, a 1200-kw unit installed 
at Angostura Dam in South Dakota. This Francis turbine- 
driven unit operates with a hydraulic head of 89 to 122 ft, and is 
supplied from a 54-in-diameter outlet tube 600 ft long, which also 
serves as a penstock in this case. Needless to say, with such a 
small unit, refinements such as a surge tank or automatic pres- 
sure relief valve, had been omitted for economic reasons. 

Initial performance of the unit seemed quite normal in the in- 
tended manner of operation in parallel with a fairly large system. 
Performance for synchronizing had been sufficiently stable, by 
virtue of the usual needle valve on the governor compensating 
dashpot, and the small amount of energy in the water column 
which has low velocity for such small turbine gate opening. 

To engineers familiar with the subject, it will not be surprising 
that the governor became extremely unstable when line relay ac- 
tion isolated it from the main system but left it serving a small 
load representing about 80 per cent of rated capacity. The part 
that may be a little surprising is that, after this condition was 
demonstrated, it was still impossible to stabilize the governor 
with the full range of the compensating adjustments incorpo- 
rated in the standard design governor. The speed versus time re- 
lation in Fig. 1 illustrates the degree of instability encountered. 

At this point, a number of things combined to give this simple 
problem of instability the appearance of more than ordinary com- 
plexity. For one, when the instability proved to be stubborn, the 
fact that less than the usual proportion of WR? had been provided 
was regarded with more concern. For another, an erecting engi- 
neer, upon being unable to stabilize the governor with the maxi- 
mum compensation adjustments afforded by the standard design 
governor, recommended that additional WR? be provided (which 
would be necessary if minimum frequency band width had to be 
adhered to). The governor manufacturer considered reducing 
ball head sensitivity, but decided to provide more compensation 
with a special dashpot in order to keep the rest of the mechanism 
standard. On field trial of this special dashpot, the governor 
was reported to be more unstable than ever (due to an undetected 
leak in the dashpot cylinder). Circumstances at this stage were 
sufficient that increasing the WR?, providing a surge tank, or 
providing a synchronous bypass, all began to be reconsidered. 


Contributed by the Prime Movers Subcommittee of the Hydraulic 
Division and presented at the Gas Turbine Power & Hydraulic 
Conference, Houston, Texas, March 6-9, 1960, of THe AMERICAN 
Society oF MECHANICAL ENGINEERS. 

Note: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those of 
the Society. Manuscript received at ASME Headquarters, December 
9,1959. Paper No. 60—Hyd-9. 
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Altention ts directed to the components or effects principally responsible. 


The practical need for a more immediate solution, not to men- 
tion a less expensive one, released the problem to field engineering. 
It was thus determined that a 120 per cent change of floating 
lever ratio was sufficient to stabilize the governor. However, 
the combination of 40 per cent stiffer ball head springs and 63 
per cent increase of the compensation ratio was considered pref- 
erable, as it afforded proportional response from the ball head 
over a wider speed range. Final characteristics are illustrated in 
Fig. 2. 

It is worth cautioning that increasing the ratio of the floating 
lever beyond the normal range may also necessitate altering the 
limit bypasses of the dashpot secondary piston in order to still 
permit gate travel at the maximum selected rate at or beyond 
limits of ball head control. 

Fig. 3 is included because it illustrates the effect of water 
starting time, an especially interesting characteristic of a unit with 
a long penstock (without surge tank). There is an optimum rate 
for opening the turbine gates. If they are opened faster than the 
optimum rate, power output momentarily drops rather than in- 
creases (because opening the gates faster than the quantity of 
flow in the penstock increases simply causes a decrease of entrance 
velocity into the water wheel). 

Some factors of significance demonstrated by this small but 
interesting case are: 

a. That instability must result when it is attempted to reg- 
ulate to a finer degree than over-all characteristics of the prime 
mover and load are suited. 

b. Some applications can exceed the range of compensation 
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GATE POSITION 


Initial instability of small hydrounit illustrated by speed versus 
gate position relation (plotted from observations for comparisons with 
oscillograms in Figs. 2 and 3). Amount of gate travel is limited by period 
of oscillation. Speed oscillation becomes progressively worse. 
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Fig. 2 Oscillogram of final performance of hydrounit of Fig. 1 after 
stabilizing alterations. Connection to main power system is opened at 
left of oscillogram constituting a rejection from 92 to 62 per cent of rated 
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Fig. 3 Test with isolated load, showing effect of gate travel too fast for 
hydraulic characteristics. Speed has been depressed by gate limit which 
is lifted at A’. Gate opening increases from A’ to B’ at maximum rete 
allowed by relay valve limit. Speed drops over corresponding interval. 


normally provided in standard design governors. In such cases, 
specifications’ emphasis on stability rather than on sensitivity 
may avoid difficulty and delay on initial operation. 


Control of Compensation 


Quite another source of difficulty with prime mover control has 
developed through the occasional need to alter compensating 
adjustments of the governors to suit various operating conditions. 

Stability requirements can be analyzed, but in general the 
dashpot leak-off time constant may decrease as the connected 
inertia increases. In a manner of speaking the dashpot must 
compensate for any lack of stability or inertia of the load. Thus, 
for each change of the system which changes the connected in- 
ertia, a new adjustment of the governor dashpot would be needed 
if optimum damping were to be continually maintained. Un- 
fortunately, the changes which affect inertia of the system usually 
take place somewhere remote from the governing units. There 
is no direct means available to sense the new adjustment needed. 
This is the crux of the problem. 

Fortunately, considerable latitude in the adjustment can be 


1N. G. Dennis, “Water Turbine Governors,” Proceedings, I. 
Mech. E., vol. 1b, series b, 1953, pp. 379-408. 
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tolerated under most conditions. The difference between 
optimum and average adjustment can scarcely be detected 
under normal system and load conditions. Hence, an ad- 
justment suitable for average normal system conditions may 
suffice as long as there is appreciable system and load con- 
nected to the unit. However, when the unit is isolated from the 
system, its own inertia being widely different than that of the 
system to which it may have been synchronized, requires a con- 
siderably different compensating adjustment in order to be stable. 
Hence, the simplest and most common automatic compensating 
adjustment is activated by gate position. If the gates are not 
above speed-no-load position, the governor is made very stable 
for ease in synchronizing. If the gates are above speed-no-load 
position, the governor is permitted to be appropriately active for 
system operation, gate position being used as a measure of the 
compensation needed, even though it is an indirect measure. 

Obviously, gate position is not always a reliable criterion of the 
compensation needed. Under abnormal! system conditions, a 
unit may become isolated from the main system, but still have 
appreciable isolated load. When this happens, the unit may be 
undercompensated and unstable. Worse still, two or more units 
may become isolated from the stabilizing effect of the main system 
and become seriously unstable. 

Still another need for adjustment beyond capabilities of the 
gate position criterion is occasioned by spinning reserve operation, 
in which a unit is kept motoring on the system but with gates 
closed and usually with the water depressed from the wheel to 
reduce losses. The objective is to have the capacity of the unit 
immediately available in case of loss of other generation from the 
system. To satisfy this the unit must respond immediately to 
drop of speed or change of the speed level adjustment. How- 
ever, if the governor is overcompensated, until the gates are 
above speed-no-load position it can respond only very slowly, 
losing many precious seconds in the emergency it is supposed to 
relieve. 

This shortcoming gave rise to our first use many years ago 
of remote control of governor compensation. This was by means 
of a solenoid rigged to replace the mechanical arm for lifting the 
dashpot bypass rod. The solenoid could be energized by the 
operator at will regardless of gate position. It permitted main- 
taining the governor in readiness for prompt action when the unit 
was placed in spinning reserve. It served this purpose quite well, 
so well in fact, that this means of control has also been employed 
for other conditions and with some success too, but with occasional 
embarrassment also. 

The use of solenoid operated bypasses on the governor com- 
pensating dashpots has also been used to speed up the response 
of governors to the impulses from load, frequency, or tie-line 
control. It has not been uncommon that units which have been 
compensated to be stable under the worst operating condi- 
tion, respond so slowly to the control impulses that excessive 
travel of the speed level control accumulates before the needed 
load change comes about, resulting in some slow hunting of the 
tie-line control system and unsatisfactory control. However, 
completely opening the bypass usually yields good response to 
the control impulses, permits holding tie-line loading within much 
closer limits, and still does not cause signs of instability if the 
connected system is large. The embarrassment does not come 
until some automatic switching operation, such as fault clearing, 
separates the unit or several units together from most of the 
inertia of the system. Instability may result at this rather 
critical time and, unless the operator can manually trip the so- 
lenoid bypass within a few seconds, the governors may oscillate 
until the servomotor oil supply is exhausted. The units then shut 
down and restoration of service must be delayed until the oil 
supply is restored and the units restarted and synchronized. 

Various means have been used in special cases to trip the sole- 
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noid bypass automatically, such as with line circuit breaker trip- 
ping or with tie-line control tripout. Other possibilities include 
overfrequency, underfrequency, and low governor oil pressure, 
but these are also indirect measures of system inertia and there- 
fore are not infallible indications of need for this stabilizing func- 
tion. Allowing the bypass to be opened only momentarily by 
each tie-line control impulse, closing again after a short time 
delay, offers the most assurance against accidental instability 
when employing the solenoid bypass but this arrangement has 
not been popular because of the added complication of the time 
delay relay and associated circuitry. 

Need for these automatic means of tripping the dashpot bypass 
solenoids is aggravated by the all too common practice of operat- 
ing with less compensation than actually yields any benefit. 
There being no simple, obvious means of determining the opti- 
mum compensation, the solenoids are often allowed to remove 
compensation altogether. No adverse effect is apparent under 
normal conditions but instability becomes more common and 
more severe under system emergencies. The tendency toward 
operating with the dashpot bypass wide open is also contributed 
to by the illusion that system frequency regulation is always 
aided by this practice. Actually, the increased activity of the gate 
positioning mechanism is no proof of benefit. The change of 
power output of a hydrounit lags behind the initiating disturb- 
ance because of inertia of the water column and gates. When the 
lag is 90 deg, the regulating efficiency becomes zero even though 
the regulating effort may still be high. For greater lags, the 
contribution to system frequency regulation is detrimental 
rather than helpful. Hence while the solenoid controlled dash- 
pot bypass is a useful device, its existence has been responsible for 
undercompensation not only on some of our plants but on those of 
some interconnections. 

One direct first step which would help relieve this problem 
would be to allow compensation to be reduced no more than to 
the optimum point, as may be determined by the frequency re- 
sponse technique.? This would avoid some of the cases of insta- 
bility against which the artificial safeguards are otherwise 
needed, but there still would remain those cases where additional 
compensation would be needed due to system sectionalizing. A 
compromise which is usually acceptable is to operate with no less 
than the latter degree of compensation. Stability under emer- 
gencies is insured and the sacrifice of response for system regula- 
tion is usually negligible. 


Load Rejection Overspeed 


The overspeed characteristics also, of hydrogenerating units, 
have given rise to some interesting control features. 

The overspeed following rejection of full load from our Francis 
wheel-type hydrounits with reasonably good hydraulic charac- 
teristics reaches about 140 per cent of normal speed. Full run- 
away speed on our units of this type is about 185 to 190°per cent 
of normal speed. 

On units of this category our experience has shown simple 
overspeed backup protection to be the most satisfactory. This 
backup protection is comprised of a speed switch set to operate at 
about 145 per cent of normal speed, that is, about 10 per cent 
above the speed which develops following rejection of full load. 
In the event this speed switch operates, it trips the shutdown 
solenoid, a hand reset device which requires attention of the 
operator before operation can be resumed. 

This scheme relies upon the ball head to effect prompt closure 
of the turbine gates under any normal condition. Following a 
rejection the ball head naturally returns the unit close to normal 

2V. Oja, ‘‘Frequency-Response Methods Applied to the Study of 
Turbine Regulation in the Swedish Power System,” Trans. ASME, 
vol. 76, 1954, pp. 1325-1333. 
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speed where it remains in readiness to be resynchronized by the 
operator. Only in the event of failure or slow operation of the 
ball head would the speed rise to the speed switch setting. This 
indication of abnormality or malfunction is reason enough to 
bring about shutdown such that checking and resetting by the 
operator is necessary to resume operation. 

In some cases, auxiliary control comprising speed-no-load 
solenoids operated by speed switches set to operate at about 115 
per cent of normal speed have been used. This seems to have 
been with the idea that gate closure would be brought about 
quicker or more reliably. Actually at speeds above about 110 
per cent normal the governor ball head initiates closure of the 
turbine gates at their maximum rate and functioning of the 115 
per cent speed switch and speed-no-load solenoid is rather super- 
fluous, the jeb having been already accomplished by the ball 
head as illustrated in Fig. 4. This prompt functioning of the 
ball head is by virtue of the normal dashpot ports which bypass 
the delaying action of the compensating dashpot when the speed 
has departed from normal setting by 5 to 10 per cent, depending 
upon the linkage ratios used. 
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Fig. 4 Rejection of 55,000 kw from Unit 5, Davis Power Plant, Arizona. 


Speed-gate position relation and speed switch operation plotted from 
oscillograms. Note that governor ball head starts gates closing at maxi- 
mum rate before operction of speed switch and speed-no-load solenoid. 


In this discussion of experience it seems appropriate to remark 
that we have found the governor ball head, in general, to be highly 
reliable. Primary dependence upon it so far as possible has 
therefore seemed quite logical. Speed switches have been less 
precise. 

On units with less ideal hydraulic characteristics and with 
correspondingly more complicated overspeed characteristics, the 
policy of placing reliance upon the ball head as far as possible, and 
providing essentially backup protection is still applied, but 
application of the policy is by no means as simple. Units with 
less ideal speed characteristics are those with long penstocks, 
with no pressure relief valve or synchronous bypass and with no 
surge tank. Still more complicated for speed control are the 
installations where two or more units may be served by one long 
penstock. Installations of this sort are not numerous, but have 
been utilized where economically advantageous. Their charac- 
teristics are varied, but in general the magnitude of the overspeed 
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following a load rejection on these units is not a reliable criterion 
of whether the condition is abnormal or not. If the penstock is 
long enough, the transient overspeed following a rejection may 
be as high as the sustained runaway speed. Where more than one 
unit is served from the same penstock, the transient overspeed 
depends upon whether one or more units have rejected load. In 
the latter case, the transient overspeed is usually very nearly as 
high as the sustained runaway speed. 

We have tried the same type of control as applied to the more 
ideal units with the idea of merely accepting the nuisance of re- 
starting the units any time the overspeed exceeds the nominal 
speed switch setting, but this has usually proved unsatisfactory. 
Consequently, for such units, duration of overspeed has been 
employed as the criterion to distinguish whether equipment has 
malfunctioned and whether the shutdown solenoid should be 
tripped. This is accomplished by interposing a time delay relay 
between the speed switch and the shutdown solenoid. If the 
speed switch resets before the relay has timed out, the shut- 
down solenoid is not tripped. 

Choice of setting for the time delay relay can become quite 
involved. Theoretical calculation is quite possible but is 
lengthy and expensive; experimental determination is also ex- 
pensive in terms of disturbance, manpower, and recording equip- 
ment. The approximate method shown in Fig. 5 has been satis- 
factory and may be of practical interest. The shape of the speed 
versus time curve is taken from the test usually made for other 
purposes on a single unit. The magnitude of overspeed when two 
(or more) units are rejected is taken from design calculations. 
The duration of overspeed to the value at which the speed switch 
resets is then determined as shown in the figure. 


Conclusions 


On small, low cost hydropower units, specification emphasis 
should be placed upon stability rather than upon governor sensi- 
tivity. 

Solenoid control of governor dashpot bypasses is a useful de- 
vice, but one which may lead to instability under system emer- 
gencies unless adjustments are made for this condition. De- 
pendence upon manual closing of the bypasses under emergencies 
is undesirable. 

For load rejections, dependence upon the governor ball heads 
as the normal gate closure actuating element with speed switch 
actuated shutdown solenoid for backup protection affords a 
simple and reliable means of control. For long penstock applica- 
tions duration of overspeed is a useful criterion for effecting 
emergency shutdown. 

The technology evolved and consolidated, mainly within the 
past decade, affords analysis of performance of hydrogenerating 
unit control, but the wide range of operating conditions from large 
interconnections to isolated or local area operation still affords 
interesting problems and occasions for exercise of ingenuity. 
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Fig. 5 Determination cf duration of overspeed following rejection of 
two 18,950-kva units (served by common penstock 1500 feet long). 
Available data were 1-unit rejection test and calculated maximum speed 
for 2-unit rejection. Alcova Powerplant, Wyoming. 
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DISCUSSION 
C. L. Avery® 


The authors are to be commended for presenting a very frank 
and informative paper. 

NEMA publication No. HT-4 recommends certain minimum 
values of mechanical! inertia for hydro-electric generators as de- 
termined by the unit hydraulic condition. These recommenda- 
tions are based on the use of standard governor designs. If the 
recommended values are not met, difficulties, which require 
special treatment, such as experienced at the Angostura Plant, 
may be expected. 

It seems evident that the solenoid control of the governor dash- 
pot provides a satisfactory method for quickly disarming the 
dashpot bypass which is convenient for automatic operation by 
suitable control relays and switches. In this respect it is an im- 
provement over the mechanical bypass. Also, the solenoid 
operated dashpot bypass is adjustable to suit system conditions 
whereas the mechanical bypass requires substitution of different 
bypass rods for change in adjustment. Both are subject to the 
same limitations respecting stable control. 

Overspeed switches are designed as back-up protection. If the 
governor does not hold the speed to normal overspeed values 
then the cause for malfunctioning should be determined. Other- 
wise, it is desirable and often imperative to get the unit back 
on-line as soon as possible. The use of a partial shutdown solenoid 
for overspeed protection seems superfluous. It may be used in 
the starting sequence in which case it is set above speed-no-load. 
It cannot be used for both purposes concurrently. 


¢ Engineer, Woodward Governor Company, Rockford, Ill. Mem. 
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Authors’ Closure 

Mr. Avery’s comments are all pertinent and valuable. He has 
wisely called attention to the fairly recent (1958) NEMA publica- 
tion which conveniently indicates acceptable ranges of hydraulic 
and mechanical characteristics as well as the range of adjustment 
in the standard design governors. This should simplify imple- 
menting the first conclusion in the paper. 
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There is no question that the solenoid control of the dashpot 
is much more convenient and flexible than the mechanical bypass. 
It is the latter advantage that gives rise to the comments of 
the paper. Limiting the amount of bypassing to no more than is 
really needed can reduce the amount of reliance which otherwise 
must be placed upon means of disarming the device in emer- 
gencies. 
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Effects of Brine Dispersion in the 
Allen Salt-Velocity Method 


This paper presents information concerning the effect of brine dispersion on the Allen 
salt-velocity method of measuring flow. The data are presented from a variety of 
sources, theory, laboratory investigation, and field tests. 


Theory and tests show that 


brine injected into a conduit disperses as the square root of the distance traveled. Asa 
result of this relationship, salt-velocity curves occupy a smaller percentage of the passage 
time as the test length increases. This results in smaller test errors for long test sections. 
Also indicated in this paper are the results of further tests which show the minimum 


Introduction 


HE salt-velocity method devised by C. M. Allen and 
reported in 1923 [1]! has been used both in the laboratory and the 
field for a number of years, and during this period a considerable 
amount of experience has been amassed. It is proposed in this 
paper to discuss primarily the effects of dispersion in relatively 
straight pipes. 

Professor Allen appreciated from an experimental point of view 
that turbulence in mixing the salt solution was quite necessary for 
the successful application of the method. Not only was it im- 
possible to measure by means of the salt-velocity method in 
laminar flow but it also became very difficult to measure by this 
method in small pipes. Furthermore, two of the important 
effects of dispersion were known very early in his experience. 
The Searsburg test in 1923 provided one of the longest test sec- 
tions either in terms of pipe diameters or in feet that has ever 
been used. This particular test section was provided with 
several sets of electrodes and the growth of the curve was easily 
established at this time. This test showed beyond a doubt that 
the brine slug grew in length as the square root of the distance 
traveled along the conduit. 

The second piece of evidence was the velocity traverse which 
was reported in his original paper and which was made in the 
40-in. pipeline at the Alden Hydraulic Laboratory. In this case 
it was demonstrated that the time to the center of gravity of a salt 
curve measured by a point electrode was very nearly the same at 
the center line of the pipe as at the wall of the pipe. This 
demonstrated quite vividly that a particle of salt must be mi- 
grating radially so that it was part of the time in fast-moving 
water and part of the time in slow-moving water if the transit 
time was to become as uniform as that indicated in the tests. 
This experience led to the simile of a salt slug moving down the 
pipe in much the same fashion as a flock of sheep moving along a 
country road with good grazing on either side. 


1 Numbers in brackets designate References at end of paper. 
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quarters, August 1, 1960. Paper No. 60—WA-138. 


limits for length of test sections with equipment commonly used. 


Theory of Dispersion 

In 1954, G. I. Taylor [2] presented a discussion of the dispersion 
of matter in turbulent flow where he determined the following 
relationship: 


S? = 437aX (=) (1) 


Using the diameter instead of the radius and dividing both sides 
of the equation by the diameter (D)?: 


S? xX 

S 


It has been shown that V,/U is equal to theY. f/8, Reference [3], 
where f is the friction factor in the Darcy-Weisbach-Fanning 
formula. 

This allows a simple evaluation of the friction velocity and 
therefore of the remaining factors of the equation. 

It is seen that in (2) the dispersion equation can be arranged so 
that the terms are dimensionless. Lengths are expressed as pipe 
diameters. This has the further advantage of facilitating com- 
parison between various test sections. The effect of friction has 
been expressed as a “friction velocity’’ and then as a ratio to the 
pipe velocity. It will be seen in Fig. 1 that this factor does not 
vary widely. Finally and most significantly it is seen that the 
dispersion of the salt solution varies as the square root of the 
distance traveled. 


or 


Field Checks of Dispersion 


Table 1 compares the predictions of Sir Geoffrey Taylor with 
what has been found in various test conditions. Pertinent data 
from a large number of field tests have been used to compute the 
constant in the G. I. Taylor equation (2). In this tabulation an 
“average” value of friction factor was used based upon the ex- 
perience curve of Russell [4]. It will be found that the value of 
V_/U as shown in Fig. 1 used in the determination can vary only 


Nomenclature 


= radius of pipe, ft = friction factor in Darcy-Weisbach- U 
C = coefficient of dispersion in equation 
(2) S= 
= 218 theoretically in conduit, ft 
= pipe or conduit diameter, ft t = time, sec 
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Fanning formula y 
longitudinal length of a slug of brine ° 


average pipe velocity, fps 


friction velocity, fps 


distance traveled by a slug of brine 
from point of introduction, ft 
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Fig. 1 Relation between U/V, and log Reynolds number. A—smooth 
pipes, B—pipes of average roughness, C—rough pipes. 


over rather narrow limits and that other factors are causing much 
larger variations so it is believed that this assumption is of sec- 
ondary importance. 

Referring to Table 1 it will be seen that there are some values of 
growth which are reasonably close to the theoretical value of 218, 
such as Hetch Hetchy, Davis Bridge, Serra No. 1, and Sorocaba 
No. 4. But the deviations away from the theoretical value are 
more the rule than the exception. It is seen that Searsburg and 
Tepuxtepec, for instance, have much larger rates of dispersion 
while there is a great number of tests with much smaller coeffi- 
cients of dispersion. 

The reason for these divergencies are for the most part obvious 
from the physical conditions under which the tests were made. 
G. I. Taylor commented that the effect of bends in a pipeline 
considerably increased the rate of dispersion and the resulting 
growth of the curve. The 8-ft conduit at Searsburg was a wood- 
stave pipe which followed the valley of the river so that there 
were very few tangent sections in the plan but a number of long- 
radius bends. The test section at Tepuxtepec similarly had a 
bend of 10 deg between pop valves and the electrodes. 

On the other hand, there is a large number of tests where the 
coefficient of dispersion was approximately one quarter of the 
expected theoretical value. It will be noted that most of these 
tests are associated with the note “turbulator.” 

A turbulator was a device installed in the pipe to create a 
large degree of turbulence immediately ahead of the test section. 
So far as is known it was first applied on the Metropolitan Water 
District tests in Southern California. In this case large amounts 
of salt solution of great concentration had to be introduced 
into the pipe to obtain suitable curves with the recording am- 
meters in use at that time. Density effects were manifest and the 
turbulator was devised and installed to provide sufficient energy 
for mixing so that the dense salt solution would not damp out 
the turbulence already existing in the pipeline. This was reported 
in a previous paper [5]. Since that time the turbulator has been 
used frequently, particularly for tests performed for the Bureau 
of Reclamation. The rate of growth of the curves is slower 
where the turbulator is used resulting in a more compact curve. 
This is of course borne out by the small coefficients of dispersion 
found in the table. 

The turbulator frequently used by the Bureau is illustrated 
in’ Fig. 2. It is constructed with several rigid plates across the 
pipe, causing a contraction effect of the order of 50 per cent of 
the pipe area. The resulting pressure difference causes a re- 
arrangement of the velocity distribution in the expansion zone 
behind the bars so that the normal pipe distribution is completely 
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destroyed and a practically uniform velocity distribution is set 
up across the pipe for a short distance. Due to the head loss 
across the turbulator the turbulence just downstream is con- 
siderably higher. It is seen, therefore, that there are two 
effects of the turbulator (and for that matter any disturbing 
influence in the pipe). One is the effect on the velocity distri- 
bution, in this case being made much more uniform than normal 
pipe flow. The other effect is the increased amount of energy 
that is available for mixing the brine or maintaining a more 
uniform distribution across the pipe cross section. Both of 
these effects work together with a turbulator, resulting in a lower 
coefficient of dispersion as seen in Table 1. 

As normally installed the turbulator is located upstream of the 
pop valves a sufficient distance so there will be no separated flow 
at the salt-injection station. The primary effect is to secure, 
by means of the increased turbulence, a more rapid mixing of the 
brine solution following its injection at the introduction station. 
The effect of this increased turbulence upon the test results was 
checked at Flat Iron Pump Turbine Unit No. 3 test when the 
unit was being tested generating. In this case the turbulator 
was used throughout the test but it was removed for the last 
three tests of the series. Comparing the last three tests with 
the three tests immediately preceding, which were made with very 
nearly the same flow and pressure conditions, showed no change 
in the discharge measurements as indicated in the following 
tabulation: 


Fiat Iron Unit No. 3 generating 
Net head, Discharge, Average discharge, 
ft cfs cfs 
402.4 
399.2 
400.5 


401.4 
401.1 
400.5 


400.7 with turbulator 


401.1 without turbulator 


Fig. 2 Turbulator in 12-ft pipe 
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20 
> 
S 
10 
Run 
20 
21 
22 
23 221.9 
24 221.7 ee 
25 221.7 
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PIPE DIAMETER 
FT. 


NAME 


TEPUXTEPEC 6.9 


SEARSBURG 8.0 


DAVIS BRIDGE 9.0 


DAVIS BRIDGE TUNNEL 14.44 


PARAHYBA NO. 2 9 


SOROCABA NO. 4 5.74 


GRAND COULEE 12.00 


FLAT IRON NO. 2 


IRON MOUNTAIN 10 

CORNER BROOK 4.5 
AGUA FRIA NO. 1 6.0 
Cc. BD. HOWE 4.0 
CHANDLER NO. | 10.0 
ROZA 12.0 
GLENDO 21.0 
HIWASSEE 18.0 


HETCH HETCHY 8.69 


SERRA NO. | 5.08 
Ss FLAT IRON NO. 3 8.00 
TURBULATOR 

FLAT IRON NO. 3 8.00 


NO TURBULATOR 


196 / 1961 


13.35 449 196 6.14 6.788 24.2 38500 3090 302 10° BEND BETWEEN 


10.6 57.5 70 13.4 7.127 25.0 4900 1090 112 


12.6 257 133 4.81 6.682 23.8 14400 1475 232 VENTURI METER 


14.81 


4.97 


7.60 80 27.5 2.28 6.358 23.1 75% 360 49 INJECTION 5D 
6.54 108 45.0 2.62 6.418 23.2 2025 648 72 INJECTION 5D 


6.58 4 19.9 1.76 6.246 22.8 3% 160 37 INJECTION 12 D FROM 


7.86 40 26.2 5.25 6.720 24.0 686 400 4) TURBULATOR 
130 53.1 2822 1300 52 
13.48 63 39.8 10.78 7.032 24.8 1584 758 52 TURBULAT OR AND 
94.3 8885 2550 86 BEND 


4.0 55 56.4 
145 101 
17.7 53.4 46.6 21.3 7.328 25.5 2172 960 58 
156. 64.6 4174 2810 38 CONVERGING 
3.99 80.1 85.5 2.31 6.34 23.1 7310 696 242 10" BEND IN 
187.6 152 23090 1630 327 SECTION 


12.91 158. 89.7 4.3 6.642 23.8 8055 802 240 BEND AFTER 
671 172.6 29770 310 208 POP VALVES 
8.00 4 31.8 4.27 6.6304 23.8 1010 320 75 TURBULATOR 
140 64.2 4120 1120 88 
8.00 “0 35.9 4.277 6.6304 23.8 1288 320 % NO TURBULATOR 
140 73.0 530 1120 114 


Assume V = 1.5x 1075 sq. ft/sec (46° water) 


Experimental check of dispersion constant 


FT. /SEC. 
FT. 


$2 
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x 
REMARKS 


ob 


PIPE VELOCITY 
POP TO ELECTRODE 
CURVE LENGTH 


LOG Rn 


POPS AND ELECTRODES 


LARGE RADIUS BENDS 


STRAIGHT 


STRAIGHT 


CONVERGING 


1175 252 54900 6750 193 STRAIGHT 
BEND 


TURBULATOR 


TURBULATOR 


TURBULATOR 


BELOW BEND 


BELOW BEND 


PUMP MANIFOLD 


6.0 6.799 24.2 3182 1150 67 TURBULATOR 
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SLUG LENGTH / PIPE DIAMETER - S/D 


20 


30 


J DISTANCE -X Fi PIPE DIAMETER -D 
Fig. 3 Rate of spread of salt solution without turbulator 


Variation of Spread With Distance 

The dispersion equation indicated that a straight-line relation- 
ship should be found when the rate of spread or dispersion 
(S/D) was plotted against the square root of the distance traveled 
VV X/D, both distances being expressed in terms of pipe diameter. 

In Fig. 3 the rather long test sections without a turbulator 
have been plotted. A reasonable agreement with a simple 
straight-line relationship was found. As noted in the foregoing, 
the basic rate of dispersion is modified by the pipe-wall roughness, 
(V,/U) and by the presence of bends or changes in section. 

In Fig. 4 the initial portion of Fig. 3 is plotted to a larger scale 
so that the test results on the shorter test sections can be seen 
more clearly as well as including the test data for turbulators. 
A first glance at Fig. 4 seems to indicate a rather large dispersion 
in the results of the spread of the salt solution. Referring to 
Table 2 and Fig. 3, it will be seen that the tests made in conduits 
without a turbulator tend to start from zero and go off at a 
constant rate as indicated by the straight-line plot at a constant 
slope. Viewing the large number of tests included where a 
turbylator has been used it will be seen that this group of tests 
starts off on one slope and continues with a reasonable agreement 
for a distance of approximately 25 diameters (when referred to 
the over-all distance between the turbulator and the electrode). 
In this region a few tests indicate that the rate of dispersion 
seems to change and the slope of the new curve follows parallel 
to the original curve arrived at without artificial turbulence 
being added. In other words the effect of the turbulator seems 
to die out at approximately 25 diameters and the rate of growth 
then returns to the more normal rate as determined by the friction 
factor for the straight pipe. 
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The test section at C. D. Howe was located approximately 13 
diameters below a manifold for several centrifugal pumps and 
it is seen in this case that the effect of turbulence of the pumps is 
much the same as for a turbulator located in the pipe only to a 
somewhat lesser degree. In this case also it is seen that the 
second curve which was located 75 diameters beyond the pop 
valves has a coefficient of dispersion of 194 which is not far from 
the theoretical value of 218. 

In the Agua Fria tests the injection station was located 5 
diameters below a 90-deg elbow. Flow table tests had shown 
that for this arrangement the velocity distribution was approxi- 
mately that for a normal pipe. Reference to Fig. 4 and to the 
test value of coefficient of dispersion shows that the first curve is 
relatively compact whereas the second curve is more nearly 
normal for a straight pipe. It is seen therefore that the effect 
of an elbow is twofold in that the first effect is to create a large 
difference in velocity from one side of the pipe to the other in the 
plane of the elbow. This effect results in lengthening the curve 
artificially as noted by Taylor and seen in these tests. The 
other effect that we see in the Agua Fria tests is to create an 
increased turbulence downstream from the elbow which tends 
toward a more compact initial curve in this case with a lower 
coefficient of dispersion. 


Electrode Location 


It is essential that the curves of the brine passing two sets of 
electrodes should not overlap for an accurate test. The rate of 
growth of the brine slug determines, therefore, how closely the 
electrode stations can be located ina conduit. Table 2 and Fig. 5 
have been prepared for design purposes establishing the length 
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Fig. 4 Rate of spread for short test section with and without turbulator 


TO CENTER OF GRAVITY- % 


DISTANCE FROM START OF CURVE 
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Fig. 5 Position of center of gravity in an electrode curve 


of the brine slugs and the center-of-gravity position in the curve far apart that the curves do not overlap. In using Fig. 5, there- 
for various test arrangements. fore, the length of the curve after the center of gravity would be 

Continuing the practical aspects of design, Fig. 5 shows the taken from the lower limit for the first curve and the length of 
position of the center of gravity with respect to the length of the the curve before the center of gravity would be taken from the 
curve for the various conditions of mixing that have been tabu- upper limit for the second curve. This is equivalent to assuming 
ated. This information can only be used to establish limits. the worst possible conditions so as to insure a safe prediction for 
In application the two electrodes must be located sufficiently performance. 
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f 
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m C. D. HOWE 
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Fig. 6 Velocities in 40-in. line as measured by pitometer and salt-velocity 
method. A—pitometer traverse, B—salt-velocity, traverse-point elec- 
trode 22D from injection station, C—salt-velocity, traverse-point electrode 
104D from injection station. 


Radial Dispersion 

The effect of radial dispersion of the brine was brought out 
very well in Professor- Allen’s paper [6]. These tests dem- 
on:trated that salt may be introduced at any point in a diameter 
and it will relatively quickly find its way to all portions of the 
pipe. 

Other test information showing the effect of radial dispersion 
is presented in Fig. 6. There the data of three traverse tests in 
the 40-in. line at the Alden Hydraulic Laboratory are plotted. 
The first curve is the instantaneous velocity diagram as measured 
by a pitometer. It is seen that the pipe factor, which is the 
ratio of the mean velocity to the center line, is 81 per cent. 
The second traverse curve was measured by a small point 
electrode which was set up to traverse the 40-in. pipe at a distance 
of 22 diameters from the injection station which consisted 
of four pop valves. In this case the pipe factor was found to be 
89.5 per cent. The third curve is another point-electrode 
traverse taken at a distance of 104 diameters. In this case the 
pipe factor was found to be 98 per cent. It is found that the 
pipe factor varies approximately as the square root of the distance 
between the introduction section and the electrode cross section. 

It is also seen that the particles do not occupy the same stream 
tube throughout their journey between an introduction station 


’ and a recording or electrode station. There is actual material 


transport from one point in the radius to another in the turbulence 
process. This means that if the motion of one particle is studied 
it may be found in one instance at the mean velocity point, at 
some distance downstream it may be traveling at the center of 
the pipe at the maximum possible velocity, and again farther 
downstream it may be found in the very slow water near the 
wall of the pipe. The longer the test section, the more op- 
portunity there is for the migration of each individual particle 
to all parts of the conduit and the more nearly the passage time 
of each individual point will approximate the average passage 
time for the group. This is indicated by the salt-velocity 
traverse curve at 104 diameters indicated previously, where the 
difference between the center-line salt velocity and the wall 
velocity, both as measured by the transit time from the injection 
station, is found to be only 4 per cent. 

As a consequence, the design of electrodes in very long test 
sections can be quite approximate indeed, since all parts of the 
traverse curve that is recorded by the salt-velocity method 
differ from one another by very small amounts, and slight un- 
balances in the conductivity weighting of the electrode will have 
little difference upon the result. Of course in short test sections 
this is not true, and in all instances electrodes are designed to 
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give as closely as possible equal conductivity weighting to equal 
areas in the conduit. 

As near as can be determined this factor of radial dispersion 
was not appreciated by deHaller [7]. In paragraph 13 of his 
paper where the flow cenditions are described preparatory to 
setting up the equation of distribution he states, ‘Turbulent 
flow in a pressure conduit differs from laminar flow in that the 
path of the fluid particle is no longer a straight line but an 
irregular line oscillating back and forth in all directions without 
in the meantime going astray very far.’’ Again, in paragraph 
14 the statement is made, ‘This front (of the brine slug) has a 
tendency to deform because of the difference of velocity between 
the adjacent points on the front. At some instant, ¢, the dis- 
placement of P can be represented by the equation: X = 

Both of these statements are contrary to what has been expressed 
in the foregoing, relating to the observed motion of particles 
in turbulent flow. In the first place the particles are actually 
transported from one part of the pipe to the other and do not 
remain in the same general location. In the second place 
the front is not deformed at the rate of velocity multiplied by 
time but rather some velocity multiplied by the square root 
of time as indicated by observations recorded. 

The assumption that the curves would lengthen linearly with 
the distance travelled in the pipe applies to laminar flow. How- 
ever, there is complete agreement that the salt-velocity method of 
testing cannot be applied to laminar-flow conditions. There- 
fore the analysis presented by deHaller would probably be true 
for laminar conditions but certainly not for turbulent condi- 
tions, especially with radial dispersion which obviously takes 
place. 

One of the effects that the combined turbulence, both lon- 
gitudinal and radial, has on the method is that the total length 
of the electrode curve occupies a less and less percentage of the 
total length of the passage time between the introduction station 
and the final station as the total length of the test section in- 
creases. It is fair to assume that electrode errors can only occur 
during the passage or registration of the curve. The error 
therefore might be expressed as some constant multiplied by the 
length of the curve, which, in turn, can be expressed as a constant 
multiplied by the square root of the length of the test section. 
The error would therefore become the error in time divided by 
the passage time which is proportional to the length of the test 
section. Whence the expression for error becomes €./L/L. 
This would indicate that a long test section would have a smaller 
error. 

Table 3, showing interchecks by the salt-velocity method, 
gives a number of instances of long test sections where this 
thinking is borne out by experience. The test of Unit No. 5 
at the Serra Development of the Sao Paulo Tramway Light and 
Power Company is a good case in point. This test was reported 
by O. H. Dodkin and reference to his paper may be made for 
further details [8]. 

Referring to the Serra No. 5 test section between the pop 
valves and the first electrode (54 diameters), it is seen that the 
casual and systematic errors were both relatively large varying 
from +1 to —1.41 per cent for individual runs. Extending the 
test length to 234 pipe diameters by using the third electrode 
reduced the errors to +0.22 and —0.33 per cent. 

In this connection it is known from laboratory tests that the 
pop-valve introduction has a consistent error differing in nature 
from the electrode error. It is known that it can be minimized 
by keeping the distribution piping size small and the test section 
long. At the present time the pop-valve-introduction error 
cannot be eliminated. 

For this reason testing between electrodes is preferred for 
precise work since experience shows that the errors do cancel 
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Table 3 Interchecks with Allen salt-velocity method 


FT. 


NO. of RUNS 
NO. of SHOTS 


TEST LENGTH 


x 


SEARSBURG 


38 


351 
281 
632 


DAVIS BRIDGE 


DAVIS BRIDGE 
TUNNEL 281 


12900 


44) 
746 - 


TEPUX TEPEC 


281. 3m 
84.2m 
223. 5m 
365. 5m 
139. 2m 
142.0m 


SERRA NO, 5 


WITHOUT SIGN 
% 


WITH SIGN 
% 
TEST SECTION 


x/D 
MIN, DIFF. 
AVE. DIFF 
AVE. DIFF. 


89 
91 


= Pop Valve Injection Station - 4 valves except Davis Bridge Tunnel and Searsburg 


where | valve only was used. 


= Electrode Station - Double Bellied Electrodes in every case. 


reasonably well. Reference to the results of tests between the 
electrodes at Serra No. 5 will bear this out so far as consistency 
is concerned at least. 

It would seem that since the effect of a turbulator is to provide 
more compact electrode curves it also should result in reduced 
electrode errors for a given length of test section. The tests 
cited in the foregoing at Flat Iron No. 3 show no change in the 
measured flow whether or not the turbulator was used. Some 
laboratory tests with and without turbulators have been made in 
a 12-in. line where again no improvement—or loss of accuracy— 
was found. It is hoped to repeat this series of tests on turbulators 
in the 40-in. line to see if any improvement in operating conditions 
will result there. 

Dr. Martin Mason in his doctoral thesis showed that the 
instantaneous distribution of salt solution in a pipe could be 
represented by a bell-shaped or Gauss curve. He further showed 
that the recording of such a curve which was undergoing lon- 
gitudinal growth would come out as a distorted Gauss curve with 
the center of gravity displaced along the time axis. From this 
analysis the displacement of the center of gravity of the distorted 
curve would be the same amount whether a long curve or a 
short curve was being recorded. Furthermore it signifies that 
the electrode errors due to curve growth should essentially 
cancel. 

Table 3 shows that experimentally this seems to be the case. 
The difference in timing between a relatively short test section 
and a very long test section gives the same answer within the 
limits of experimental accuracy. < 


Short Test Sections 

In all of the theoretical discussions a Gauss type distribution 
curve is assumed which has been distorted by the method of 
recording the salt concentration. As the test section becomes 
shorter the rate of growth of the curve becomes relatively much 
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more rapid and the distortion of the curve becomes extreme. 
For this reason a series of tests was made in the 40-in. line at the 
Alden Hydraulic Laboratory to determine how short the test 
section could become before serious errors were introduced when 
measuring between electrodes. 

Four */,-in. pop valves were used to introduce the salt solution 
into the 40-in. line as shown in Fig. 7. Two 9-bar grid electrodes 
were used for each electrode station constructed as indicated in 


Fig. 7 Injection station in 40-in. line 
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Fig. 8 Nine-bar grid electrode in 40-in. line 


Fig. 8. The pop-valve introduction station remained fixed and 
the location of the electrodes within the pipe was changed from 
test to test. The recording ammeter used for this work had a 
range from zero to 11/2 amp. 

In the first series of tests the second electrode was located at 
a reasonably large distance and the first electrode was moved 
progressively closer to the pop valves until serious errors were 
produced. The results are shown in Table 4 and Fig. 9. It was 
found that the electrode error remained unchanged until a 
distance of about three pipe diameters was reached. The error 
then increased rather rapidly in a direction to measure too small 
a quantity of water when measuring between the two electrodes, 
reaching 1 per cent at a distance of 2.3D. 

The cause for the error in this case was quite obvious from the 
curves themselves inasmuch as they became very jagged showing 
that the salt solution introduced by the pop valves had not 
reached a reasonably uniform mixture in the pipeline. 

The first electrode was then fixed at a distance of 4.5 pipe di- 
ameters where its error was small and the second electrode was 
moved progressively closer to the first electrode. The results 
are shown in Table 5 and Fig. 10. The error reached 1 per cent 
with an electrode separation of 1.2 pipe diameters. Reference 
to the recorded curves showed that there was no change in error 
until the curves started to overlap, that is to say, the second curve 
started to come through before the first curve was entirely cleaned 
up. As soon as this overlapping occurred the error started 
to increase as shown in the plotted results. 


Conclusions 


From a consideration of theory and experience given in the 
paper it is believed: 

1 The length of the curve varies as the square root of the dis- 
tance from the introduction station to the electrode. The length 
of the curve is also affected to a lesser degree by pipe roughness, 
bends, and initial turbulence. 

2 The G. I. Taylor equation shows that test sections may be 
compared on a model basis when their linear dimensions are ex- 
pressed in pipe diameters. 

3 Computing the passage time between centers of gravity of 
electrode curves provides the most reliable short test section with 
present knowledge. 
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ERROR IN Q-% 


DISTANCE BETWEEN POP VALVES AND THE 
FIRST ELECTRODE IN PIPE DIAMETERS 


Fig. 9 Effect of length of mixing section 


Table 4 Effect of mixing section length pop valves to first electrode, salt- 
velocity method 


Distance 
P-1 P-2 1-2 Q Error Ave. Dist. 
Error P-1 
Fr. Fr. Ft. c.f.s. % % Pipe D 
15.0 64.0 49.0 10.71 -0.09 
21.59 0.14 
32.91 0.33 +0.12 4.51 
14,93 28.30 13.37 22.04 -0.6 
34.27 +1.2 
9.78 -0.9 -0.1 4.49 
8.64 17.46 8.82 23.26 -1.2 
35.3% +1.2 
10.09 +0.7 +0.2 2.59 
8.28 17.46 9.18 20.58 ka 
41.0 2.4 
8.90 -1.1 +1.0 2.49 
4.00 17.46 13.46 22.85 a 
3.87 10.6 
9.93 7.0 +9.2 1.20 
6.72 17.46 10.54 35.80 0.6 
20.34 0.7 
8.81 27 1.3 2.08 


P = Pop Valve Station 
1 = First Electrode Station 
2 = Second Electrode Station 


Error Based on Discharge Measured by Calibrated Venturi Meter 


4 Theerrors both due to the introduction and the electrode sec- 
tion tend to become less as the length of the test section increases. 

5 With four pop valves in a cylindrical pipe it is reeommended 
that four diameters of pipe be allowed before the first electrode 
section. It is probable that with more pop valves this distance 
can be reduced. 

6 When the test section is taken between two electrode sec- 
tions the recorded salt-velocity curves must not overlap. 

7 While theoretically a turbulator should give reduced errors, 
the fact has not been established experimentally. 
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Fig. 10 Effect of length of test section between electrodes 


Table 5 Effect of length of test section between electrodes, salt-velocity 
method 


Distances 


Test Section 
Pipe Diam. 


Pop Valve Section 
1 First Electrode Section 
2 = Second Electrode Section 


Error in $.V.M. Discharged compared to calibrated Venturi Meter 
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DISCUSSION 
E. B. Strowger? 


Professor Hooper’s paper dealing with the effects of brine dis- 
persion on the salt-velocity method is a welcome contribution to 
the Hydraulic Division of ASME. It goes a long way in defining 
the conditions required for an accurate application of the method. 

The writer wishes to discuss the use of the method in the upper 
part of a penstock having a well-rounded entrance leading from 
a large reservoir. In this case, a combination of laminar and 
turbulent flow would take place. The work done by Donald Ross 
at the Ordnance Research Laboratory, in connection with the 
working section flow studies for the 48-in. water tunnel at The 
Pennsylvania State College,* showed as per Fig. 11 that both types 
of flow exist for perhaps 15 to 20 diameters. Also, tests made by 
Stanley Barnes at Stanford University,‘ using small smooth tubes, 
showed that turbulent flow throughout the entire pipe cross 
section did not take place until about 30 diameters from the 
intake had been reached. 

This is explained as follows: When a fluid flows from a reser- 
voir through a rounded entrance into a long tube or pipe, the 
velocity distribution is nearly uniform at the entrance. Wall 
friction soon begins to act on the boundary layer retarding it. The 
shearing force is gradually transmitted radially as the fluid flows 
down the pipe, resulting in an increased thickness of the boundary 
layer until finally the entire flow is incorporated in the bound- 
ary layer. 


20d + to completely turbulent flow 


Uniform 
‘entral Loyer 


x \ TURBULENT FLOW 
TRANSITION 


LAMINAR FLOW 


Fig. 11 illustrates the growth of the boundary layer which is 
laminar up to the transition point and which then becomes tur- 
bulent. If Reynolds’ number is greater than 10°, it may be as- 
sumed that the boundary layer acts as though it had been turbu- 
lent from the start, the effective point of zero thickness being at 
a distance Xo from the start of the straight pipe. At high Reyn- 
olds’ numbers the distance X») may be negative. The velocity in 
the central core is uniform and is turbulent in the annual area 
around the core. 

Of course, the installation of racks and intake piers will cause 
the length of pipe to the point of completely turbulent flow to de- 
crease. On the other hand, many pipes today are of welded con- 
struction and are extremely smooth with a Mannings’ friction 
coefficient of 0.011 and even 0.010. At any rate, the point where 
the flow is turbulent in the whole pipe cross section may be at a 
considerable distance downstream from a good bellmouth intake. 

With this situation a more or less uniform distribution of veloc- 
ity would exist at and near the intake and nonuniform distribu- 
tion would take place some 15 to 20 diameters downstream where 


2 System Project Engineer, Hydro Stations, Niagara Mohawk 
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3 Donald Ross, “Turbulent Flow in the Entrance Region of a 
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4Stanley Barnes, ‘‘An Experimental Investigation of the Flow 
Conditions in the Inlet Length of aSmooth Round Tube With a Bell- 
mouthed Entrance,”’ graduate thesis, Department of Civil Engineer- 
ing, Stanford University, Stanford, Calif. 
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20 
1.0 
P-1 P-2 T-2 Q Error Ave. 
Ft. Fr. Fr. cfs % Error 
15.0 64.00 49.0 10.71 -.09 
21.59 +.14 
32.91 +.33 +0.12 14.70 2 
14.93 28.30 13.37 22.04 -0.6 -0 
34.27 +1.2 __ 
9.78 -0.9 -0.1 4.01 
14.93 20.59 5.66 21.15 0.8 
34.22 1.5 Bi 
9.82 -1.8 +0.2 1.70 
4 
14.93 16.67 1.74 21.00 15.3 
34.20 3.9 
10.41 4.2 4.0 0.52 
Fig. 11 
14.25 17.46 3.21 9.05 0.8 
22.30 2.0 
35.70 2.0 1.6 0.96 


turbulent flow is fully developed. In between these two points, 
the stream line flow in the core would not be accompanied by 
turbulent mixing but outside the core turbulent mixing would take 
place. 

Dr. P. deHaller, in his analysis of the salt-velocity method, 
first made an evaluation of the error due to the deformation of the 
salt cloud as it passes from the injection station to the first elec- 
trode. He made the assumption that it does not continue to de- 
form as it travels from the first to the second electrode station and 
further, that the effect of dispersion is negligible. He concluded 
that for n = 7 in the distribution equation given below, the error 
would amount to —4.8 per cent. He then analyzed the method 
for various electrodes symmetrical with respect to the axis of the 
conduit but not exposing equal pipe areas for corresponding areas 
of electrode. Thirdly, he made a theoretical study of the effect of 
dispersion which he claimed would not modify his conclusions as 
to the accuracy of the method. 

Professor Hooper’s experience, however, indicates that the 
effect of dispersion is to reduce the theoretical error of —4.8 per 
cent quite appreciably. 

In his analysis, deHaller did not attempt to evaluate the error 
which might result from a different distribution at the two elec- 
trode stations. He used the following general equation for the 
distribution of velocity over a diameter of the pipe: 


in which v represents the velocity at any radius r, up the velocity 
at the pipe axis, n an exponent which defines the degree of non- 
uniformity of velocity across a pipe diameter, and FR the radius 
of the pipe. He assumed n to be the same for the two electrode 
stations but pointed out that n may be different at these two 
points, say, ~ at the upper station if it is near the entrance and 7 
at the downstream station if it is well downstream. He stated 
that this constituted another error. The following is a quotation 
from one of his articles: 

“Directly at the pipe inlet, n has very high values corresponding 
to an approximately uniform distribution of the velocity, with a 
thin boundary layer of friction at the walls. (n = o implies 
perfect uniformity.) After passing through a stretch equivalent 
to some 20 to 30 pipe diameters, the boundary layer has gradually 
become so thick that also the central parts are affected by wall 
friction and the typical turbulent distribution of velocity sets in. 
In between, all conceivable transitions exist, so that with shorter 
pipes it may happen that the first electrode is in a zone of velocity 
distribution, whose exponent n is considerably larger than at the 
point where the second electrode is situated. This is, of course, 
another source of error.”’ 

To apply the method to pipes with bellmouth intakes and with 
a total length less than, say, 15 diameters, it would seem that the 
use of a turbulator to create fully turbulent flow would be 
advisable. 


G. I. Taylor® 


Professor Hooper compares measurements of the dispersion of 
brine in turbulent pipe flow with a theoretical value given in one 
of my papers.’ He finds, as I found that bends in the pipe give 
rise to increase in the length of the plug of brine beyond that cal- 
culated, using my assumption that the pipe is straight and is 


5 P. de Haller, ‘Theoretical Considerations of the Allen Method of 
Measuring Quantities of Water,” Escher Wyss News, January-March, 
1930. 

‘Cavendish Laboratory, University of Cambridge, Cambridge, 
England. 

7“The Dispersion of Matter in Turbulent Flow Through a Pipe,” 
Proceedings, Royal Society, vol. 223, series A, 1954. 
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long enough for a steady regime to establish itself. This last assump- 
tion was legitimate in discussing dispersion in the pipes used for 
transporting oil over long distances, the problem to which I ap- 
plied the results of calculation, but is not legitimate when applied 
to cases where extra turbulence is produced by a turbulator. 
The effect of a turbulator is shown experimentally to be to shorten 
the plug. This is, in fact, exactly the effect predicted theoretically 
in a paper of mine® to which Professor Hooper did not refer. Since 
the experiments described in that paper were made with fluid in 
laminar flow it is perhaps natural that attention was not called 
to the contrast between the results in the two papers.” 

In both papers the dispersion was taken to be due to variation 
in the longitudinal convection with radius, the dispersion being 
proportional to Umax — U where Umax is the maximum velocity 
and U the mean. In both papers the effect of transverse dis- 
persion was shown to be to reduce the convective longitudinal 
dispersion and the analysis described the reaction beween the two 
effects. When the transverse dispersion is due to a cause like 
molecular diffusion which is not related to the flow velocity, an 
increase in transverse dispersion necessarily leads to a decrease in 
plug length. The effect of increasing the resistance of a turbulator 
is analogous to that of increasing molecular diffusability and 
therefore should shorten the plug.’ Increase in the hydraulic re- 
sistance of a turbulator may also be expected to have another 
effect if the resistance is uniformly distributed over the cross 
section. It will reduce the value of (Umax — U)/U. This will 
also reduce plug length. 

The analysis of my 1953 paper,® therefore, shows why a turbula- 
tor would be expected to decrease plug length. The question re- 
mains why my 1954 paper’ predicts an increase in plug length 
when the turbulence is increased, not by a turbulator, but by in- 
crease in the roughness of the wall, ie., by increase in V,/U. 
The analysis of my 1954 paper’ is based on the same principles 
as my 1953 paper,’ but in order to calculate the transverse dis- 
persibility I had to make use of Reynolds analogy between the 
transfer of matter and of momentum in turbulent flow. The re- 
sult of the analysis was that, for given values of U and the pipe 
diameter, an increase in roughness of the pipe, which causes an 
increase in V,/U, causes an increase in longitudinal dispersion, 
i.e., an increase in plug length. The reason for this is that increase 
in V,/U is associated with an increase in (Umax — U)/U and, 
therefore, with a direct increase in longitudinal convecture dis- 
persion and this is greater than the decrease which results from 
the increase in lateral turbulent dispersion. 


Author's Closure 


The discussions by Mr. Strowger and Professor Taylor are 
greatly appreciated and an effort will be made to discuss the 
points raised particularly by Mr. Strowger. 

First, there is described the mechanics of flow at a bellmouth 
entrance to a straight pipe, the flow entering without previous 
turbulence of any sort. It is believed that there exists no differ- 
ence in our viewpoints on this condition, namely, that where 
there is no turbulence as in laminar flow the salt velocity method 
will not work. The only difference in viewpoint is the next step 
as to whether or not this situation is actually found in the testing 
of water wheels. It is the author’s belief that the practical situa- 
tion and the laboratory situation are different in one important 
respect, namely, the turbulence given to the flow by racks and 
supporting beams in any normal power plant. By way of illustra- 
tion it may be pointed out that it is customary procedure in 
towing tanks never to take the first run in the morning until a 
preliminary run has been made to stir up the water. Obviously 
this disturbance is in the nature of turbulence and as compared 


§ “Dispersion of Soluble Matter in Solvent Flowing Slowly Through 
a Tube,” Proceedings, Royal Society, vol. 219, series A, 1953, pp. 186, 
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to the energy of turbulence in penstocks it must be very modest 
indeed. However, this turbulence allows the performance of 
the ship models which are being tested to approach a standard 
operating condition and it is for this reason that it is introduced. 

A number of tests have been made in the past with various 
intakes but the most accessible is that given in Reference [8], 
“Field Checks of the Salt Velocity Method,’’ by Dodkin, who 
includes a test at Parahyba of Unit No. 2. This unit was tested 
with two different arrangements of the salt velocity equipment 
and check tested with the Gibson Method. 

In both cases a converging intake connected a canal to the 
cylindrical penstock. The brine valves were located at the 
emergency stop log slot about 4 meters downstream of the racks. 
The electrodes were located in the cylindrical penstock and just 
below a converging elbow. The performance measurements 
by the salt velocity method in a converging intake agreed very 
well with the measurements made by the Gibson method. This 
would indicate to the author that while the arrangements indi- 
cated by Mr. Strowger are theoretically possible, it is a 
difficult situation to meet in practice. 

At a later point Mr. Strowger makes reference to the paper by 
Dr. P. deHaller [7]. In this case several conclusions are drawn 
from an analysis made by deHaller. In all fairness it should be 
pointed out that these conclusions were reached in 1930. Stating 
the case briefly Dr. deHaller made an empirical approximation 
of a set of salt-velocity test curves and based his salt-cloud-growth 
analysis on this situation. The electrode analysis referred to was 
based on the resistance of the fluid contained between the bars 
as for a solid whereas in the practical case the resistance of elec- 
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trodes in an electrolyte is a three-dimensional pattern and is not 
as simple as assumed. Finally, the dispersion of the brine was 
assumed to extend only axially and to vary linearly with distance. 
As noted in this paper this last assumption in particular is be- 
lieved to be the most important and leads to the same sort of a 
conclusion that would be reached for an analysis of laminar 
flow since no exchange of material radially was permitted. From 
this analysis fairly large errors (—4.8%) were deduced. In 
contrast Dr. Mason [9] assumed that the brine was injected into 
a cross section in a very short slug and showed that this brine 
was dispersed into Gaussian shaped curves downstream in the 
pipe as a result of the action of turbulence in the flow. The 
analysis indicated that the rate of growth varied with the one 
half power of the distance as has been shown by experience. 
Further analysis of the distortion caused by recording such a 
curve with a fixed electrode was made with the final result that 
the center of gravity of such a distorted curve was displaced a 
constant amount for a given pipe. This would indicate practical 
cancellation of errors between electrodes which is in good agree- 
ment with experimental evidence. 

Prof. G. I. Taylor in his discussion has presented a very 
clear explanation as to why the effect of a turbulator is to tem- 
porarily shorten a curve whereas the presence of elbows and sim- 
ilar dislocations in the pipe alignments results in an excessive 
elongation of the curve. The author is in complete agreement 
with the discussion of the mechanics involved. 

In conclusion both discussions are a help to a clear under- 
standing of the problems and the author is grateful for both con- 
tributions. 
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Supercharging of a Large Two-Cycle, 
Loop-Scavenged Diesel Engine 


The paper covers the mechanical and thermodynamic aspects of increasing the power 


output and thermal efficiency of an existing nonsupercharged two-cycle, loop-scavenged 
diesel engine through supercharging. The power output was increased by 40 per cent 


and the thermal efficiency by approximately 4 per cent. 


Much of the thermodynamic 


work was done on a UNIVAC computer and a brief description of the theoretical ap- 
proach and results appear in this paper. Field test results from a Nordberg two-cycle, 


Introduction 


HE ever-existing demand for increased output, im- 
proved economy, and reduced weight in power vehicles presents 


Table 1 Engine data 


Non- 
super- uper- 
charged charged 

21'/2 X 31 2117/2 X 31 
Engine speed, rpm.................. 240 240 
Bhp (continuous), hp................ 3050 4300 
Bhp (2-hr overload), hp............. 3355 4730 
Bmep (continuous), psl.............. 74.5 105 
Bmep (2-hr overload), psi............ 81.95 115.5 
Number of cylinders................. 6 6 
Maximum firing pressure, psig........ 850 1050 
Combined exhaust temperature, deg F 600 625 
Scavenge manifold pressure, psia .... . 18.01 29.00 
Scavenge air flow, lb/sec............. 18.4 25.1 
Auxiliary-blower shp................ 300 350 
Wristpin bushing max press, psi ..... 2430 2680 
Crankpin bearing max pressure, psi.... 1265 1742 
Crankshaft main bearing max pressure, 

Lube-oil pump, gpm................. 300 660 
Lube-oil pressure at engine, psig...... 30 45 


20 


Cooling-water pressure at engine psig. . 20 
Forged one piece 


Crankshaft construction............. 


Crankpin diameter, in............... 143/, 143/, 
Crankshaft-journal diam, in.......... 151/2 
Engine eff volumetric compression ratio 12.6 9.15 


Contributed by the Oil and Gas Power Division and presented 
at the Winter Annual Meeting, New York, N. Y., November 27- 
December 2, 1960, of THe American SocreTy oF MECHANICAL 
Enoineers. Manuscript received at ASME Headquarters, August 
12, 1960. Paper No. 60—WA-185. 


Nomenclature 


21'/:-in. X 31-in. supercharged engine are also included in this paper. 


all engine manufacturers with a problem, the successful solution 
of which is the key to their prosperity and continued growth. 

Supercharging has long been the diesel industry’s recognized 
answer to this problem and, through ceaseless efforts in research 
and development, the industry is continually providing the market 
with more powerful engines having improved efficiency and 
greater reliability. 

This paper will attempt to describe some of the problems en- 
countered and methods of solution in the process of converting a 
large nonsupercharged two-cycle, loop-scavenged diesel engine 
to supercharging, thereby increasing the power output by 40 per 
cent and improving the fuel economy by over 4 per cent. 

The engine that will be discussed in this paper is a two-cycle, 
21'/;-in. bore, 31-in. stroke, 240-rpm, 6-cylinder engine. Com- 
plete data on engine as nonsupercharged and supercharged are 
given in Table 1. 

This engine was originally erected as a test vehicle in the plant 
of the author’s company, Milwaukee, Wis., where tests were 
conducted from April, 1957, through July, 1959. 

The engine is presently operating as part of the Municipal 
Power Plant in the Town of Hudson, Mass., where acceptance 
tests were successfully conducted in May, 1960. 


Design Changes 


Due to a number of improvements, such as stiffening of bed- 
plates, A- frames, and connecting-rod bearing box, and so forth, 
having been made over the preceding few years, it was possible 
to convert an existing engine to supercharging for test purposes 
without any changes to the existing design. 

Following a complete mathematical analysis of all relevant 
parts, it was decided that mechanically the maximum cylinder 
pressure could be raised from 850 psi to 1050 psi without any 
changes whatsoever. However, as the engine is intended for 


bhp brake horsepower 
ahp = air horsepower 
bmep = brake mean effective pressure, la 
psi P, = 
rpm = revolutions per minute AP 


cfm = cubic feet per minute 
V, = scavenging air, cfm at 29.92 in. 


Hg and 68 deg F T. 
K, = engine equivalent orifice, nondi- 
mensional Tn 
k = adiabatic exponent 1.4 
P, = atmospheric pressure, psia ie 


P, X 144 psf 
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P,, = scavenging manifold pressure, 


psia 
preturbine exhaust pressure, psia 


Neom) = turbocharger combined  ef- 


turbine back pressure, psia in. 
= P,, — P,, psi H conversion factor, in-lb/lb fuel 

preturbine temperature, deg F F = fuel rate, Ib/min 

ambient temperature (deg F +- D engine piston displacement, cu 
460) deg R in. 

scavenging manifold tempera- d 
ture (deg F + 460) deg R N camshaft speed, rpm 

preturbine exhaust temperature * constant 
(deg F + 460) deg R b constant 


ficiency, per cent 
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EFFECT OF LUBE OIL HEADER PRESSURE 
UPON WRIST PIN LUBE OIL PRESSURE 
AT 240 R.P.M. ENGINE SPEED 


| | 
SUPERCHARGED, 42 L.O.HEADER PRES 


NON-SUPE RCHARGE 


8 


PRESSURE PSIG 
8 


LUBE OIL HEADER 


10 20 30 40 50 


LUBE OIL FLOW THROUGH 
PISTONS, GPM PER CYL. 


pressure, ps / 


Table 2 Heat balance based on high value of 19350 Btu/Ib and 100 
per cent load ie} 180 270 360 


LUBE OIL PRESSURE AT WRIST PIN, PSIG 


Nonsupercharged ogwetenns DEGREE CRANK ANGLE 


Btu/ = Beu/ Fig. 2 Effect of lube-oil header pressure upon wristpin lube-oil pressure 


at 240-rpm engine speed 


Water vapor 
Intercooler, Ist.... 
Intercooler, 2nd.... 


Jacket water CONNECTING ROD 


Nore: Gain in fuel economy by supercharging: 


7250 
= 6950 = 1.0432 or 4.32 per cent 


both marine and stationary service, it was decided to increase 
the crankshaft main-journal diameter from 15!/2 to 16'/2 in., 
Table 1, so as to satisfy the requirements of The American Bureau 
of Shipping and other similar inspection agencies for any increase 
in load that could be anticipated. 

Thus, except for the foregoing change, the engine was built as 
a test unit with standard parts only. 

During the entire shop tests only two changes having any 
effect on mechanical and thermal stresses were made. 

During the initial stages of the shop test the engine was run as 
a nonsupercharged unit, and heat-balance tests, Table 2, were 
conducted. Based on these tests an increase in piston-cooling 
oil flow was decided upon. This was accomplished by increasing 
the piston-cooling outlet orifice and also by increasing the lube- 
oil header pressure from 30 to 45 psig. The total effect of these 
changes upon lube-oil flow through the pistons may be observed 
in Fig. 1. 4 


The increase in lube-oil header pressure from 30 to 45 psig was -5 


not undertaken solely to increase flow through the pistons. Pres- 
sure pickup tests at the wristpin bushings revealed that with 30 
psig lube-oil header pressure the pressure at the wristpin bushing 
was negative for approximately 30 per cent of the complete cycle, 
and that only by raising the header pressure to 40 psig and over Fig. 3 Connecting rod, old and new design 


NEW NEW OLD 
DESIGN DESIGN DESIGN DESIGN 
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could a positive pressure at the wristpin bushing be assured for 
the entire cycle, Fig. 2. 

It should be noted here that for the last 1200 hr of shop test 
whereof 800 hr of operation were at 117 bmep, no changes at all 
were made to the engine. 

Although no distress in any parts could be noted at the com- 
pletion of the shop tests, mathematical and static studies had 
indicated that the connecting rod was close to the borderline of 
safety stresswise in the smallest cross section of the eye. This 
condition was found to be most critical at extreme low loads where, 
due to the very low scavenging pressure, the inertia forces of the 
moving parts would override the compression pressure during 
part of the compression stroke. 

Accordingly, the afore-mentioned cross section was increased to 
the maximum allowable within the limitations of the existing 
piston space, thereby increasing the critical area by approximately 
50 per cent. Although not required for the present rating of the 
engine, but with future increases in power in mind, the upper radii 
and the thickness of the foot of the rod were also increased sub- 
‘stantially, Fig. 3. 

It was considered further that, should the lube-oil header pres- 
sure for some reason drop below the recommended 45 psig at a 
time when the engine could not be shut down, there would be 
very little margin of safety in the wristpins. It was therefore de- 
cided to redesign the wristpin lubrication system to safeguard 
against such an eventuality. An illustration of the old and new 
designs may be seen in Fig. 4. As will be noted, the main dif- 
ference in the two designs lies in the location at which the lube oil 
is fed te the bearing and to the piston. 

In the old design the oil is fed through a hole in the bottom of 
the bushing and forced up to the six longitudinally located oil 
grooves through a centrally located groove, and to the piston 
through the bore in the wristpin. 


NEW DESIGN 


Fig. 4 Wristpin lubrication, old and new design 
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In the new design the lube oil is led around the bushing through 
a groove centrally located in the eye of the connecting rod and is 
fed to the bushing at two points located 75 deg from the bottom 
on either side of the bushing. This bushing also has six longi- 
tudinal groovesand one centrally located feed groove. The piston- 
cooling oil in this design is fed over the top of the bushing and 
through the bore in the wristpin. 

The advantages of the new design over the old design are two- 
fold: (1) Should the lube-oil pressure for some reason drop to such 
a degree that negative pressure at the wristpin might be ex- 
perienced through part of the cycle, a suction effect resulting in 
air being sucked into the bearing may be experienced with the 
old design, whereas in the new design this not only cannot happen, 
but as an added measure of safety there is a built-in oil reservoir 
to keep the bearing lubricated during such a period of negative 
oil pressure as mentioned previously. , 

(2) By relocating the oil-feed holes in the bushing, valuable 
bearing surface was gained. 

The increase in the width of the connecting-rod eye, Fig. 3, 
also added considerably to the bearing area. 

It is interesting to note that during the shop tests at 100 per 
cent load the wristpin-bushing maximum pressure based on net 
bearing area was 3500 psi. This pressure with the new design 
has been reduced to 2680 psi which compares with 2430 psi for 
the old design as a nonsupercharged unit. Thus, while increasing 
the output by 40 per cent and the maximum cylinder pressure by 
25 per cent, the wristpin-bushing maximum unit loading in- 
creased by only 10 per cent. 

No further changes have been made to the original design. 


Thermodynamics 


From preliminary studies it was apparent that very little re- 
liable information was available on the problem of matching and 
sizing turbochargers and auxiliary blowers for supercharging of 
two-cycle engines. 

Thus, from the very beginning, it was realized that, in order to 
achieve success, a rational and definite approach to this problem 
had to be developed. 

From preliminary tests two useful values were found to have 
very good correlation. 

The engine equivalent orifice K, is constant for a given two- 
cycle engine regardless of load, temperature, and pressure condi- 
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tions. This was found to be the same for nonsupercharged and 
supercharged conditions: 


(1)? 


For this particular engine K, was established to be 1600. 

Further, based on numerous tests a straight-line equation was 
established, empirically, as a suitable definition of the relationship 
between preturbine exhaust temperature, bhp, and scavenging 
air, cfm, Fig. 5: 

t= 192.7 + 2020.5 CBP (2) 
efm cfm 

It is further known that for the successful operation of any 
given diesel engine, the following equation must be satisfied in 
order to supply sufficient air for combustion: 


Adiabatic Air Hp 


For a turbocharged diesel engine the foregoing condition may 
be satisfied by either one of the following methods. 


(a) Turbochargers only—self-sustained. 
(b) Turbochargers with partial assist of an auxiliary blower. 


In either case, the following equations may be written for the 
turbocharger turbine: 


Turbine Adiabatic Hp 


k-1 
k PoV 41.027, k 
(4) 
In the foregoing it has been assumed that the weight of fuel 


burned has increased the intake-air weight by 2 per cent. The 
turbocharger, of course, has a net output of + zero hp. 


Compressor shaft hp = turbine shaft hp (5) 
or 


Compressor adiabatic hp = turbine adiabatic hp X teomb (6) 


Depending upon engine requirements, available energy to the 
turbine and turbocharger combined efficiency, the turbine out- 
put may be greater or smaller than the required adiabatic hp, 
equation (3). 

Additional adiabatic hp + 
= required adiabatic hp (equation 3) 


—turbine adiabatic hp (equation 4) X tfeom» (7) 


Additional adiabatic hp + 
k-1 
( P, 


Additional adiabatic hp + 


1R. Birmann, ‘‘Aerothermodynamic Considerations Involved in 
‘Turbocharging Four and Two-Cycle Diesel Engines,’’ Proceedings 
of the ASME Oil and Gas Power Division, 1955. For definition of 
values, see Nomenclature. 
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k 


P, 1.02 


T, 
4 X Neomb (9) 


Equation (9) when programmed and run on a UNIVAC com- 
puter gave all desired values in their optimum conditions. 

Fig. 6 illustrates a typical plot of computed values based on 60 
per cent turbocharger combined efficiency at 105 bmep and 240 
rpm. 


PRE-TURB 


PRE- TURB 
650 F 
PRE-TURB 
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TURBOCHARGER ADIABATIC HORSEPOWER 
DEFICIENCY /CYL. WITH TURBO- 
CHARGER COMBINED EFFICIENCY OF 
60% AT 105 BMEP AND 240R.P.M. 
DEFICIENCY TO BE MADE UP BY 
AUXILIARY BLOWER 


Fig. 6 Air-horsepower requirements as related to pressure ratio and 
preturbine exhaust temperature 


The plot shows turbocharger adiabatic horsepower deficiency 
to be made up by an auxiliary blower as related to pressure ratio 
and preturbine temperature. 

The plot demonstrates particularly well the importance of ar- 
riving at the optimum pressure ratio since obviously both too 
great and too small a pressure ratio will greatly increase the re- 
quired auxiliary-blower horsepower for any given preturbine 
exhaust temperature. 

The equation may be computed for self-sustained operation for 
any given combination of exhaust temperature, pressure ratio, 
and turbocharger combined efficiency by simply considering the 
additional adiabatic horsepower requirements as zero. 

Of course the selection of preturbine exhaust temperature can- 
not be governed by the natural desire to reduce auxiliary-blower 
horsepower requirements, but by the thermal loading that the 
engine can tolerate and to what extent the scavenging air flow 
may be reduced without increasing the fuel economy. 

For this particular engine the selection of 625 F preturbine 
temperature, Table 1, was determined by the fact that the engine 
was sold as a dual-fuel engine and the scavenging air flow was 
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based on requirements for detonation-free operation on gas and 
not for requirements on oil operation. 

The reliability of equation (9) is best demonstrated by com- 
paring the auxiliary-blower adiabatic horsepower requirements as 
calculated for operation at the power plant at Hudson, Mass., and 
the value measured during the acceptance test of the engine. 

The values are given at 100 per cent load and pressures and 
temperatures as specified in Table 1: 


It should be noted that the calculated value was based on 
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Fig. 7 Fuel economy as related to turbocharger combined efficiency 


ATMOSPHERE 
EXHAUST LOUVERS 
SILENCER, 
= FILTER 


625°F AUXILIARY 
BLOWER 
= INTERCOOLER 


3 


FOR 
AUX. BLOWER 
GUIDE VANES 


\ 
PROTECTIVE 
SCREEN 


EXHAUST HEADER 


SCHEMATIC ARRANGEMENT 
OF ENGINE, AUXILIARY 
BLOWER,AND TURBOCHARGERS 


Fig. 8 Schematic arrangement of engine, auxiliary blower, and turbo- 
chargers 
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maximum specified pressure drops through piping, filters, inter- 
coolers, and so on, whereas the measured pressure drops were 
somewhat smaller. 

In order to check turbocharger combined efficiency during tests, 
measured values may be substituted in equation (8) and solved 
for combined efficiency. However, it should be remembered that 
pressure losses in piping, and the like, are charged against the 
turbochargers in this equation and that the actual efficiency, de- 
pending on where the pressure losses are most severe, may be 
greater or smaller than the measured efficiency. 

The importance of the turbocharger efficiency is best demon- 
strated in Fig. 7 where this efficiency is plotted against fuel con- 
sumption. 
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Fig. 10 Mechanical-efficiency determination by Willan line (fuel-rate 
method) 
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Field-Test Results A schematic arrangement of engine, auxiliary blower, and tur- 
bochargers is shown in Fig. 8. Controlled temperatures such as 
preturbocharger compressor, scavenging header, and preturbine 
exhaust are noted in this arrangement. 

Due to this being the first field application of an engine of this 
model and to be able to compare shop tests with field tests, a 
complete laboratory instrumentation was used during the tests. 


As mentioned earlier, the engine is presently operating as part 
of the municipal power plant in the town of Hudson. Also men- 
tioned earlier, the engine was sold as a dual-fuel engine, but due 
to gas not being available, only operation on oil has been ex- 
perienced at the time of writing. 

The acceptance tests were conducted May 9 through May 12, 
1960, in accordance with DEMA Field Test Code. 
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Fig. 15 Arrangement of turbochargers and air intercooler-engine 
It is the author’s pleasure to note that within the limitation of 


accuracy of this instrumentation the correlation between shop 
tests and field tests was complete. 
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Fig. 14 View of engine from air- 
manifold side 


Fig. 9 shows the engine air-fuel ratio as obtained during these 
tests. The air was measured through an ASME type flow nozzle. 

In Fig. 10 the Willan line method for determination of the 
mechanical efficiency has been used. It is realized that this 
method is not entirely reliable since by this method the friction 
loss is measured at zero load. Thus the value for mechanical 
efficiency as represented in this graph is a little higher than the 
actual value. 

The curves in Fig. 11 show bmep, brake thermal efficiency, in- 
dicated thermal efficiency, and differential thermal efficiency as 
plotted against specific input. These curves represent a particu- 
lar method of evaluating engine ratings by fuel-air ratio effects. 
This method was presented in a previous paper,? and the authors 
claim that, at any value of the differential thermal efficiency 
above 30 per cent, the fuel-air ratio is such that the smoke limit of 
the engine has not yet been reached. In Fig. 11 it wiil be seen 
that the lowest value of the differential thermal efficiency is 40 
per cent, and thus it appears that there is a comfortable margin 
before the smoke limit is reached. It also will be seen that the 
characteristic line is completely straight and that the brake 
thermal-efficiency curve is still on the rise, both of which sub- 
stantiate this theory and indicate that the engine has very good 
load-carrying ability. 

Table 2 shows a comparison of heat-balance data taken on the 
nonsupercharged (74.5 bmep) and the supercharged (105 bmep) 
versions of the engine. 

Fig. 12 contains engine, turbochargers, and auxiliary-blower 
performance curves. 

Figs. 13, 14, and 15 are views of the installation in the municipal 
plant in the town of Hudson, Mass. 


Future Developments 


Based on the encouraging results so far obtained, further in- 


2 R. J. Cramer, Jr., and K. F. Freohlich, ‘‘Rating Engines by Fuel- 
Air Ratio Effects,’’ Proceedings of the ASME Oil and Gas Power 
Division, 1956. 
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creases in load and improvement in fuel economy are contem- 
plated and are under study at the present time. 


DISCUSSION 
T. M. Robie* 


The dotted line curve of Fig. 2 is interesting because of its 
rapid variations. Instrumentation to secure these curves could 
be somewhat difficult. Would the author please tell us in more 
detail the kind of pressure pickups which were used and the 
point at which the pressures were taken. 


3 President, Eibor Services, Inc., Beloit, Wis. Mem. ASME. 
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Author’s Closure 


In answer to Mr. T. M. Robie’s questions regarding instru- 
mentation used in order to obtain the curves on Fig. 2, it is the 
author’s pleasure to furnish the following information. 

The pressure pickup, a 0-100 psi Consolidated Electronics 
Corporation Type 4-313 was located in the connecting rod one 
inch below the wrist pin bushing. The pressure signals were 
fed through a Consolidated Electronics Corporation 1-127 Carrier 
Amplifier to a Minneapolis-Honeywell 906 Visicorder. An 
Electro Products Angular-Synchronizer Model 3850 was used to 
obtain 5 degree crank angle and top dead center. 
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SOREN HANSEN 


Chief Diesel Design Engineer, 
Burmeister & Wain's, Copenhagen, Denmark 


Introduction 


HE diagram treated in the following emanated from 
the work in a committee appointed by CIMAC for drafting 
recommendations for the testing of diesel engines. 

During this work it proved desirable to establish that the out- 
puts stated in contracts and specifications must always be the 
outputs which the engines are to develop in their places of 
service—that is, for example, on the site at high altitudes for 
stationary engines and under tropical conditions in the case of 
engines for ocean-going ships. 

By using the diagram, the question of a general and common 
reference condition, about which it will hardly be possible to reach 
international agreement, may be dropped. In the case of small 
engines manufactured for stock, for which the site is not 
known at the time of the test, the supplier may, for example, 
state the output under the atmospheric conditions in his work- 
shops and moreover refer to the diagram. 

To.a given atmospheric condition corresponds a single reference 
line in the diagram, and the use of tables and graphs prepared for 
each particular atmospheric condition is thus eliminated. 

By means of the diagram it is easy to convert a specified engine 
output into an output under the conditions existing during the 
testing in the engine builder’s workshops, and thus achieve 
that the output stated for the engine is always the output which 
the engine, according to the contract or specification, is to de- 
velop in its place of service. 


Conversion Formulas 


The following considerations apply to normally aspirated 4- 
stroke engines and the 2 and 4-stroke engines in which charging 
and pressure-charging, if any, are effected by means of me- 
chanically driven blowers (piston pumps, Roots, rotary, and cen- 
trifugal blowers) for which the ratio between the speed of the 
charging blower and of the engine is fixed and remains un- 
changed in case of changed atmospheric conditions. 

The following symbols are used: 


N,; = indicated output 

effective output 

absolute ambient temperature 

absolute ambient pressure 

partial pressure of saturated water vapor 
relative humidity 

mechanical efficiency 


For the afore-mentioned engine types, the following relation 
may be assumed to apply between the indicated outputs at two 
different atmospheric conditions, designated by index 1 and 2: 


No (2) = 2 


K 
Na P; OP 


a) 


Contributed by the Oil and Gas Power Division and presented at 
the Winter Annual Meeting, New York, N. Y., November 27—Decem- 
ber 2, 1960, of Tue American Society oF MECHANICAL ENGINEERS. 
Manuscript received at ASME Headquarters, March 29, 1960. 
Paper No. 60 —WA-162. 
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Diagram for Conversion of Diesel Engine 
Qutputs for Changed Atmospheric Conditions 


The vapor pressure is only deducted in case of engines whose 
load limit is determined by excess air, while it is not deducted in 
case of engines whose load limit is determined by thermal load, 
because the specific heat of water vapor is of the same order of 
magnitude as that of the other gases present in the combustion 
space. 

We consider the latter group of engines first, and assuming the 
exponent n = 1 we have: 

T, 


2 
TP, P; (2) 


K, = 


If index 1 applies to a given condition, e.g., the condition of the 
test site, we see (Fig. 1) that, for a given value of 7, K is deter- 
mined by the ordinates to a straight line with point of origin K 
= 0, P, = 0, and with tan Y = 71/(72P:).. The K-value for a 
point corresponding to simultaneous P; and 7: values is thus 
given by the ordinate of the point and is conveniently read off a 
suitable scale. 


Fig. 1 


If it is desired to refer the values 7)’, and P,’ instead of 7, P1, 
tan y is altered in the ratio (7,’/71)(P:/P1’) =). 

If the values 7;/P; have been chosen so that 7,'/P;’ > 71/P, 
we have j > 1, and the alteration of tan Y to 7,'/(T2P,’) is 
equivalent to a reduction of the length of the K, scale in the pro- 
portion 1/y = cos A, i.e., the K,’ values are read on the original 
scale turned an angle Z from the ordinate axis. 

The equation K = (7;/T2)"P2/P, is represented by a similar 
graph. 

If the air humidity is to be considered, equation (1) may also 
be graphically represented, as the abscissa for K,; = 0 is P2; = 
and tan = (7:/T2)"/(P: — (Fig. 2). 

The scale lines are hereinafter called reference lines. 


The Effective Output 


The reduction factor for the effective engine output: 


a= Na = K — a(1 K) 1) (3) 
N el Nmi 

where a is that part of the lost indicated output Na = Na(l — 
Nmi) Which is independent of the atmospheric conditions, and 
(1 — a) that part of the lost indicated output which is propor- 
tional to K. 

Nm = Na/Nx is the mechanical efficiency corresponding to the 
ambient condition under index 1. 
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Assuming a((1/7m) — 1) = C, and substituting K from equa- 
tion (2) we have: 
Ps 


P, (1+ C,) (4) 


a= K(1+() -CQ,= 
which is a co-ordinate system with the P2 values as the abscissas 
are represented by lines with tan Y = (7,(1 + ())/(T:P;) and 
point of origin P; = 0; a = —C, (Fig. 3). 

In the same way as shown in Fig. 2, lines can be drawn which 
correspond to given air humidities, and as stated in the foregoing, 
a change in the condition to 7'’P,’ corresponds to reading the 
a-scale turned an angle \ from the ordinate axis. 


Fig. 2 


It appears immediately from Fig. 3 that a = 1 on any reference 
line corresponds to unchanged C; value, for which reason ,, is 
also unchanged for all @ = 1, but changed for all other conditions 
that are referred to this reference line. 

Fig. 4 shows the diagram in 2 convenient form for practical use. 
It is based on equation (4) with a C; value corresponding to nm 
= 0.8 and a = 0.7. As basic ambient condition has been chosen, 
T; = 273 + 20° C and P; = 760 mm Hg, because higher air 
densities are normally not used for reference. The diagram in- 


4. 


cludes air pressures down to about 400 mm Hg, and besides the P 
scale, the abscissa has a scale for the altitude above sea level. 
The latter scale is not linear. Further, circular arcs are shown for 
@ = 1.0 and 0.5. 

If it is desired to draw a reference line for another condition 
T;', P;', this point is sought out in the diagram and a horizontal 
line is drawn through it to intersect the circle for @ = 1, and the 
distance | is marked = 0.57L, whereupon the new reference line 
is drawn and graduated (see also Fig. 3). 

It would result in a good many lines if the influence of the air 
humidity were to be included in the main diagram, and it has 
therefore been deemed most expedient to include this influence 
as a special correction diagram. 

On the main diagram have been drawn correction lines for 91 
= 0.7, 0.75, 0.85, and 0.90. These correction lines are used in the 
following way: Subtract the correction for air humidity y 
from a. From a point a-y on the (760 mm Hg, 20°C) scale draw 
a horizontal line to the 20 C temperature line. Against the 
atmospheric pressure thus determined, the correction for al- 
teration of the efficiency nm is read. 
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Fig. 6 Consumption correction for nonturbocharged engines 


As appears from Fig. 4, the correction in 7,, for humidity will 
in most instances be so small that it may be ignored, because the 
correction for humidity will rarely exceed 8 per cent. 

The use of the diagram will later on be illustrated by examples. 


The Meshanical Efficiency and 
the Specific Fuel Consumotion 


A change of the engine output from the value at which the 
mechanical efficiency 7,,: is given changes the efficiency in the 
proportion: 


8 N Nak 
which by means of (2) is converted into 


8 K 


Cc 
1+— 
a 


The curve for 8 in relation to @ is in Fig 6 shown for dif- 
ferent values of 7,,,; and witha = 07. 

For engines that are limited by thermal loading and for 
which the outputs are converted according to the corresponding 
formula, the ratio between the charge weights and the outputs 
remains constant, and the consumption per IHPh of the engine 
may also—provided the adjustment of the engine remains un- 
changed—with sufficient accuracy be regarded as unchanged. 

Thus the consumption per BHPh depends only on the change 
in the mechanical efficiency and can be expressed by: 


be = 


For engines that are limited by air excess and for which the 
outputs are converted according to the corresponding formula, 
the ratio between the charge weights and the outputs will increase 
due to the reduction in output for water vapor. This will— 
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provided the adjustment of the engine remains unchanged— 
generally cause a reduction in fuel consumption per IHPh. 

In calculating the modification in fuel consumption it is 
therefore justified to leave the reduction for water vapor out, 
wholly or part!y dependent on the fuel consumption characteristic 
of the engine. . 


Examples 


(1) Draw reference lines for a test bed for which the average 
ambient conditions are 745 mm Hg, 20° C, @ = 0.6. 


Solution 

From the point (1), Fig. 4, which represents the mentioned 
test-bed conditions, a horizontal line is drawn to the circle a = 1. 
L is measured to 47 mm, | = 0.57L = 27 is marked and the 
reference line is drawn. A horizoutal line through 20° C and 
¢ = 0.6 is drawn in the correction diagram for humidity. 

(2) Under the afore-mentioned test-bed conditions, the follow- 
ing has been determined: 


Continuous output 
Mechanical efficiency 0.85 
165 g/BHPh 


excess of air 


The engine site has the following conditions: 2200 m above 
sea level, 30° C, and @ = 0.85. 

It is desired to determine output, mechanical efficiency, and 
fuel consumption. 


Solution 
Horizontal line from point 2 (Fig. 4) gives a = 0.71. 
Correction for 
Correction for air humidity 


a, = 0.71 + 0.01 — 0.02 = 0.70 
Ne = 1000-0.70 = 700 BHP 
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85 
= —— = 0.81 
1.045 


B = 1.045 Nm2 


be = 1.045-165 = 172 g/BHPh 


(3) Under test-bed conditions like those in example 1, the 
following has been determined: 


The engine site has the following conditions: 
O m above sea level; 40° C, and @ = 1.00 


It is desired to determine the output, mechanical efficiency, and 
fuel consumption, with correction in efficiency due to the air 
humidity correction. 
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Solution 
Horizontal line from point 760 mm Hg, 40° C gives a = 0.95. 


Correction for air humidity 
Correction for efficiency 


as 0.95 — 0.05 = 0.90, to which corresponds P, = 695 mm Hg. 
Qrotar = 0.95 — 0.05 — 0.01 = 0.89 
Ne = 20°0°89 = 17.8 ~ 18 BHP 


.70 
1.03 = 0.68 
be = 1.03-200 = 206 g/BHPh 


If the content of water vapor is not to be considered, 8, is 
determined from an @ value of 0.95 — 0.01 = 0.94. 


Be = 1.015; be = 200 X 1.015 = 203 g/BHPh 


8B = 1.03 


= 


aprit 1961 / 207 


a 
—0.05 
Continuous output. 20 BHP 
Mechanical efficiency.................... 0.70 
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